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THE LATTICE SPACINGS OF NICKEL SOLID SOLUTIONS! 


W. B. PEARSON AND L. T. THOMPSON? 


ABSTRACT 
New lattice spacing measurements in terminal nickel solid solutions with 
Cr, Mn, Co, Zn, Ga, Ge, In, and Sn solutes are reported and discussed, and 
the relationship of lattice distortion and the depression of the solidus tem- 
perature is examined in solid solutions of the First Long Period elements 
in nickel. 


INTRODUCTION 


Measurements, most of which have been carried out in recent years, of the 
lattice spacings of nickel solid solutions containing as solutes the transition 
metal elements of the First Long Period* provide one of the most systematic 
sources of information on lattice spacing variation in transition metal alloys. 
These measurements are particularly interesting because of the definite 


kinks found in the lattice spacing — composition curves of vanadium, chromium, 
manganese, and iron dissolved in nickel. Taylor and Floyd (1952) associated 
the change of slope in the curve of lattice spacing vs. composition in the Ni-Cr 
solid solution with the occurrence of ordering in the Ni;Cr region. The curve 
(Fig. 1) obtained by a combination of their measurements with those of 
Pearson, carried out in 1950, but reported for the first time here, shows that 
the bend occurs at about 20 at. % Cr. It does not, therefore, appear to be due 
to the superlattice ordering at Ni;Cr which takes place at temperatures below 
570° C. Similarly in the lattice spacings of the Ni—V solid solution (Fig. 1), 
although the bend occurs in the region of the compound Ni;V, the slopes 
of the lattice spacing curves at higher and lower compositions are considerably 
different and nearly linear suggesting that the effect is independent of the 
ordering process which takes place on formation of Ni;V. The large change 
of slope found in the iron solid solution at about 60 at. % Fe (Fig. 3) is asso- 
ciated with the formation of the alloy “Invar.’’ Bradley et a/. (1937) also 


found indications of a slight bend in the lattice spacing curve at about 22 at. 


1Manuscript received November 5, 1956. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
Canada 

Issued as N.R.C. No. 4245. 

2Now with Defence Research Board, Ottawa, Canada. 

*Notably, for instance, NiTi: Taylor and Floyd (1952); Ni-V: Pearson and Hume-Rothery 
(1952); Ni-Cr: Taylor and Floyd (1952) and Pearson (unpublished, reported here for the 
first time); Ni-Mn: Tyler and Hume-Rothery and Pearson (unpublished, reported here for 
the first time); Ni-Fe: Owen and Sully (1941), Bradley et al. (1937), Jette and Foote (1936); 
Ni-Co: Taylor (1950); Ni-Cu: Owen and Pickup (1934), Coles (1956). 
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% iron in the ‘permalloy’ region; however, this bend was not apparent 
in the work of Jette and Foote (1936) or of Owen and Sully (1941). 

In contrast to the work on transition metal solutes in Ni, very little has 
been done recently on the behavior of the B group elements dissolved in Ni. 
On the other hand such an examination may well be important, for investiga- 
tion of these solutes dissolved in Cu, Ag, and Au has been largely responsible 
for an understanding of the distortion caused by solutes of different valencies 
in solid solutions (cf. Raynor 1949; also Owen 1947). Previously measure- 
ments of lattice spacings of the Ni-Zn solid solution have been made by 
Tamara and Osawa (1935), but the lattice spacing of a = 3.503 kX given 
for pure nickel (accepted value a ~ 3.5167 kX) indicates that the results 
are not too reliable. Measurements of several alloys in the Ni-Ga solid solution 
were made by Hellner (1950), but he failed to notice the formation of a Ni;Ga 
superlattice in the solid solution region (Pearson 1954). Weibke and Pleger 
(1937) measured the lattice spacings of alloys in the Ni-In solid solution, 
but they found a rather low lattice spacing for pure nickel (a = 3.514 kX). 
Measurements of the lattice spacings of the Ni-Sn solid solution made by 
Mikulas et a/. (1937) and by Jette and Foote (1935) are in rather good agree- 
ment with each other. Finally lattice spacings of the Ni-Sb solid solution 
have been measured by Osawa and Shibata (1940). 

We have thus undertaken anew systematic measurements on the lattice 
spacings of the B group solutes in nickel with three main interests in view: 
(i) to gain further understanding of the electronic behavior of Ni as solvent 
in relation to that of the succeeding monovalent metals Cu, Ag, and Au; 
(ii) to see if any further systematic information can be gained about the 
bends which occur in the lattice spacing -composition curves of certain Ni 
solid solutions (somewhat similar bends have also been found in some Cr and 
some Fe solid solutions); and finally (iii) to examine the relationship of 
macroscopic lattice distortion (as determined by lattice spacing variation) 
and the depression of the solidus temperature caused by a solute, for we 
expect that the solidus depression should show a certain dependence on the 
distortion of the crystal lattice. Nickel, as we shall show, is a rather interest- 
ing solvent from this point of view, as all of the transition metal and B group 
metal solutes are (nominally at least) of larger size than the solvent atom. 


EXPERIMENTAL METHODS 

The alloys have been prepared from Mond nickel shot (typical analysis: 
99.95% nickel, 0.02% carbon, 0.02% iron) kindly supplied by the Inter- 
national Nickel Co. and electrolytic hydrogen-treated Cr (99.9%), Mn 
(99.9%), and Co (99.9%) from Johnson, Matthey & Co. Ltd., 99.999% Zn 
from the New Jersey Refining Co. Ltd., 99.9% Ga, from Johnson, Matthey & 
Co. Ltd. or Mackay Inc., 99.999% Ge from various sources, ‘spec pure’ 
Cd and Sn from Johnson, Matthey & Co. Ltd., and 99.9% In from the same 
firm or from the Indium Corporation of America. With the exception of the 
Ni-Zn and Ni Cd alloys they were all melted in vacuo in an induction furnace 
in pure alumina crucibles and then homogenized for several days at 1100° C. in 








PEARSON AND THOMPSON: LATTICE SPACINGS 351 


alumina collars sealed in evacuated silica tubes. Filings were prepared and 
strain-relieved, generally at 800° C. The composition of the actual X-ray 
powders was determined by chemical analyses of both components and also 
of iron. The first series of powders prepared was found to contain generally 
as much as 0.5% Fe which was picked up during the filing process. The 
measurements of these samples were therefore discarded, and subsequently 
the iron was carefully removed from the powders with a magnet by heating 
them in vacuo to a temperature above their Curie point so that the iron 
remained ferromagnetic but not the alloys. The iron content of most powders 
examined then generally lay in the range 0.02 to 0.25 wt. %. 

The Ni-Zn and Ni-Cd alloys were prepared by interdiffusion by heating 
together mixtures of filings of the pure metals in sealed silica tubes for a 
sufficient time to attain equilibrium. 

Lattice spacings were determined from powder photographs taken in a 19 em. 
diam. Unicam camera in a thermostatically controlled room. CuK, radiation 
was generally used throughout the work. Scattered radiation from Ni-Co 
alloys was absorbed with aluminum foil placed between the specimen and the 
X-ray film. Lattice spacings were calculated using the extrapolation method 
of Nelson and Riley. The X-ray measurements on any specimen were repro- 
ducible to 0.0001 kX, but factors involving the alloys (impurities, chemical 
analyses, etc.) made the general reproducibility of the measurements from 
different alloys somewhat less, particularly it seems in the Ni-Ge svstem. 

EXPERIMENTAL RESULTS 

The lattice spacing of Mond nickel shot (a = 3.5169 kX at 18°C.) was 
found to be higher than that of Johnson and Matthey spectroscopically 
pure Ni (a = 3.5165 kX at 18° C.), which was absolutely cobalt-free and 
contained only trivial amounts of other metallic impurities. This agrees with 
the observations of Coles (1956). 

The lattice spacing variations in the systems Ni-Cr, Ni-Mn, Ni-Co, 
Ni-Zn, Ni-Ga, Ni-Ge, Ni-In, and Ni-Sn, together with measurements by 
others on Ni-Ti, Ni-V, Ni-Fe, and Ni-Cu solid solutions, are shown in a 
series of diagrams (Figs. 1-5). We should notice the sudden change in slope 
of the lattice spacing curves in the Ni-V, Ni-Cr, Ni-Mn, and Ni-Fe solid 
solutions at about 24, 20, 60, and 60 at. % solute respectively. Smaller changes 
of slope may also exist in the Ni-Zn, Ni-Ga, Ni-Ge, and Ni-In systems. 
Measurements on the Ni-Sb solid solution by Osawa and Shibata (1940, 
1942) apparently indicate a steplike increase in the lattice spacing — composi- 
tion curve at 6 wt. % Sb which they attributed to the occurrence of a super- 
lattice. 

It is possible to compare the present work with that of other investigators 
Taylor and Floyd (1952) 


in several of the systems. The agreement with 
on the lattice spacings of the Ni-Cr solid solution is excellent, and geod 
agreement is also found with Taylor (1950) in the Ni-Co system and with 
Weibke and Pleger (1937) in the Ni-In system. Because of the Ni;Ga super- 
lattice only one of Hellner’s measurements on the Ni-Ga solid solution can 
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Fic. 1. Lattice spacings of terminal nickel solid solutions. 
O Ni-Ti Taylor and Floyd (1952) A Ni-Cr Taylor and Floyd (1952) 
O Ni-V_ Pearson and Hume-Rothery (1952) V Ni-Cr Pearson 


Fic. 2. Lattice spacings of the Ni—Co solid solution. 


A Taylor (1950) O Present work 
Fic. 3. Lattice spacings of terminal nickel solid solutions. 
A Ni-Mn_ Pearson O Ni-Cu_ Coles (1956) 
V Ni-Mn_ Tyler O Ni-Cu Owen and Pickup (1934) 


xX Ni-Fe Owen and Sully (1941) 


be compared with our results and the agreement is reasonable. In the Ni-Sn 
system our results lie slightly higher than the measurements of Mikulas ef 
al. (1937) and of Jette and Foote (1935), which agree with each other. How- 
ever, it appears that only in the present investigation were the actual X-ray 
powders chemically analyzed, and as the analytical totals deviate from 
100°, by no more than 0.09%, we have confidence in these new measurements. 

The solubility of Cd in Ni appears to be rather small at 800° C. We pre- 
pared a number of alloys by the method which was successful for the Ni-Zn 
solid solutions; however, as the alloys regularly showed a loss of weight 
during sintering and the X-ray films frequently showed extra reflections 
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due to some impurity phase we have rejected the results. The lattice spacing 
of heterogeneous alloys quenched from 800° C. was found consistently to be 
a= 35180;KX. 
DISCUSSION OF RESULTS 

(a) Lattice Distortion 

Fig. 6 shows the total initial lattice distortion Aa, per atomic per cent 
solute for the elements of the First Long Period dissolved in nickel, together 
with the sizes of the solute atoms relative to nickel.* It is perhaps remarkable 


tT Vv Cr Mn Fe Co WNi Cu Zn Goa _ Ge 


in 107! K units 


+40 


+30 4a 


in 107'X units 
+20 





Ti Vv Cr Mn Fe Co Ni Cu Zn Go Ge 
SOLUTE 


Fic. 6. (a) Difference in sizes of solute atoms relative to Ni (see text), plotted as the 
increase in lattice spacing expected per atomic per cent solute dissolved in nickel. 

(6) Total observed initial lattice distortion per atomic per cent solute in terminal nickel 
solid solutions. 


(c) Initial solidus depression in ° C. per atomic per cent solute in terminal nickel solid solutions. 
how well the relative variation in lattice distortion (except with Ge as solute) 
follows the relative variation of solute atom size. In the Ni-Ge solid solution 
it is apparent that specific electronic or valence effects contribute to the 
observed lattice distortion. 








*The relative atom sizes for a hypothetical face-centered cubic lattice have been calculated 
for the body-centered cubic metals from the interatomic distances, which were increased by 
3% to allow for increase in coordination Following the lead of Raynor (1949) the interatomic 
distances of Zn and Ge were taken as the closest distances of approach of the atoms in the 
pure metals and the same procedure was followed for hexagonal Ti. The atomic diameter of 
Ga was calculated from the observed lattice distortion in Cu, Ag, and Au as solvents using 
Raynor's relationships for the data of the Group II and IVB solutes in the same solvents. 
The size of the Mn atom was obtained by taking a face-centered cubic unit cell which had 
the same volume as the face-centered tetragonal y-Mn unit cell at room temperature. This 
is the most satisfactory method of estimating its size as Christian and Basinski have shown 
that cubic y-Mn transforms to tetragonal 7-Mn without change of volume. 
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At present the only apparently successful quantitative analysis of lattice 
spacing variation in terminal solid solutions is that given by Raynor (1949) 
for the B group solutes dissolved in Cu, Ag, and Au. In this analysis the 
total observed distortion was divided into parts due to a size effect and to a 
valence effect and, for solvent and solutes which lay in the same period, 
the valence distortion was found to be proportional to the square of the 
valence difference of solvent and solute. It can be seen from Fig. 6 (a) and 
(6) that a similar analysis of solid solutions in nickel (taking the same sizes 
as were used by Raynor for Zn, Ga, and Ge atoms) would not lead to this 
same dependence of the electronic component of lattice distortion on the 
square of the valence difference of solvent and solute. 

With our present picture of the nature of a solute ion in solid solution 
(derived largely from an understanding of electrical measurements) it is 
perhaps hard to understand the nature of a size effect which can be separated 
entirely from an electronic effect in the general way found by Raynor. It is 
therefore of particular interest to note that the factors which control the 
sizes of the transition metal atoms in the pure state also appear to play a 
large part in determining their apparent sizes when they are dissolved in 
nickel. The number of s-electrons does not vary greatly throughout the 
series of transition metals, but the d-electron content varies more or less 
regularly from one element to the next and, although it is known that these 
electrons contribute to the cohesion of the lattice, the cohesion is certainly 
not a linear function of the d-electron content of the metal (cf. Hume-Rothery, 
Irving, and Williams 1951). The correspondence between the observed dis- 
tortion in nickel and the relative solute sizes obtained from the pure metals 
is therefore perhaps rather closer than might have been expected, since the 
electronic state of the solutes in nickel is presumably different to that in the 
pure metals. However, judging by known data, the same observation appears 
to be generally true for solid solutions of most of the transition metals in 
each other, although there are some obvious exceptions to this rule. In an 
attempt to gain more knowledge of the relative electronic states of these solutes 
in nickel, measurements of the resistivity of these alloys are being undertaken. 
(6) Bends in Lattice Spacing Curves 

In the terminal iron solid solutions of the Fe-Co and Fe-Cr systems sharp 
bends are found in the lattice spacing - composition curves at the same com- 
positions as bends also occur in the curves of Curie temperature and saturation 
magnetic moment as functions of composition. This correspondence does 
not appear to hold in nickel solid solutions and this is perhaps not surprising 
if the magnetic interactions are weaker. In the Ni-Mn solid solution, for 
instance, the saturation magnetic moment first increases for additions of up 
to 8 or 10 at. % Mn after which it decreases very rapidly (Sadron 1932); 
however, the few lattice spacing measurements which we have obtained 
in this region certainly indicate that the lattice spacing changes linearly as a 
function of composition up to about 50 at. % Mn. 

Both Owen and Yates (1936) and Esser ef a/. (1938) have shown that the 
lattice spacing of pure nickel as a function of temperature passes through a 
small anomaly in the region of the Curie point, and it would therefore be 
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expected that the lattice spacing-composition curve of alloys should show 
an anomaly at the composition at which the Curie point passes through room 
temperature. The lattice spacing — composition curves shown in Figs. 1 and 
3 indicate that this anomaly must be rather small and not bigger than the 
general scatter of the experimental measurements, for it has not been possible 
to observe it, except perhaps in the Ni-Cu system where Coles’ (1956) measure- 
ments showed a change in the rate of deviation from Vegard’s law in this 
region. 

It does not seem possible to find any explanation of the large bends in the 
lattice spacing — composition curve of the Ni-V, Ni-Cr, and Ni-Mn systems 
in terms of Brillouin zone overlaps, and even if the bend in the Ni-Fe solid 
solution were in part due to magnetic interactions which changed rapidly 
with composition in that region, there seems to be no clear evidence that 
this should also be the case in the other solid solutions. It is, however, 
probably significant that in the Ni-V and Ni-Cr solid solutions the rate of 
lattice distortion beyond the bend increases, whereas, when Mn or Fe is 
solute there is a decrease in the rate of distortion of the lattice as a function 
of composition beyond the bend in the curve of lattice spacing vs. composition. 
This change thus follows changes in the sequence of other physical properties 
which occur between Cr and Mn and which probably result from a change 
in the distribution of electrons in the d-band as the second half of the 


d-band starts to be filled. 


(c) Correlation of Lattice Distortion and Solidus Depression 

The solidus depression, A7m,,* in a terminal solid solution may be expected 
to show some correlation with the total lattice distortion, Aa, as determined 
from lattice spacing measurements, for the energy stored in the lattice as the 
result of a positive distortion (i.e. an increase in lattice spacing) should lessen 
the total thermal energy required to disrupt the lattice in the melting process. 
In other words we would expect the solidus temperature to depend very 
largely on the properties of the crystalline lattice of the solid phase. Indeed, 
when we compare (Fig. 6 (b) and (c)) the total lattice distortion and solidus 
depression for nickel solid solutions we find an excellent qualitative correlation. 
Such a correlation is lost if we compare estimated size or electronic com- 
ponents of the distortion with the solidus depression, showing as expected 
that it is only for the fofa/ distortion that the relationship exists. The relation- 
ship is perhaps more remarkable when it is remembered that the lattice spacing 
measurements were made at room temperature and there is no a priori reason 
to expect that the same relative distortions should be found in the various 
solid solutions in the vicinity of the melting point of the alloys. 

When we compare similar data for solid solutions in Cu and in Au we find 
that, although many of the features are similar, a general correlation of 
solidus depression and lattice distortion does not exist. This may be due to the 
relative lattice distortions near the melting point being different in these 
solvents from those observed at room temperature, or it may in part be due 


*°C, per atomic per cent solute. 
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to the fact that some solutes cause positive distortions and others negative 
distortions of the Cu and Au lattices, while nickel is unique as a solvent in 
having a smaller diameter than all of the solute atoms and each solute causes 
a positive distortion of the lattice. If a negative distortion of the solvent 
lattice (i.e. a contraction) occurs, it seems unlikely that this should necessarily 
lead to an increase in the solidus temperature for, although the cohesion of the 
solvent atoms may increase, any strain energy in the lattice would still be 
expected to contribute towards a lowering of the melting point. 

A thorough discussion of ‘“‘Melting and Crystal Structure’’ was given by 


Ubbelohde (1950). 
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LATTICE DISTORTION IN TERMINAL SOLID SOLUTIONS OF THE 
FIRST LONG PERIOD TRANSITION METALS AND COPPER! 


W. B. PEARSON 


ABSTRACT 


Collected data of the gross lattice distortion in solid solutions of the First Long 
Period elements in each other and in copper are given together with similar 
data for B group solutes in Cu, Ag, and Au. It is shown that there is no correlation 
between the size of the lattice distortion and the solid solubility of the First Long 
Period transition metals in copper. 


The whole of the X-ray work on metals and alloys has recently been reviewed 
and a compilation of the most reliable structural and lattice spacing data is 
in the course of publication (Pearson 1957a). Using these collected data we have 
now calculated the initial lattice distortion (expressed as change of lattice 
spacing per atomic per cent solute) in terminal solid solutions of the transition 
metals of the First Long Period in each other and with the B group elements 
as solutes. This information is of particular interest in relation to the other 
physical properties of transition metal alloys such as electrical resistivity and 
thermoelectricity and to recent considerations of their electronic nature 
(e.g. Friedel 1956). 

Information about the local, microscopic distortion caused by a_ sub- 
stitutional solute in a solvent lattice can be gained from examination of the 
diffuse X-ray scattering accompanying the spectral reflections (cf. Averbach 
et al. 1953, 1954; Warren and Averbach 1953); on the other hand lattice 
spacing variation gives a direct measurement of the gross or macroscopic 
distortion of the solvent lattice due to the solute. In Fig. 1 we show all the 
available data for the initial* macroscopic lattice distortion per atomic per 
cent of solute in the terminal solid solutions of the transition metals of the 
First Long Period plotted as a function of the atomic number of the solute, 
and in Fig. 2 we show the same information as a function of the atomic number 
of the solvent. In each case we have shown the actual distortion (Aa) rather 
than the proportional distortion (Aa/a) as this can more conveniently be 
compared with the relative sizes of the solvent or solute atoms, which are 
also plotted in Figs. 1 and 2. The method used to obtain the sizes of the 
atoms and to compare them in terms of a (hypothetical) face-centered cubic 
lattice has been described by Pearson (1957d). 

The factors causing lattice distortion are still not too well understood and 
in the only successful quantitative interpretation, which is due to Raynor 
(1949), the total distortion caused by the B group solutes in solid solution 
'Manuscript received November 5, 1956. 

Contribution from the Division of Pure Physics, National Research Council, Ottawa, 


Canada. 
Issued as N.R.C. No. 4241. 
*Details of references from which the data were collected can be found in Pearson (1957a). 


We have taken the best average of the initial slope of the lattice spacing - composition 
curves. 
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Fic. 1. Initial lattice distortion per atomic per cent solute in terminal solid solutions of the 
transition metals and copper plotted for the various solvents as a function of the atomic 
number of the solute, together with the difference of the solute sizes relative to nickel, plotted 
as the increase in lattice spacing expected per atomic per cent dissolved in nickel. 

? V-Fe: the two points mark the probable limits of the somewhat uncertain size of the 
lattice distortion. 

? Cr-Ti: the value of the initial lattice distortion is not too certain. 

? Cr-Co: obtained from very old measurements made in the region of 25 wt. % Co. 

? Cu-Ti: the results of two recent investigations differ rather widely. The value plotted 
is a mean value weighted slightly in favor of the work which appears to be the most reliable. 


in Cu, Ag, or Au was separated into components due to the relative sizes 
and relative valencies of the solutes and solvent. The occurrence of other 
factors causing lattice distortion has also been recognized (Raynor 1949; 
Axon and Hume-Rothery 1948) but it is not possible to calculate them a 
priori. With our present physical picture of a screened solute ion in these 
solid solutions—a concept largely derived from an understanding of electrical 
measurements—it is perhaps hard to see how a size effect can exist indepen- 
dent of the electronic or valence effect as found by Raynor. It is therefore 
interesting to note in Figs. 1 and 2 that the lattice distortion in transition 
metal alloys generally follows the relative variation in the solute or solvent 
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Fic. 2. Initial lattice distortion per atomic per cent solute in terminal solid solutions of the 
transition metals and copper plotted for the various solutes as a function of the atomic number 
of the solvent, together with the difference of the solvent sizes relative to nickel, plotted as the 
increase in lattice spacing expected per atomic per cent dissolved in nickel. For comments 
on ? V-Fe, ? Cr-Ti, ? Cr-Co, and ? Cu-Ti see caption to Fig. 1. 


atom sizes.* This is certainly not the case in the solid solutions of the B 
group solutes in Cu, Ag, and Au (Fig. 3), for the sizes of the solutes of the 
First and Second Long Periods decrease in the elements following Zn and Cd 
respectively, while the distortion which they cause in Cu, Ag, and Au con- 
tinues to increase. However, it can also be seen from Fig. 3 that although the 
solute size is not the major factor in causing the lattice distortion in these 
solid solutions, it does exert some influence. 

Although we do not intend in the present note to discuss these lattice 
distortion data in any detail, we do, however, examine them in relation to 
one of the most interesting unsolved problems of metallurgical chemistry; 


*There are of course notable exceptions to this, such as some solid solutions in which Mn 
is the solute. 
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Fic. 3. Initial lattice distortion per atomic per cent solute in solid solutions in Cu, Ag,”or 
Au plotted as a function of the atomic numbers of the B group solutes of the First and Second 
Long Periods. ? indicates some uncertainty in the observed lattice distortion values. 


that is to the relative variation of solid solubility of the transition metals 
in copper. In Table I we compare the total measured lattice distortion and 





TABLE I 
Aa, Aa,, Aa,, 
total measured pure metal apparent ‘elec- 
Solute Solid solubility distortion per at. sizes relative to tronic’ distortion 
% solute in 10-' Cu.expected of Cu lattice per 
X units distortion per at. % solute in 
at. % solute 107? X= units 
in 107! X units 
Ni 100% —10 ~* wel 
Co Very small at low tem- ? —- 7 - 
perature 
Fe Very small at low tem- + 5 0 +5 
perature, ~5% at 
1090° C. 
Mn 100% +37 + 8 +29 
Cr Very small +28 + 2 +26 
V Probably very small ? +21 
Ti About 6 at. % ~+16 +48 ~ —32? 
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also the apparent ‘electronic’ distortion* of the copper lattice with the solid 
solubility limits of the various transition elements of the First Long Period. 
This shows, surprisingly perhaps, that there is no correlation between a small 
solid solubility and a large lattice distortion, and the factors which restrict 
the solid solubility of Co, Fe, Cr, and V in Cu do not at the same time lead 
to a large distortion of the copper lattice, so that the behavior of the transition 
metal solutes is again seen to be different to that of the B group solutes. 
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THE VELOCITY DEPENDENCE OF CERENKOV 
RADIATION FROM u-MESONS! 


C. H. MILLAR AND E. P. HINcKs 


ABSTRACT 


The intensity—velocity relationship of the Cerenkov light emitted by u-mesons 
passing through Plexiglas has been studied for the velocity range 0.83c to 
0.99970c (84 Mev. to 4.2 Bev.). No departure of the observed from the theoretical 
intensity was found in the high velocity region, in disagreement with the results 
of Bassi, Bianchi, and Manduchi (1952). In the low velocity region slightly less 
radiation was observed than expected, but this is attributed to an optical effect 
in the apparatus rather than to a departure from theory of the Cerenkov light 


emission. 
INTRODUCTION 
The classical theory (Frank and Tamm 1937) of Cerenkov radiation by a 
moving, charged particle gives for the energy radiated in a distance x in a 
medium of refractive index n 


(1) W = (e’x/c’) J, , (1-1 /Bn') » dy, 


where 8 (=v/c) and e are the velocity and charge respectively of the moving 
particle, v is the frequency of the emitted light, and c is the velocity of light 
in vacuo. Ginsburg (1940) has obtained by a quantum calculation a formula 
closely coincident with the above. From (1) it will be seen that after the 
particle velocity exceeds a threshold value (Btnresnoa = 1/n) W at first increases 
rapidly with particle velocity, but as 6 approaches unity, W asymptotically 
approaches a constant maximum or ‘plateau’ value. 

Experimental evidence has been reported that the velocity dependence of 
the intensity of Cerenkov light from mesons passing through Plexiglas deviates 
significantly from that given in equation (1). Bassi, Bianchi, and Manduchi 
(1952) reported a plateau intensity of Cerenkov light from fast cosmic-ray 
u-mesons whose relative value was about 25% too high when compared with 
the intensity of light from mesons in a velocity range corresponding to a 
steeply rising portion of the theoretical curve. Winckler, Mitchell, Anderson, 
and Peterson (1955), on the other hand, using artificially produced 2-mesons 
observed good agreement with the steeply rising part of the expected intensity 
curve, but obtained a slight indication that the plateau commenced at a lower 
velocity than theoretically predicted. Since some Cerenkov counters in use 
at Chalk River depend on this intensity—velocity relationship to measure 
particle energies, we have carried out an experiment using cosmic-ray u-mesons 
to investigate these reported anomalies. 


‘Manuscript received December 10, 1956. 

Contribution from Physics Division, Atomic Energy of Canada Limited, Chalk River, 
Ontario. 

Issued as A.E.C.L. No. 410. 
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EXPERIMENTAL PROCEDURE AND OBSERVATIONS 

The general method of this experiment was to select a vertically-directed 
u-meson beam by means of a coincidence telescope and to study the light 
output from a Cerenkov counter placed in this beam. The energies of the 
mesons were determined by their ranges after leaving the Cerenkov counter. 

Since the experiment was carried out at Chalk River (altitude 475 ft.), and 
since a considerable thickness of absorber was always present either within 
or above the telescope, about 99% of the particles selected by the telescope 
were yu-mesons. The contamination of the beam by the remainder of the 
particles was estimated to have a negligible effect on the experimental results. 
'n addition, an auxiliary experiment was performed in which showers were 
detected by recording coincident pulses in trays of Geiger counters placed 
nearby but out of the telescope beam. The results indicated that in the region 
of the most probable Cerenkov light output by mesons the effect of showers 
was also negligible. 

The Cerenkov counter consisted of a 13 in. X13 in. x5 in. block of Type II 
Plexiglas viewed by two RCA No. 5819 photomultiplier tubes centrally 
positioned on opposite edges of the block. The Plexiglas was surrounded by 
MgO reflecting powder and a thin aluminum light shield. 

It should be noted that the photomultiplier tubes used had an ‘S-4’ spectral 
sensitivity so that they responded only to visible and ultraviolet light in a 
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Fic. 1. Arrangement of counters and absorber for Part I of experiment. 
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restricted wavelength band. Thus the term “Cerenkov light” as used in this 
work refers to light in the region from 3200 to 6000 A. 

The experiment was carried out in two phases. The first phase was designed 
to look for any anomalous intensity variation of Cerenkov radiation that 
might occur in the high velocity or plateau region. The arrangement of appara- 
tus for this part of the experiment is shown in Fig. 1, and is described below. 

Vertically directed charged particles were selected by a triple coincidence 
between the two top liquid scintillation counters (Scl, Sc2) and the Cerenkov 
counter (Cv). For a selected event the pulse height from the Cerenkov counter 
was measured by a 30-channel kicksorter and automatically recorded on an 
IBM punch card. On the same card was recorded which, if any, of the three 
lower scintillation counters (Sc3, Sc4, Sc5) had detected a particle in coinci- 
dence with the selected event. By sorting the cards according to this latter 
coincidence information it was possible to obtain pulse height distributions 
from the Cerenkov counter measured simultaneously for particles of four ad- 
jacent range bands. 
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Fic. 2. Pulse height distributions from Cerenkov counter measured in Part I of experiment. 
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A u-meson required a minimum kinetic energy of 1.5 Bev. when passing the 
central plane of the Cerenkov detector in order to penetrate the absorber and 
produce a signal in Sc2. The corresponding energies to reach counters Sc3, 
Sc4, and Sc5 were 2.4, 3.3, and 4.2 Bev. respectively. Thus Cerenkov light 
distributions were obtained for yu-mesons in three 0.9 Bev. energy bands 
between the above limits, and also in a fourth band comprising all energies 
above 4.2 Bev. These four distributions (V-VIII) are shown in Fig. 2. (The 
bands are numbered in order of increasing energy over both parts of the experi- 
ment.) The experimental points have been normalized so that each distribution 
contains the same total number of counts and one may be compared with the 
the other through the solid curves, which are identical in each plot. All four 
experimental distributions are of the same shape, and in particular have the 
same most probable light intensity as far as can be determined by visual 
inspection. 

Assuming that the shapes of the four distributions are, in fact, identical, 
a quantitative estimate of the degree of equality of the peak positions may be 
made by comparing average values of light intensity calculated for the region 
covered by the kicksorter. These average values—denoted by A—are thus 
for a restricted range from zero to about 24 times the most probable pulse 
height. The comparison of individual values of A for each curve with the A 
for all data combined is shown in Table I and it may be seen that within the 
indicated standard errors all the values are identical. We conclude that the 
most probable values are constant to the same accuracy (within +~1%) as 
are these A values, and this conclusion is strengthened by a comparison of the 
r.m.s. deviations from A for each of the four distributions. These values, 
also listed in Table I, themselves agree to within +~1% and prove that the 
calculated values of A are not influenced significantly by any changes in 
shape of the wings of the curves. 





TABLE I 
SUMMARY OF NUMERICAL RESULTS (PART I) 











‘Average’ Normalized r.m.s. % pulses 
Energy pulse ‘average’ deviation greater 
Band range Velocity height pulse from than 
number (Bev.) range (A) height ‘average’ 34.5 
Vv 1.5-2.4 0.9978-0.9991c 16.77+0.12 0.998+0.006 0.281 1.4+0.3 
VI 2.4-3.3 0.9991-0.99952c 16.72+0.16 0.995+0.009 0.276 2.3240.5 
VII 3.3-4.2 0.99952-0.99970c 16.98+0.22 1.010+0.012 0.275 2.9+0.8 
VIII >4.2 >0.99970c 16.83+40.11 1.001+0.005 0.275 3.3+0.4 
V-VIII > 


1.5 >0.9978c 16.81 +0.07 1.000 (0.276) 2.5+0.3) 


It appears then that the most probable Cerenkov light output remains 
constant to within 1% for y-meson energies from 1.5 Bev. up to above 4.2 
Bev. This conclusion does not, however, apply to the true mean Cerenkov 
intensity since a few large pulses formed a tail to each distribution which was 
not measured. These pulses, whose size fell above the range of the kicksorter, 
were of course not included in the calculations of the A values used above to 
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compare the peak positions. Their number did, in fact, show a significant 
increase with meson energy as shown in the last column of Table I. This 
suggests that the mean Cerenkov intensity may be increasing, perhaps as a 
result of secondary Cerenkov radiation from knock-on electrons. However, 
lack of knowledge of the size distribution of these large pulses as well as the 
lack of adequate protection against showers in our apparatus prevents any 
quantitative determination of the variation of mean Cerenkov light output 
with meson energy. 

In the second phase of the experiment the energy bands studied were lowered 
by reducing the amount of lead absorber, the arrangement of apparatus being 
as shown in Fig. 3. Here the initial particle selection was by a triple coincidence 
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Fic. 3. Arrangement of counters and absorber for Part II of experiment. 


between liquid scintillation counters Scl and Sc2, and the coincidence output 
of the crossed Geiger counter trays Gl and G2. The absorber thicknesses 
which the mesons had to penetrate to reach scintillation counters Sc2, Sc3, 
and Sc4 were adjusted to correspond roughly to those used by Bassi ef al. 
(1952). Additional absorber and the scintillation counter Sc5 were included 
in order to overlap the energy bands studied in the first part of this experiment 
and thus to enable normalization of the one set of results to the other. Because 
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of the small thickness of absorber required for the lowest energy group, 
an additional 86.5 g./cm.? of lead was placed above the concrete roof to heip 
filter out the soft component of the cosmic radiation. The boundary values of 
the meson energy bands in this second experiment were 0.084, 0.27, 1.02, and 
1.85 Bev., and four simultaneous pulse height distributions were obtained as 
previously. The results are shown in Fig. 4, where again the four curves have 


been normalized to equal areas. 
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Fic. 4. Pulse height distributions from Cerenkov counter measured in Part II of experiment. 
Note that the arbitrary pulse height scale differs from that in Fig. 2. 


In this instance there are obvious peak shifts as well as shape changes among 
the various curves so that we are not justified in using average values to 
compare peak positions as was done previously. Instead the peak positions 
and their uncertainties have been estimated from individual curves fitted to 
the sets of experimental points. 

The results of the second part of the experiment are summarized in Table 
Il. The most probable pulse height values for each energy group are normalized 
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to make that for the >1.85 Bev. group equal to unity in accordance with the 
results of the earlier runs. To compare our results directly with those of 
Bassi et a/l., the 1.02-1.85 Bev. and >1.85 Bev. energy groups were added to 
give a > 1.02 Bev. group which corresponds exactly to the energy band reported 
by Bassi to give an anomalously large Cerenkov intensity. This comparison is 
shown in the last column of the table. 
TABLE II 
SUMMARY OF NUMERICAL RESULTs (Part IT) 











Peak Results 
Energy Peak pulse of 
Band range Velocity pulse height Bassi 
number (Bev.) range height (normalized) et al.* 
I 0.084-0 . 27 0.83-0 .960c 12.30+0.30 0.769 +0.020 0.88+0.03 
II 0.27-1.02 0.960—-0 . 9956 15.15+0.17 0.947 +0.013 0.97 +0.03 
III 1.02-1.85 0.9956-0 .9985c 15.90+0.17 0.994+0.014 — 
IV >1.85 >0.9985¢ 16.00+0.13 1.000T — 
III+1V >1.02 >0.9956¢ 15.96+0.10 0.998 +0.005 1.26+0.03 











*Results of Bassi et al. (1952) normalized to theoretical curve at energy band II. : 
tEnergy band IV has been taken as the plateau value consistent with the normalization 
made in Part I. 

It will be seen that our observed value of the most probable Cerenkov light 
output remains constant to within the estimated uncertainty of 1% down toa 
meson energy of 1.02 Bev. Below this energy the most probable intensity 
decreases to somewhat lower values than reported by Bassi et al. 


DISCUSSION OF RESULTS 
If equation (1) is normalized to the plateau intensity (Wimax), assuming ” 
to be constant,? the expression for the relative intensity / of Cerenkov radia- 
tion becomes 


at eee (1-1 Brn’) 
9 =e ead can 
( ) I W, WW max (1 1 n°) 


For convenience in comparing the radiation from particles of different 
masses, the quantity 7 may be plotted as a function of (y—1), the kinetic 
energy of the particle in units of its rest mass, where 
(3) y= (1-8). 

The solid curve in Fig. 5 shows such a plot for Plexiglas, whose index of 
refraction is 1.49 in the visible region. 

The normalized peak position for each of the four high energy groups 
studied in the first part of this experiment is shown as a black bar whose 
width is twice the standard error for the group. The cross-hatched bars show 
the normalized peak positions obtained in the second part of the experiment 
with energy groups III and IV added together. The vertical lines in groups 
I and II show the positions at which the bars should correspond to the theo- 
retical curve when one takes into account the finite thickness of the counter 


2It is only necessary for a comparison with experiment that be independent of y over that 
part of the spectrum where the detector has a significant response. It may be shown that 
within the sensitivity region of the photomultiplier the effect of the variation of m is negligible. 
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Fic. 5. Comparison with theory of various experimental observations of relative Cerenkov 
light intensity as a function of meson energy. Note that the results from different experiments 
have been independently normalized. 


and the sea level meson spectrum corrected for the effect of the absorber 
above the telescope. 

The dotted blocks show the results of Bassi ef al. normalized by fitting the 
effective mean value of their middle energy group to the theoretical curve. 
Mean positions for the various groups are taken from Bassi’s paper. These 
authors call their values ‘‘average’’ values, but since they describe no method 
for measuring very large amplitude pulses, it is probable that these ‘‘average”’ 
values are calculated for a limited region about the most probable pulse height 
as was the case with the A values used for comparison purposes in the first 
part of our experiment. Thus the variation of their ‘“‘average’’ values should 
be consistent with the theoretical intensity—velocity relationship and therefore 
directly comparable with our relative data on most probable intensities. 

The various experimental points in the figure are from the work of Winckler 
et al. (1955) and are for z-mesons of specific energies as shown. We have 
renormalized intensities and have made an approximate energy adjustment in 
order to fit these points to the ‘thin counter’ curve of Fig. 5. This procedure is 
not entirely valid since the variation in Cerenkov emission as the particles 
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slow down in passing through the counter becomes very large at the lower 
energies, an effect that is properly taken into account in the calculated curve 
of Winckler’s paper. On the other hand the shapes of Winckler’s curve and of 
the thin counter curve are very similar except near the threshold, and so 
Fig. 5 demonstrates fairly well the agreement between his experimental points 
and theory. 

From Fig. 5 it will be seen that we do not observe any anomalously large 
Cerenkov light output from high energy u-mesons corresponding to that found 
by Bassi et a/. The apparent constancy to within 1% of the peak positions 
above 1.5 Bev. (bands V, VI, VII, VIII) which has already been noted is in 
agreement with the expected increase of 0.3% in the same energy range. 
The least squares fit of a linear relationship between J and 1/6? to these four 
experimental points gives an increase of (+0.5+1.1)%, which indicates the 
accuracy to which the theoretical curve is verified. 

On the other hand, our low energy groups give slightly less light than 
predicted, group I being 7% low and group II about 23% low. Though this 
may, in fact, indicate a true departure from the theoretical curve in the energy 
region where Winckler ef a/. find some hint of a discrepancy, we believe that 
this apparent disagreement is due to an optical effect in our particular 
Cerenkov counter which is discussed below. 

Cerenkov radiation is emitted at an angle to the particle direction of 
cos~'!(1/8n). For high energies 6 approaches unity and the maximum angle 
becomes cos~!(1/n). For Plexiglas this maximum angle is 47.8°, which is about 
6° greater than the critical angle (sin-'(1/m) = 42.1°) for a Plexiglas—air inter- 
face. Thus in our experiment all the Cerenkov light from the selected high 
energy u-mesons (which enter the counter within 6° of the vertical) is totally 
reflected at the bottom surface of the Plexiglas, and continues outwards to 
the edges of the counter by a series of total internal reflections at the top and 
bottom surfaces. However, as the meson energy decreases, the angle of the 
Cerenkov cone decreases so that, for the oblique mesons first and ultimately 
for all mesons, more and more of the light strikes the Plexiglas—air boundary 
at less than the critical angle and depends on the diffuse reflectivity of the 
MgO to reflect it out to the photomultipliers. This transition region commences 
at the upper end of energy band II and continues throughout groups II and 
I until at the lower end of group I all the initial reflection is diffuse. 

Inaccurate knowledge of the parameters involved prevents an exact calcula- 
tion of the effect of this transition on the observed signal. However, since 
there will be a few per cent loss at each diffuse reflection as compared with 
practically no loss for total internal reflection, and since, on the average, the 
light will be reflected several times before reaching the photomultipliers, it 
would appear that this effect could easily be large enough to account for the 
difference between our observed results and the theoretical curve. The observa- 
tion that in our experiment the fractional discrepancy increases with decreasing 
energy, while in Bassi’s experiment—-where no such transition effect should 
occur—this discrepancy is not found gives qualitative support to this hypo- 


thesis. 
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CONCLUSION 


On the basis of our observations we conclude that there is no high energy 
anomaly in the intensity of Cerenkov light produced by y-mesons passing 
through Plexiglas up to energies somewhat above 4.2 Bev. (v = 0.99970c). 
An apparently low output in the 100 to 1000 Mev. region is attributed to 
optical effects in our counter rather than to a disagreement between theory 
and experiment. 

We believe that our results taken together with those of Winckler eft a/. 
(1955) establish the validity of the theoretical Cerenkov intensity relation to 
within a few per cent from threshold up to a velocity 6 = 0.99970, well into 
the saturation region. 
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HIGH RESOLUTION RAMAN SPECTROSCOPY OF GASES 


VIII. ROTATIONAL SPECTRA OF ACETYLENE, DIACETYLENE, 
DIACETYLENE-d,, AND DIMETHYLACETYLENE! 


J. H. CALLomMon? AnpD B. P. STOICHEFF 


ABSTRACT 


The pure rotational Raman spectra of acetylene (C:H2), diacetylene (C,H2 and 
C,D2), and dimethylacetylene (C,sHe) have been photographed with a 21-ft. 
grating spectrograph. The rotational constants were found to be Bo(C:H2) = 
1.1769 cm™, Bo(CyH2) = 0.14689 cm=, Bo(CyD2) = 0.12767 cm™, and Bo(C4He) 
= 0.1122 cm™. An analysis shows that 7o(C—H) in diacetylene is shorter than 
in acetylene and cyanoacetylene; in dimethylacetylene at least one of the carbon- 
carbon bonds is longer than in the closely related molecule methylacetylene. 


A. INTRODUCTION 


The present paper is part of a program of bond length determinations in 
which the dependence of various bonds on their environment is studied 
(Herzberg and Stoicheff 1955), especially of bonds in non-polar molecules. 
A detailed discussion of the results obtained from the rotational Raman 
spectra of diacetylene (CyH2 and C,yD2) and dimethylacetylene (C4H¢) is 
given here together with a brief report on the rotational Raman spectrum of 
acetylene (C2H»). 

There have been two determinations of the moment of inertia of diacety- 
lene: one from bands in the photographic infrared (Jones 1952) and another 
from the near infrared spectrum (Craine and Thompson 1953), both in 
good agreement. From this value and assumed values of the C—H and C=C 
bond lengths, it was shown that the central C—C bond was unusually short. 
However, in order to establish the validity of the assumed bond lengths, 
the moments of inertia of two additional isotopic molecules must be deter- 
mined. Since a sample of diacetylene-d; was available, the pure rotational 
Raman spectrum was photographed. At the same time, photographs of the 
pure rotational Raman spectrum of the ordinary diacetylene were also 
obtained. 

Although the infrared spectrum of dimethylacetylene has been studied 
several times, the rotational structure of a parallel band has not yet been 
resolved. In the present investigation photographs of the pure rotational 
Raman spectrum were obtained showing a closely-spaced rotational structure 
which is just resolved. From an analysis of this spectrum, a value for the 
large moment of inertia was obtained. This result is discussed in terms of 
probable dimensions of the dimethylacetylene molecule. 

B. EXPERIMENTAL PROCEDURE 


‘Pure’ grade acetylene was obtained from a cylinder supplied by the 
Matheson Company, and sublimed in vacuo before using. The diacetylene 
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was prepared from 1,4-dichloro-2-butyne, and the diacetylene-d2 was pre- 
pared by exchange from heavy water and diacetylene (Callomon 1956). The 
sample of dimethylacetylene was supplied by the Farchan’s Research Labora- 
tories. 

The apparatus used to photograph the spectra was that described in I 
of this series (Stoicheff 1954) with some modifications. Four lamps of in- 
creased length replaced the previous two, resulting in an intensity gain of 
approximately two. A two-lens condensing system (Callomon 1956) was used 
to match the apertures and aperture stops of the 21-ft. grating with the mirror 
system in the Raman tube. 

The spectra of acetylene and dimethylacetylene were photographed at a 
pressure of 30 cm. Hg in exposure times of 10 hours. The gas pressures of 
diacetylene and diacetylene-dz were 15 cm. Hg and the exposure times were 
6 hours. 

Several more intense plates of the diacetylene spectra were obtained in 
30 minutes with the well-known method of placing a cylindrical lens in front 
of the photographic plate of the 21-ft. grating spectrograph (Humphreys 
1903; Bowen 1938). The plano-convex cylindrical lens used in this work is 
made of fused quartz and has a focal length of about 20 cm. Its length is 
15 cm., its height 3 cm., and its maximum thickness about 2 mm. The lens 
is mounted with independent translation and rotation relative to three 
mutually perpendicular axes, the most critical adjustment being the align- 
ment of the cylinder axis at right angles to the rulings of the grating. 

With this lens, the final image height at the photographic plate is reduced 
from 3 cm. to about 2 mm., resulting in an intensity gain of almost 15. The 
length of the lens covers a wavelength region of about 200 A in the second 
order of the 21-ft. grating, and in this region a resolving power of approximately 
100,000 has been obtained with the lens, which is adequate for the present 
investigations. As a check, and in order to detect any errors which may have 
been introduced by the use of the cylindrical lens, the spectra of a number 
of substances, including CyHz and CyDs, were photographed with and without 
the cylindrical lens. These gave identical results. 

C. ROTATIONAL ANALYSIS 

Photographs of the pure rotational Raman spectra of acetylene, the two 
diacetylenes, and dimethylacetylene are shown in Fig. 1. 

(a) Acetylene 


The mean values of the Raman displacements |Av| of the rotational lines 
obtained from the Stokes and anti-Stokes branches are listed in Table I. For 
a linear molecule, the displacements are given by 

|Av| = (4Byo—6Dy)(J+3/2) —8Do(J+3/2)?. 


Values of the constants were obtained by plotting the values of |Av|/(J+3/2) 
against (J+3/2)* in the usual way. The graph is shown in Fig. 2. The inter- 
cept is 4B, (neglecting 6Do) and the slope 8Dp, giving 

By = 1.1769+0.0003, Do = 1.940.4X 10—-* cm“!. 
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Fic. 1. Pure rotational Raman spectra (Stokes branches) of: (a) acetylene, CoHe; (6) 
diacetylene, C,H»; (c) diacetylene-d2, C,D2; (d) dimethylacetylene, CsHg (this photograph 
was made by moving the photographic paper parallel to the spectral lines during enlargement). 
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TABLE I 
OBSERVED WAVE NUMBER SHIFTS OF THE ROTATIONAL RAMAN LINES OF ACETYLENE (C2H2) 














I [Av |, cm7 | St |Av|, cm! 
0 —_-— 13 68 . 206 
] 11.84; 14 72.923 
2 16.44; 15 77.605 
3 21.16, 16 82.32, 
4 25.91; 17 87.01, 
5 30. 58s 18 91.74, 
6 35.305 19 96 . 38; 
7 40.00, 
8 44.73; 20 101.065 
9 49.405 21 105.74; 
22 eee 
10 54.13; 2 115.15; 
11 58.81; 24 —_—_—— 
12 63.53; 25 124.465 
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Fic. 2. |Av!/(J-+3/2) graphs for the ground state of acetylene, diacetylene, diacetylene- 
dz, and dimethylacetylene. 


This value of Bo is in good agreement with values reported from the infrared 
spectrum of C,H», but favors slightly the value 1.1769 cm ~'! obtained by 
Herzberg and Spinks (1934) and Bell and Nielsen (1950) rather than the 
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lower values of about 1.1766 cm~! obtained by Wiggins, Shearer, Shull, and 
Rank (1954), Allen, Tidwell, and Plyler (1956), and Christensen, Eaton, 
Green, and Thompson (1956). For comparison, a value of Do was obtained 
from the simple expression Dy = 4B*/w? (w = 1974 cm~') and also from a 
valence force field calculation. These yielded 1.67 X10-* cm! and 1.60X10-° 
cm! respectively, which are slightly lower than the observed value. 
(b) Diacetylene and Diacetylene-dz 

The displacements of the rotational Raman lines on plates taken with the 
cylindrical lens are listed in Table I] and the plots of Av|/(J+3/2) against 
(J+3/2)? are shown in Fig. 2. The constants obtained were: 

for CyHo: By = 0.14689+0.00004, Dy = (3.2+0.6) X 10-5 cm—', 

for CyDe: By = 0.12767 +0.00005, Do = (2.30.8) K10-° cm. 


TABLE II 


OBSERVED WAVE NUMBER SHIFTS OF THE ROTATIONAL RAMAN LINES OF DIACETYLENE (C,H2) 
AND DIACETYLENE-d2 (C4De) 








Av|, em™ |Av|,cm=! 2 Av|,cm™ 
J J 
CaHe CaDe CaHe CsD. CsHe CsDe 
_ — — —_ \} — a _ — = = ee _ — 
7 4.953* — j} 31 19 .09s* 16.587 | 55 33.153 28.80;* 
8 5.551* 4.875* 32 19. 689* 17.108 \} 56 33.741 29.302 
9 6.166 5.38:* | 33 20.253 17.616 i| 57 34.315 29.829 
10 6.749 5. 882* | = 20.849 4 ae | > a oo 
| 3é 21.425 18.624 | 9 35.462 30.821 
11 7.348 6.385 | 36 22.012 19.117* | 60 36.076 31.350 
12 7 .973* 6.906 || 37 22.595 19.640 | 
13 8.521 7.415 \| 38 23.197 20.130 i} 61 36.648 31.869 
14 9 .076* 7.928 | 39 23.778 20.633 | | 62 37 .27:* 32.364 
15 9 .670* 8.437 | 40 24.365 21.154 1| 63 37. 813* 32.881 
16 10. 286 8.960* || | 64 38 .425* 33.375 
17 10.873 9.445* | 41 21.666 1] 65 39.023 33. 88s 
18 11.464 9 .98.* | 42 22.173 66 39. 582 34.387 
19 12.053 10.456 43 22.687 | 67 40.147 34.893 
20 12.634 10.969 44 i i} 68 om 35.397 
45 23.722 i] 69 41.33; 35.931 
21 13.215 11.470 46 : 24.225 | 70 ——— 36.418 
22 13.795 11.993 1 28.475* 24.727 | 
23 14.391 12.493 48 29 .063 25.243 1] 71 42.503 36.926 
24 15.000 13.006 || 49 29.639 25.756 | 72 ee 37.419 
25 15.566 13.52s* 50 30.225 26.254 | 73 43.676 37 .94s* 
26 16.173 14.034 | 74 38. 460* 
27 16.749 14.537 51 30.813 26.76; || 75 38.979 
28 17.336 15.048 52 31.408 27.271 1} 76 39.481 
29 17.919 15.581 53 31.975 27.772 77 39.982 
30 18. 489* 16.08: 54 32.575 28.297 | 78 40.480 


*Lines blended with grating ghosts. 


The standard deviation of observed line positions (of lines not blended with 
grating ghosts) from the values calculated with the above constants was 
0.02 cm~'. The uncertainties quoted for By and Do refer to random errors 
only. As a check, the spectra of both molecules photographed without the 
cylindrical lens were measured and evaluated. These spectra were somewhat 
weaker and less extensive than those photographed with the cylindrical 
lens. A graphical analysis yielded the following values of the rotational 
constants: 

for CsHe: Bo = 0.14692+0.0001, Do = (4.041) X10-° cm™', 

for CysD2: Bo = 0.12752+0.0001, Do = (3.841) X10-§ cm™', 


confirming the values of By obtained above, but with slightly larger values of 
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Dy. As already mentioned, the spectra photographed with the cylindrical 
lens were more extensive and, for this reason, are considered to give somewhat 
better values of the constants. 

The value of By for CsHe determined here is slightly larger than the values 
obtained by Jones (Bo = 0.14641+0.00015 cm ~') and by Craine and 
Thompson (Bo = 0.1463s+0.001 cm~'). The difference is most likely due 
to the effects of superimposed ‘‘hot”’ bands in the earlier investigations which 
did not permit the evaluation of the centrifugal distortion constant Dp. 

The observed values of Dy were compared with the values obtained from 
the simple formula 4B*/w?, and also from a simple valence force field calcula- 
tion: 


C,D2 


CyHe 


Do (observed) 3.21078 cm“ 2.3X10-§ cm“! 
Do (calculated) = 4B3/w? 1.66 (wo = 874 cm“) 1.13 (a = 858 cm“) 
Using a valence force field 1.64 1.14 





The observed values are approximately twice the calculated ones. Such a 
large discrepancy has also been noticed in the Do value of cyanogen (Mller 
and Stoicheff 1954), where the ratio of observed to calculated values was 
1.7. The cause is not yet understood but a systematic experimental error can- 
not be ruled out altogether. 

In diacetylene, as in cyanogen, there is a low-lying vibrational level (v9) 
at about 220 cm™', and at ordinary temperatures, about 40% of the molecules 
are thermally excited to this level. The B value for this level (B’) is probably 
slightly different from By and the centrifugal distortion constant is most 
likely nearly equal to Do. The rotational Raman spectra of molecules in the 
different vibrational levels are therefore superimposed with a slight relative 
displacement (but in this case not resolved) according to the value of ay. 
An approximate value of ay was calculated according to the method of Shaffer 
and Nielsen (1941) giving —0.00010 cm~—'. The true value of Bo should there- 
fore not differ through this cause from the value derived here by more than 
—0.00004 cm~!. Also, the values of By for CysHz and C,D2 would be affected 
almost equally, so that their difference, which is important in bond length 
determinations, is not affected. 


(c) Dimethylacetylene 

The observed spectrum consists of a series of finely spaced lines which 
are just barely resolved, extending to about 30 cm~' on either side of the 
exciting line. Unfortunately, this region is also overlapped by many grating 
ghosts. Nevertheless, it was possible to measure approximately 80 Stokes and 
anti-Stokes lines, and their wave number shifts are listed in Table III. These 
rotational lines arise from transitions with AJ = 2 and hence are part of the 
S branches. Since the dimethylacetylene molecule is a symmetric top, the 
displacements from the exciting line are given by the equation (in simplified 
form) 


|Av] = 4By(J+3/2) -8D,(J+3/2)3. 
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TABLE III 
OBSERVED WAVE NUMBER SHIFTS OF THE ans RAMAN LINES OF DIMETHYLACETYLENE 
(CyHe) 
J |Av|, cm™ J |Av|, cm7! 
21 10.09 41 acest 
22 10.53* 42 19..53* 
23 10.98 43 19.95* 
24 11.41 44 20.43 
25 11.88 45 20. 86* 
26 12.32* 46 21.30* 
27 12.79 47 21.75 
28 13.25 48 22.20 
29 13.78" 49 22.65 
30 14.14 | 50 23.10 
31 14.57 dl 23.55 
32 15.03 52 23.99 
33 15.48 oa 24.45* 
34 15.93 54 24.88 
35 16.38* 55 25.35 
36 16.83 56 25.82* 
ae 17.27 | 57 26.24 
38 42 58 26.68 
39 18.20* | 59 27.12 
40 18.68* 60 27.58 





*Lines blended with grating ghosts. 


A plot of |Av|/(7+3/2) against (J+3/2)?, shown in Fig. 2, gave the rotational 
constants: 
Bo = 0.11229+0.0002 cm—! and Dy, = 1.381078 cm—!. 
D. BOND LENGTHS 

In using the moments of inertia of isotopic molecules to determine bond 
lengths, the assumption is made that the bond lengths are the same in the 
various isotopic species of a molecule. As is well known, this assumption 
is strictly correct only for the equilibrium bond lengths, 7,, but is often 
applied to effective bond lengths, 7. The latter approximation breaks 
down because of zero-point vibration effects and, as might be expected, these 
are largest for bonds involving hydrogen atoms. Thus, in several investigations 
in which deuterated isotopes were used to determine values of ro(C—H), 
it was found that consistent results could only be obtained by assuming 
differences between r(C—-H) and »(C—D). Differences of 0.002 A were 
found in HCN (Nethercot, Klein, and Townes 1952) and of 0.005 to 0.009 A 
in methyl groups of several molecules (Miller, Aamodt, Dousmanis, Townes, 
and Kraitchman 1952; Costain 1955; Heath, Thomas, Sherrard, and Sheridan 
1955). 

Since pure rotational Raman spectra yield only values of By and not B, 
(hence 7» and not r,), the limitations of the assumption that 7) values are 
equal in isotopic molecules must be kept in mind. 

(a) Acetylene 

The effective internuclear distances ro(C==C) and ro(C—H) for the acety- 

lene molecule were calculated from existing values of the rotational constants, 
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namely Bo(C2He) = 1.1766 and 1.1769 em~!, Bo(C2HD) = 0.9911 and 0.9914 
cm~! (Overend and Thompson 1955; Herzberg, Patat, and Spinks 1934), 
and By(C2De) = 0.8479 cem~! (Saksena 1952; Overend and Thompson 1956). 
With the assumption that ~(C—H) = ™(C—D), the average values were 
found to be ro(C=C) = 1.208; A and ro(C—H) = 1.057. A. When these are 
compared with the most recent values for the equilibrium internuclear dis- 
tances 7,(C==C) = 1.204, A and r.(C—H) = 1.058, A (Christensen, Eaton, 
Green, and Thompson 1956), it is found that 7,(C—H) is slightly larger than 
ro(C—H). This is an indication that the above r»>(C—H) cannot be a true mean 
internuclear distance in the ground vibrational level and hence that ro(C—H) 
does not equal 7)(C—D). An analysis of the existing data on the three isotopic 
molecules shows that 7o(C-—-H) —r(C—D) probably has a value in the range 
0.0005 to 0.0035 A. If the mean value ro(C—H)—7(C—D) = 0.002 A is 
assumed, then ro(C=C) = 1.206; A and »(C—H) = 1.061 A. 


(b) Diacetylene 


The present rotational constants of the two diacetylene molecules give 
the following values for the moments of inertia: : 


Io(C4H2) (190.54+0.06) XK 107-*° g. cm.?, 
Io(C4D2) 219.40+0.09) X 10-* g. cm.? 


These two values together with several assumed values of ro(C=C) ranging 
from 1.203 to 1.207 A were used to calculate 7(C—C) and 7(C—H) and the 
results are shown in the graph of Fig. 3. The solid lines of Fig. 3 show the 


f)(C-H) 
oO 


1.060 
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Fic. 3. Graphical determinations of ro(C—H) and ro(C—C) in the diacetylene molecule 
assuming various values of ro(C==C). The solid lines give the relationships of these parameters 
in diacetylene and diacetylene-d2. assuming ro(C==C) = 1.205 A. The broken line is the locus 
of the points of intersection of corresponding solid lines with other assumed values of ro(C==C). 


relationship between these parameters for each molecule, assuming the value 
ro(C==C) = 1.205 A, the length of this bond in the closely related molecule 
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cyanoacetylene.’ The intersection of the two curves corresponds to the 
values 


ro(C—H) = ro(C—D) = 1.046+0.009 A, 
r(C—C) = 1.376+0.002 A, 
ro(C=C) = 1.205 A (assumed). 


The estimated errors in the bond lengths refer to random errors in By. As 
usual, the value of the CH bond length has a rather large uncertainty since 
its value is very sensitive to small changes in the rotational constants. The 
locus of the point of intersection of corresponding solid lines with other assumed 
values of 7)>(C==C) is also shown in Fig. 3, as a broken line. From this line 
it is observed that a change in the assumed value of 7(C==C) changes only 
7o(C—C) and not the CH bond length, which remains at the low value of 
1.046 A. If it is now assumed that (CH) ¥ r(CD) (shown in Fig. 3 as the 
vertical difference between the solid lines) a larger value can be obtained for 
the CH bond length. For example, if (CH) — (CD) = 0.005 A, as in 
cyanoacetylene,* then 7(C—C) = 1.374 A and r(CH) = 1.057 A. Again this 
is a rather short CH bond length. Thus it would seem that the CH bond 
length in diacetylene is shorter than that in acetylene and cyanoacetylene. 

The bond lengths of the diacetylene molecule given above are in general 
agreement with the values found for similar bonds in other molecules (see, 
for example, Herzberg and Stoicheff 1955). In particular, the carbon-carbon 
single bond, which in diacetylene is adjacent to two triple bonds, is con- 
siderably shorter than other carbon-carbon single bonds. For example, with 
two adjacent triple bonds, the C—C bond has a length of 1.380 A, while with 
only one adjacent triple bond the C—C bond has a length of 1.460 A. This 
result was first shown by Pauling, Springall, and Palmer (1939) from electron 
diffraction data and has also been confirmed from the spectra of diacetylene, 
dicyanogen, and cyanoacetylene. From the present analysis, however, it 
appears that in addition the C—H bonds in diacetylene are slightly shorter 
than in other similar molecules. 


(c) Dimethylacetylene 


The structure of the dimethylacetylene molecule is shown in Fig. 4, along 
with the designation of the structural parameters. Dimethylacetylene is a 
symmetric top molecule with the four carbon atoms in a linear chain. The 





3Westenberg and Wilson (1950) measured the microwave spectra of 10 isotopic species of 
cyanoacetylene involving the atoms H, D, C®, C3, N44, and N®, and used the data to evaluate 
four bond-lengths assuming ro(C—H) = ro(C—D). The bond lengths have been briefly 
recalculated, assuming 79(C—H) not equal to ro(C—D) and, as the values are overdetermined 
from the data, they were adjusted to give the minimum standard deviation (0.10 Mc./s.) 
between the observed and calculated values. The values found were: 


Westenberg and Wilson 





ro(C—H) = 1.069+0.006 A a 
ro(C—D) = 1.064 \ 1.057 A 
ro(C=C) = 1.204;+0.001 1.203 
ro(C—C) = 1.378,+0.002 1.382 
ro(C=N) = 1.157y+0.001 1.157 
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Fic. 4. The structure and structural parameters of the dimethylacetylene molecule. 


moment of inertia about any axis perpendicular to this line and passing 
through the center of the molecule is given by 


Ip = 2[mc(ri? + ro”) + 3my(73" + h? /2)]. 


This expression is valid whether free or hindered rotation is assumed for 
dimethylacetylene. Numerically the value of this moment of inertia calculated 


from Bo is 
Tp® = (249.46+0.45) X10-” g.cm.? 


The evaluation of one moment of inertia is, of course, inadequate for the 
determination of the dimensions of dimethylacetylene; however, since the 
structure of dimethylacetylene is closely related to that of methylacetylene, 
it seemed of interest to calculate the moment of inertia 7, using the methyl- 
acetylene parameters. The values 

ro(C—C) = 1.460+0.002 A, r(C=C) = 1.207+0.002 A, 


ro(C—H) = 1.097+0.010 A, and ZHCH = 108°14+1°, 


found by Trambarulo and Gordy (1950) from the microwave spectrum of 
methylacetylene were used in the calculation and yielded a value J,° = 247.60 X 
10-* g. cm.? This value of the moment of inertia of dimethylacetylene is 
significantly smaller than that determined here from the rotational spectrum 
and suggests that at least one of the bond lengths in dimethylacetylene is 
larger than the value found in methylacetylene. If the change is assumed to 
be wholly in the C—C or C=C bonds, the lengths are found to be 1.469 A 
and 1.221 A, respectively, the remaining parameters being those of methyl- 
acetylene. These two values are to be compared with the bond lengths 1.460+ 
0.003 A and 1.207+0.002 A as found for similar bonds in several molecules 
(Herzberg and Stoicheff 1955). Unfortunately the present data on dimethyl- 
acetylene are not sufficient to determine these two bond lengths independently. 

Other related molecules such as dimethyldiacetylene and dimethyltri- 
acetylene have been investigated by electron and X-ray diffraction methods 
(Bastiansen 1956; Jeffrey and Rollet 1952). While the results appear to be in 
general agreement with values of 1.21 A for the C=C bond length and of 1.46A 
for the C—C bond length, the limits of accuracy range from 0.01 to 0.02 A 
and would completely overshadow any small differences in bond lengths such 
as are indicated by the Raman spectrum of dimethylacetylene. 

The present results on acetylene, diacetylene, and dimethylacetylene show 
that fairly accurate values of the rotational constants of molecules having 
even relatively large moments of inertia can be obtained from their pure 
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rotational Raman spectra. There are, however, certain limitations. Firstly, 
it appears that, while the present techniques are adequate for resolving 
rotational lines having a spacing as narrow as 0.4 cm~', the practical limit 
of resolution is about 0.3 cm~'!, because of line broadening by pressure and 
finally because of the width of the exciting line itself (about 0.25 cm-!). A 
spacing of 0.3 cm~! corresponds to a moment of inertia of about 350 10-4° 
g. cm.* for molecules having nuclei with spin # 0 (and twice this value for 
those with nuclei of zero spin). Secondly, the rotational Raman spectra of 
all but asymmetric top molecules give the value of only one moment of 
inertia. Therefore, some assumptions about the molecular structure usually 
have to be made in order to determine the internuclear distances of molecules 
with several parameters. The number of assumptions can, of course, be re- 
duced by investigating the rotational Raman spectra of isotopic molecules, 
but since rather large amounts of a pure sample must be used in this work, 
one is limited to deuterated molecules at present. 
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RHEOLOGY OF THE EARTH: 
THE BASIC PROBLEM OF GEODYNAMICS'! 


A. E. SCHEIDEGGER 


ABSTRACT 


A discussion is given of the data which are available to provide information 
on the rheological behavior of the Earth’s mantle and crust. It is shown that a 
fundamental distinction has to be made between stress durations that are 
“short’’ (up to 4 hours), “‘intermediate’’ (4 hours to 15,000 years), and “‘long”’ 
(longer than 15,000 years). For short durations, the behavior of the Earth is well 
known from the results of seismological investigations. However, these results 
should not be applied directly to the longer time durations. For intermediate 
time durations, the Earth shows an elastic aftereffect. One obtains a value for the 
(Kelvin) viscosity if it is assumed that the damping of the variation of latitude 
is caused by a Kelvin effect. The same (Kelvin) viscosity is arrived at if the 
strain rebound characteristics of earthquake aftershock sequences are also 
interpreted in terms of an elastic aftereffect. For stresses of long duration, some 
information may be obtained from the uplift of Fennoscandia. It appears that 
viscous creep is the prevailing phenomenon. If this information is combined 
with some further considerations, indications are that the Earth behaves like a 
Bingham body. Then, the fact that orogenetic phenomena are concentrated in 
narrow bands may be due to the existence of solutions of the Bingham equation 
exhibiting rheological instability. 


1, INTRODUCTION 


The chief problem of the science of geodynamics is to determine deformations 
within the Earth and upon its surface. On the surface, the present-day defor- 
mations are known and the task is to explain the latter in terms of stresses 
that could be considered as reasonable. In any definite deformation theory, 
it is a matter of mathematical analysis to determine the stresses from the 
boundary conditions. If it were possible to formulate an exact deformation 
theory applying to the Earth as a whole, it would be, in principle, a straight- 
forward (though not necessarily easy) matter to calculate the stresses from 
the (known) strains and then to look for causes of the latter. Unfortunately, 
it is a fact that the deformation theory applying to the Earth is not known, 
thus leaving the matter of finding the causes of its present-day appearance 
wide open to speculation. It is thus evident that the basic problem of geo- 
dynamics is to determine the rheological conditions in the Earth. The explana- 
tion of the present-day physiography would then follow provided our mathe- 
matical resources were equal to the task. 

As stated above, it is quite clear that the principal difficulty is the problem 
of determining the proper rheological condition for the various layers of the 
Earth, including the criteria applicable for possible fracture. The difficulties 
inherent in finding the rheological conditions applicable to the Earth arise 
from two causes. First, the state of the material in all but the uppermost few 
kilometers of the Earth’s crust is not easy to envisage. Pressures and tempera- 
tures are such that it is unlikely that they can be duplicated in the laboratory 


'1Manuscript received November 26, 1956. 

Contribution from the Dominion Observatory, Ottawa, Vol. 3, No. 4. 

Published by permission of the Deputy Minister, Department of Mines and Technical 
Surveys, Ottawa, Ontario. 








383 











384 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1957 


in the near future, leaving only theoretical guessing to determine the behavior 
of the material under consideration. Second, the time elements involved are 
for the greater part such that, even if experiments involving the correct 
temperatures and pressures could be performed, the human life span would be 
millions of times too short to obtain the desired answers. This, again, forces 
one to speculate. 

The time element in the rheological equations is a concept which, in general, 
has not been sufficiently noted. Terms like ‘‘rigid”’ and ‘‘fluid’’ have a meaning 
only if it is stated what the time intervals are over which the stresses in 
question are applied. In order to survey what can be said about the rheology 
of the material within the Earth, it is therefore necessary to make a distinction 
between various typical duration-intervals for the stresses. It is convenient 
to classify stresses according to whether their duration is ‘‘short’’, “‘inter- 
mediate”’, or ‘‘long’’. We shall consider stresses as short if they have a typical 
duration of the order of 3 seconds, as intermediate if their typical duration 
is of the order of 3 years, and as long if their typical duration is of the order of 
100 million years. It will be noted that the ratio of subsequent typical stress 
durations is 3107. As limits for the interval for which stresses may be 
considered as short, intermediate, or long, one may take the geometric mean 
between the “‘typical’’ durations, i.e. about 4 hours and 15,000 years. 

The ratio of subsequent typical durations, being 3X10", is very large. 
It can therefore not be expected that the rheological behavior deduced for one 
typical duration can be extrapolated for another. The various typical durations 
have to be studied separately. 

It now appears that it is possible to make some progress in determining the 
rheology of the Earth. It is true enough that the exact rheological equation 
cannot be given, but there are a number of data known which permit one to 
set a limit to the boundless a priori possibilities. We shall first discuss some of 
the possible classifications of rheological bodies and then proceed to investigate 
the various time ranges. 


2. RHEOLOGY 

Let us discuss first briefly some of the implications of the science of 
“rheology”. 

In order to determine the motion of a continuous material, there are various 
steps that must be observed. In the first place, one must decide upon a de- 
scription of the deformation. Once this has been achieved, one must express 
various physical laws: the condition of continuity, the law of motion, and 
boundary conditions. 

Finally, one needs a connection between the dynamical quantities (e.g. the 
stresses) and the kinematical quantities (i.e. the strains). The latter is called 
“rheological condition”. 

The schemes applicable to the treatment of deformations of continuous 
matter have been discussed, e.g. by the writer (Scheidegger 1956), in general 
terms. It is, however, very difficult to operate with the general equations 
in particular problems and one is therefore forced to introduce simplifications. 
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Originally, one had only the ‘‘classical’’ materials: viscous fluids and ideally 
elastic substances. The mathematical treatment of these substances has been 
carried to a high degree of refinement; convenient monographs on it have 
been given by Lamb (1945) for viscous fluids and by Love (1934) for elastic 
bodies. 

An extension of the classical bodies to approximate some of the more 
commonly found natural bodies is usually done by, first of all, restricting 
oneself to infinitesimal strains. This, of course, is a severe limitation of general- 
ity and, in fact, usually quite an inconsistent procedure as, if any type of true 
rheological effects occur at all, the displacements become large, at least after 
the elapse of sufficient time. Nevertheless, a very approximate idea of the 
physical possibilities may be obtained by assuming the strains as infinitesimal. 
At any rate, it is customary to do this in most engineering literature, whether 
it is justified or not. A convenient summary of this type of treatment has been 
given by Reiner (1949). We shall present here the more important cases as 
well as discuss some special integrals which will be of signifieance later on. 

In terms of infinitesimal strain terminology, the part of the rheological 
equation referring to the deviatoric components can be written for the Hooke 
solid as follows (where 7,, are the stresses and e€, the strains): 


(2.1) Tr = Quer. (@ x k) 


and for the viscous (Newtonian) liquid as follows: 
(2.2) Ta = 2n€ xx (a e k). 


Combination of both gives (after Reiner) 
(2.3) éx = Tie/(2ua) +7 x / (20m) (i # k). 


A similar combination with appropriate constants can be set up for the iso- 
tropic components. The quantities yy, and ny are constants of the body; 
they are often referred to as (Maxwell) rigidity and viscosity. However, it 
should be noted that the names “‘rigidity”’ and ‘‘viscosity’’ are not very good 
ones. The ordinary viscous fluid is obtained if, in (2.3), uy, — © in which case 
nu becomes the fluid viscosity. Similarly, for ny — © one obtains an elastic 
solid with rigidity uy. It is for this reason that the two constants py and nay 
have been called ‘‘rigidity’’ and ‘‘viscosity’’. However, ‘‘Maxwell constants” 
would be a better name. 

The rheological equation (2.3) describes a liquid which shows stress relaxa- 
tion. That the latter occurs can be seen immediately if it is assumed that the 
deformation is kept constant (i.e. € = 0). Then the stress diminishes expo- 
nentially with a time constant 7, 

(2.4) T = M/uM- 


Furthermore, if a constant stress is applied, it is seen that deformation occurs 
at a constant rate. This phenomenon is called creep. 

Stress relaxation phenomena have been studied extensively by Maxwell. 
Hence the name ‘‘Maxwell liquid”’ for the material at present under considera- 


tion. 
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The equations of the Hooke solid and of the viscous liquid in their infinitesi- 
mal-strain form as displayed in (2.1) and (2.2) can also be combined in a 
different manner, viz. 

(2.5) Tix = 2uxéxt2ncéx (1 ~ k), 


where again a similar combination with appropriate constants could be written 
down for the isotropic components (1 = k). Again, the quantities ux and nx 
are constants of the body which, again, are often referred to as “rigidity” 
and ‘‘viscosity’’, respectively. Again, these names are not very good. One notes 
that a ‘‘Kelvin solid’’ becomes an elastic solid if nx = 0 and a viscous fluid if 
Ux = 0. These are the opposite limits to those required to yield the same 
classical bodies in the case of Maxwell liquids. The use of Maxwell liquids 
and Kelvin bodies has given rise to much confusion. Consequently, we should 
prefer to call ux and nx ‘Kelvin constants’’. 

The Kelvin solid is characterized by an elastic aftereffect; if a stress change 
is performed, the body will eventually reach that state which would correspond 
to Hooke’s law, but only exponentially. The time constant is again n/p. 
A loading-unloading diagram has therefore the characteristics shown in 
Fig. 1. 





LOAD REMOVEO 





Fic. 1. Loading—unloading curve of a Kelvin body. 


The elastic aftereffect characteristic of a Kelvin body also causes free 
oscillations to be damped; as this will be of importance later on, let us try to 
calculate the damping constant. The equation of motion of a Kelvin body is 


(2.6) ux+nx = —cx, 

where c is some constant of the system. Trying for a solution 
(2.7) x = exp[(atiw)f], 

one obtains 

(2.8) utnlatiw) = —c(a®—w?+2atw). 
Equating real and imaginary parts yields 

(2.9) ut+na = —c(a*—w’), 

(2.10) wn = —c2aw, 


and hence one obtains 
(2511) a = —n/(2c), 


as = lee _ a” 
(2.12) w /' qe" 
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Finally, eliminating the constant c yields 
: 0, 2H, 
(2.13) a’ + — +o? = 0. 
n 
The last equation can be used to determine w/y if @ and w are measured in a 
given system. 

Since the rheological equations for Maxwell and Kelvin bodies each contain 
two arbitrary constants, it may appear that either can be used to describe 
any particular phenomenon provided these constants are adjusted properly. 
It can be shown, however, that this is not the case as long as the phenomenon 
considered is transient. The character of the transience is completely defined 
by the type of equation that is used. In order to explain an observed motion 
by one or the other of the two equations, one has to change completely the 
physical cause of the motion. 

Finally, it may be remarked that more generalized bodies are obtained if 
one introduces a yield stress vz into Eq. (2.5). Setting « = 0 to avoid further 
complication, one then has 
(2.14) Ti = Viet2neéu, 


which is the constitutive equation of a ‘‘Bingham solid’’. Here, 7 is called the 
‘‘plastic viscosity’. Further complications can be achieved by including in 
the constitutive equations more and more terms. However, the basic short- 
comings of this type of theory outlined in the beginning of this section are 
thereby not overcome. Nevertheless, the various rheological equations discussed 
here, and commonly employed by engineers, give a useful indication of the 
tvpes of inelastic behavior that might be encountered. 


3. STRESSES OF SHORT DURATION 

According to the above, we consider stresses as of short duration if the 
typical times involved are of the order of 3 seconds—or if the upper limit 
of their duration is about 4 hours. As might well be expected, the rheological 
behavior of the various layers of the Earth is best known for this time range. 
It is found that (infinitesimal) elasticity theory gives a good description of 
the phenomena for small displacements and it is therefore natural to classify 
the findings as to whether they apply to the material being stressed below 
or beyond the elastic limit. 


(a) Stresses below the Elastic Limit 

Stresses of short duration (in the above sense) below the elastic limit occur 
in the Earth during the passing of earthquake waves. Seismology, therefore, 
vields important information bearing upon the question under consideration. 
It is found that the upper layers of the Earth, to a depth of 2900 km., behave 
as an elastic solid, below as an elastic fluid. 

Seismology, directly, yields only the variation of seismic velocity with 
depth. In order to determine the elastic parameters therefrom, one has to 
make certain assumptions regarding the density variation within the Earth; 
Bullen (1953) studied various possibilities. Characteristic values are about 
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2X10" c.g.s. for the rigidity, and 5X10" for the modulus of elasticity in the 
mantle of the Earth. 

Bullen’s estimates start at the layer just below the Mohorovici¢c disconti- 
nuity. Above this discontinuity the crust is very heterogeneous. The elastic 
constants of the various types of rocks for short-duration stresses can be 
measured in the laboratory. A typical value is Young’s modulus E = 4X10" 
c.g.s., Poisson’s number o = 0.25. 


(b) Stresses above the Elastic Limit 
(2) Crust 


A significant number of experiments have been made to determine the 
modes of disintegration of various types of rock samples. The first experiments 
along these lines were done early in the century to determine the modes of 
failure of such rocks as can be used as building materials. After Bridgman 
(1949) initiated the development of high pressure techniques, Griggs (1936) 
and others investigated the deformation of many rock samples under extreme 
conditions. All these experiments must be taken as investigations belonging 
in the present category of short-duration stresses. 

From these investigations it appears that most crustal rocks fracture in a 
manner as postulated by Mohr’s theory; i.e. one observes two preferred 
failure planes inclined at an angle of ¢ < 45° toward the direction of maximum 
pressure and parallel to the intermediate principal stress. The value of the 
stress under which the samples break would be expected to be given by 
an envelope in Mohr’s diagram. As this is rather difficult to tabulate, an idea 
of the magnitude of the stresses involved can be obtained by listing the 
maximum stress differences that the material can support in a state of pure 
shear. A typical value is 10'° c.g.s. (Benioff and Gutenberg 1951). 

A different mode of fracture is observed under pure tension. The tensile 
strength is much below the compressive strength of rocks, presumably owing 
to the presence of Griffith cracks. It is unlikely, however, that such a state 
exists in the Earth anywhere but on the very surface. 


(71) Mantle 

Stress changes of short duration above the elastic limit occur in the mantle 
at the focus of an earthquake during the actual shock. The fact that the 
phenomenon of an earthquake is at all possible indicates that the material 
of the mantle has some resemblance to an ‘‘ordinary”’ solid. The earthquakes 
seem to correspond to sudden shear fractures wherein it is very characteristic 
that most of the displacements along the fracture surfaces are essentially 
parallel to the surface of the Earth (cf. Hodgson 1955; Scheidegger 1955). 

There are indications that modes of sudden displacement different from 
shear fracture may occur in the mantle. The seismic waves emanating from the 
focus of an earthquake do not always permit an interpretation in terms of 
an ‘orange being sliced’. Thus, there are indications that, in certain cases, 
the distribution of waves is such that it would correspond to a multipole of 
higher order. The physical significance of this in terms of an envisageable 
fracture pattern is not yet quite clear. 
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Finally, because it has often been thought that the dissipation of energy 
in seismic waves might be an indication of an anelastic behavior of the Earth 
in the short time range, a remark might be added here regarding this possibility. 
The case, in fact, has been satisfactorily dealt with by Jeffreys (1952, p. 107): 
The damping of seismic waves increases with frequency. Hence it must be 
due to scattering. If it were due to anelasticity, it would have to increase 
with period, and this is not the case. 

4. STRESSES OF INTERMEDIATE DURATION 

The next task is to investigate the rheological behavior of the Earth under 
stresses of intermediate duration. As outlined above, we assume intermediate 
duration to mean from 4 hours to 15,000 years, the typical duration being 
about 3 vears. As might be expected, it is much more difficult to obtain an 
idea of the rheological behavior of the Earth in this time range than for 
short stresses. The information which we do have, thus, comes only from 
two sources: (a) from the damping of the variation of latitude, and (0) from 
the release of stresses in earthquake sequences. 


(a) The Damping of the Variation of Latitude 

An analysis of the observations on latitude variation yields the result that, 
after an annual term (probably due to the movement of air masses) has been 
taken out, a period of approximately 420 days exists. On the assumption 
that this motion is due to random impulses of unknown origin, it is possible 
to say that the damping causes the motion to decrease to a value of 1/e of 
the original between 3 (Jeffreys 1952) and 13 (cf. Bondi and Gold 1955) 
periods. It had originally been thought that this damping would be caused 
by motions in the core of the Earth (Jeffreys 1952). However, in a renewed 
discussion of the problem, Bondi and Gold (1955) claim that the motion of 
the core can not account for the observed damping. The damping would 
therefore have to be due to a Kelvin effect in the mantle of the Earth. Under 
this assumption, it is possible to obtain an estimate for the value of n/u 
occurring in the constitutive equation of a Kelvin body. Indeed, assuming 


(4.1) T = 2r/w = 420 days 


and the damping as of the order of 10 periods 
(4.2) 1/a = —4200 days 


vields for the ratio of the Kelvin constants according to (2.13) 


(4.3) nx/bx = 2 days = 1.72X 10° sec. 


rom the above considerations, the two Kelvin constants cannot be obtained 
separately. One may be tempted to use for ux the “rigidity’’ for the mantle 
as deduced from seismology and confirmed by a discussion of the Earth’s 
bodily tides (Jeffreys 1952), viz. 

(4.4) Mx = 2X10" c.g.s., 

which would vield 

(4.5) nx = 3X10" c.g.s., 
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and then to say that the “rigidity’’ and the ‘‘viscosity”’ of the Earth’s mantle 
have the ‘‘values’”’ shown in (4.4) and (4.5). However, this leads to a fallacy. 
The names “‘viscosity’’ and ‘‘rigidity”’ have no properly defined meaning per se, 
unless it is specified to which rheological equation they refer. There might be 
some justification for using the elastic » for the time range presently under 
consideration, particularly as the discussion of the Earth’s tides gave the same 
result, but then the viscosity nx obtained above in (4.4) and (4.5) is a “Kelvin 
constant” and has no direct relation with the creep viscosity occurring in the 
rheological equation of a Maxwell body. 

The above result that the Earth behaves as a Kelvin body with the constants 
(4.4) and (4.5) in the time range under consideration depends, as indicated, 
on the assumption that the motion in the core cannot account for the damping. 
This, however, is by no means established (Gutenberg 1956). Alternatively, 
it has also been tried to explain the Chandler wobble as caused by a Maxwell 
effect, in which case it leads to a higher (Maxwell) viscosity. 


(b) Earthquake Sequences 

After an earthquake occurs in some part of the world, it often happens 
that the first, ‘‘principal’’ shock is followed by a series of aftershocks of 
dimimishing magnitude. Benioff (1951) has studied a series of such aftershock 
sequences regarding the time dependence of their magnitudes. The ‘‘magnitude”’ 
of an earthquake, as defined by Gutenberg and Richter (1942), permits an 
interpretation in terms of ‘‘strain release’. The total energy of an (elastically) 
strained body is proportional to the square of the strain; hence the strain 
release in a sudden shock should be proportional to the square root of the energy 
release J. The energy release J (in ergs) is connected with the earthquake 
magnitude M by the equation 
(4.6) log J = 12.0+1.8 M. 


Benioff (1951) used the “old’’ equation of Gutenberg and Richter as given 
above for his studies; if the ‘‘new’’ (corrected) equation (Gutenberg and Rich- 
ter 1956) is used, the essence of the results is not changed, as the interest 
lies solely in the relaxation times involved. If the magnitude of each shock 
of an aftershock sequence is known, it is therefore possible to plot a quantity 
(4.7) Su= > Ji, 

which is proportional to the cumulative strain release. If this be done, one 
obtains graphs of the type shown in Fig. 2. 

Such graphs suggest the following interpretation. The principal shock 
occurs when the strain accumulation has reached a critical value. This starts 
a break along an existing or new fault. As soon as this has happened, the fault 
will be ‘‘greased”’ so that any strain accumulating can be released, apart 
for “bumps” that cause the motion to be somewhat discontinuous. One can 
thus picture an earthquake aftershock sequence to be the expression of the 
strain release of a body in which, after loading, the stress has been suddenly 


removed. 
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The resemblance of curves of the type shown in Fig. 2 to the strain-release 
curves of a Kelvin body (Fig. 1) is rather striking and would be even more so 
if the aftershock sequence had not been plotted on semilogarithmic paper. 
The maximum S reached in Fig. 2 corresponds to the unstrained position of 
the Kelvin body, the zero value of S to the loaded position. Naturally, the 
curve of Fig. 2 cannot be expected to come to a complete levelling-off as, in 
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Fic. 2. Accumulated strain release diagram of an earthquake aftershock sequence (after 
Benioff 1951). The curve shows the same pattern as the unloading diagram of a Kelvin body. 
In order to make this comparison, note that the state of zero strain corresponds to 10X10" 
ergs}; the fully loaded position corresponds to S = 0 ergs}. 


addition to the strain release, it must be expected that a constant build-up 
of stress occurs at a slow rate which, presumably, never comes to a stop. 
Thus, after a few years, one would expect the same pattern to repeat itself. 

The study of strain release in earthquake aftershock sequences thus yields 
the result that the latter occur in a fashion that can be explained in terms of 
elastic afterworking of a Kelvin body. Moreover, the relaxation time (i.e. 
the time in which the strain drops to 1/e times the initial value) turns out to 
be of the order of 2 days. This, however, is equal to the ratio n/u of the Kelvin 
constants. It thus turns out that the interpretation of earthquake sequences 
in the above fashion leads to the same value for 7/u as that obtained from the 
damping of the variation of latitude. It would be strange indeed if this coinci- 
dence should be accidental. We can thus summarize the findings of this para- 
graph as follows: Under stress of intermediate duration, indications are that 
the mantle behaves like a Kelvin body with the Kelvin constants of the order 
of magnitude indicated above. 

A final remark should be made regarding the possibility of fractures. 
The occurrence of earthquakes makes it evident that the Kelvin-type be- 
havior subsists only up to a certain limit of stress. If this limit is exceeded, a 
fracture-type process (an earthquake) occurs. What the limiting stress is 
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and what the mode of fracture is (i.e. whether it is, for example, governed by 
Mohr'’s criterion), no-one can say. Caution, however, should be used in apply- 
ing to the Earth strength values from short-duration experiments with rocks 
or such like. In other words, we can not make direct use of the results from 
laboratory experiments to draw conclusions regarding the mode of fracture 
to be expected in the Earth under intermediate-duration stresses. One direct 
piece of evidence is that most earthquakes show displacements essentially 
parallel to the surface, as was already mentioned above. However, since the 
stress state producing the earthquakes is not yet known, the field is wide open 
for speculation. Mohr’s criterion is just one of many that could be valid. 

Finally, a remark should be added explaining why we did not take into 
account any evidence resulting from tidal friction, which, at first glance, 
would fall into the time interval at present under consideration. Jeffreys 
(1952, p. 227) has shown that the damping of the Chandler wobble and the 
tidal friction cannot be due to the same cause. If the time constant obtained 
from tidal friction is applied to the Chandler wobble, assuming the two effects 
as being due to the same cause, the wobble should be damped out in a few 
days instead of in some 12 years. Since Jeffreys has also shown that friction 
in shallow seas can very well account for tidal friction, it is not necessary to 
try to resolve the discrepancy in time constants by assuming some rheological 
peculiarities. 

5. STRESSES OF LONG DURATION 

We finally come to the crucial problem of trying to elucidate the rheological 
behavior of the material in the Earth’s crust and mantle under stresses of 
long duration. By long duration we mean, as stated earlier, from 15,000 years 
upwards, with characteristic duration of 100 million years. Needless to say, 
the behavior in such time intervals is only very vaguely known. 

At the lower limit of this time interval is the rise of some ancient shields 
since the melting of the overlying ice caps which were present during the last 
ice age. It is generally assumed that the ice melted approximately 20,000 
years ago and that, for example, Fennoscandia has been (and still is) rising 
since that time owing to the principle of isostasy, buoyancy providing the 
driving force. 

A discussion of the postglacial uplift of Fennoscandia with respect to its 
possible significance regarding the viscosity in the Earth’s mantle and crust 
has been given by Gutenberg (1941) and by Vening Meinesz (1937). The 
results are summarized in an article by Benioff and Gutenberg (1951). 
Accordingly, the uplift of Fennoscandia proceeds at a rate of the order of 
centimeters per century; “‘viscosities’ (assuming Maxwell behavior) calcu- 
lated therefrom are of the order of (Haskell 1935) 

(a1) "a = 107 c.¢.s. 

(Note that the Maxwell viscosity has no relation to the Kelvin viscosity 
mentioned in the last section.) No value can be given for the other Maxwell 
constant (the rigidity), since all the processes that can be investigated are 
steady-state processes; i.e. t = 0 in Eq. (2.3). The term containing it may 
therefore possibly be neglected, i.e. wy is set equal to infinity. 
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It is also possible to interpret the uplift of Fennoscandia in terms of a 
Kelvin effect, but the physical picture is then somewhat different from that 
above: the driving force causing the uplift is the expression of elastic after- 
working (and not of buoyancy) after the ice melted. The relaxation constant 
in a Kelvin body is equal to nx/ux; if this is set equal to, say, 10,000 years, 
one has in ¢.g.s. units: 

nx/ Bk = 3X10" sec. 


No value for the rigidity can be given. The time constant obtained in this 
manner is obviously not the same as the corresponding one obtained in (4.3) 
for the intermediate time range. One is therefore faced with having to accept 
that, in the two time ranges, either the Earth has the same type of rheological 
behavior (Kelvin type), but with different constants, or else different types of 
behavior (Kelvin and Maxwell) with correspondingly unrelated constants. 
The latter possibility would appear as the more reasonable, especially since it is 
known that many bodies show increasing tendency towards Maxwell behavior 
in long time-intervals. 

The above interpretation of the uplift of Fennoscandia has been criticized 
by Lyustikh (1950), who claims that the figure of viscosity calculated in this 
manner has no meaning. This is because, according to Lyustikh, the uplift 
of Fennoscandia took place also before the last glaciation and therefore cannot 
be caused by isostatic forces. 

However, if we accept that the uplift of Fennoscandia is an expression of 
creep during the past 20,000 years, then it is more likely that there will be 
even more creep in longer time intervals. It is a general experience that creep 
becomes the more pronounced the longer the stresses apply. Thus, we may 
assume that the Earth’s mantle and crust show creep effects if the stress 
relaxation time is long compared with 20,000 years. However, this cannot 
possibly be the whole story. It is well known that mountain ranges can 
persist over several millions of years; in fact the eventual disappearance of 
mountains is not at all due to a creep effect, but to the gradual wearing down 
by the action of wind and water. This implies that the material, to a certain 
depth at least, also must have a non-vanishing yield stress. Introduction of a 
vield stress into the equation of a Maxwell liquid produces that of a Bingham 
body. Noting that, in all processes that are at present accessible to discussion, 
no indication as to the prevailing Maxwell rigidity ,, can be obtained, it 
may thus be permissible to neglect the term containing the latter and there- 
fore to use the Bingham equation as shown in (2.14). 

It remains to attempt to give an estimate for the yield strength vy (7 # k) 
occurring in the Bingham equation. Since the presencc of a yield stress may be 
deduced from the existence of mountains, it appears valid to attempt to extract 
some information regarding its va/ue from the same source. 

The stability of mounds of various materials has been discussed extensively 
by Terzaghi (1943), who postulated the following formula 

v 


(5.2) H = ai N. 


Here, v is the vield stress (to shear), /7 the critical height of the mound which 
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can just support its own weight, g as usual the gravity acceleration, and NV 
a function of the angle 6 of the slope of the mound. This function is shown 
graphically in Fig. 3. 








—> STABILITY FACTOR N 








—— SLOPE ANGLE p 


Fic. 3. Graph of stability factor versus slope angle of a mound (after Terzaghi 1943). 


If we assume, as a crude approximation, that mountains also satisfy Eq. 
(5.2), then we obtain an estimate for the yield stress by putting 7 = 8 km., 
p = 3, g = 980 cm./sec., and N & 6 (corresponding to a slope angle of 45°). 
This yields 


(5.3) vy = 4X10° c.g.s. 


This value is estimated for the surface of the Earth. It is an upper limit if it 
is assumed that mountains would be higher if they could support themselves. 
At a depth, it must be assumed that the yield strength is lower, in conformity 
with the trend of most materials to become weaker if the pressure and tem- 
perature increase. It is interesting to observe that this is of the same order of 
magnitude as the laboratory strength of rocks. 

Thus, combining the phenomenology of the uplift with that of the persist- 
ence of mountains suggests that, under stresses of long duration, a rheological 
equation similar to that of a Bingham body should apply. However, there is a 
further puzzling factor. This is that the orogenetic phenomena occur in bands 
of rather narrow width. In these bands, stresses of intermediate time duration 
are created so that the well-known tectonic effects of folding, faulting, and 
earthquakes occur. The motion of the displacements in these bands seems to 
be in no simple relation to the boundaries of the bands (Hodgson 1955; 
Scheidegger 1955) but generally parallel to the surface of the Earth. For the 
explanation of this, one can think of the following possibilities: 


(i) The bands occur because the crust of the Earth is built up of “stable 
masses”’, viz. the continents and ocean blocks, which are in ‘‘collision’’ owing 
to contraction of the Earth’s interior, continental drift, or some other cause. 
It would naturally be the boundaries of the blocks that would be subject to 
the greatest deformations. A difficulty with this explanation is that the oro- 
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genetic bands are known to extend to great depths as evidenced by deep 
focus earthquakes, whereas it is generally assumed that the “‘stable blocks”’ 
do not extend much below the Mohorovici¢ discontinuity. It thus appears 
that the ‘‘collisions’’ could have no more than a triggering effect for the stress 
concentration, but could not be its actual cause. 

(ii) The bands occur because the material is thixotropic: At the edges of 
the stable blocks a stress concentration is created in the upper regions, which 
causes motions in the intermediate time intervals. These, in turn, change the 
rheological equation so as to make the material more vielding. However, this 
“explanation” is actually nothing but an ad hoc hypothesis and therefore 
just a restatement of the observations. It explains nothing and is merely a 
tautology. 

(iii) The bands occur because the rheological equation for long-term stresses 
allows this type of instability to occur. It must be expected that an equation 
of the Bingham type has integrals in which the deformation becomes unstable 
in narrow bands. The actual position of the bands could be triggered by the 
‘stable masses’”’ as indicated under (i). This is suggested by an analogy with 
the bands of plastic instability observed by Bijlaard (1936). However, not 
sufficient information is available to make this conclusion more definite. 


6. CONCLUSION 


In conclusion, we may emphasize once more that the behavior of the 
material in the Earth’s crust and mantle is substantially different in the three 
time ranges which we have been considering. 

(a) In the short range (up to 4 hours) the material behaves as an elastic solid 
with a rigidity of about 210" c.g.s. and a modulus of elasticity E of about 
5X10" c.g.s. If the elastic limit is exceeded, the material undergoes brittle 
fractures. 

(b) In the intermediate time range (4 hours to 15,000 years) the material 
‘still behaves as an elastic solid with similar elastic constants, but with a 
significant aftereffect. The relaxation time is about 2 days, the ‘‘viscosity”’ 
(Kelvin constant) nx about 3X10"? c.g.s. It appears that there is a critical 
stress which, if exceeded, causes a fracture-type process. The limits of strength 
referring to this fracture-type process are not known. The process itself is 
probably very similar to that discussed under (a) as, once a fracture starts 
to move, this occurs in a time interval characteristic of (a). 


(c) In long-duration stresses, creep is the significant phenomenon, i.e. the 
material yields indefinitely at a constant rate if the stresses stay constant. 
The relaxation time appears to be of the order of 20,000 years upward. The 
‘“‘viscosity’’ (Maxwell constant) ny is probably about 10” to 10% c.g.s. From 
the relative persistence of mountain chains it may also be assumed that a 
low yield stress is present so that the rheological equation is essentially that 
of a Bingham body. The concentration of orogenetic phenomena in narrow 
bands may be due to rheological instability. For easy reference, the results 
arrived at in this paper are put together in Table I. 








396 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1957 


TABLE I 
Time interval Applicable to: Behavior Constants 
Up to 4 hours Seismic wave Elastic solid uw (rigidity) 2X10" ¢.g.s. 
(“‘short’’) transmission with strength E (Young’s modulus) 
limit 5X 10!? c.g.s. 
Strength limit 10! c.g.s. 
4 hours to 15,000 vears Chandler wobble Kelvin body ux (‘“rigidity’’) 2K 10" c¢.g.s. 
(“‘intermediate’’) Earthquake with strength nx (‘‘viscosity’’) 
aftershocks limit 3X10" c.g.s. 
Faulting Strength limit: present, 
Earth’s tides value unknown 
Longer than Large scale orogenesis Bingham body uy (‘‘rigidity’’) term = con- 
taining it may possibly be 


15,000 years (“‘long’’) Uplift of 
Fennoscandia neglected 

Polar wandering nu (‘viscosity’) 102? c.g.s. 

v (yield stress) 4X 10° c.g.s. 
at surface 


Any investigation of geodynamic problems should take the above facts 
into consideration. In general, the creation of an orogenetic system will take 
place during millions of years; it is therefore not permissible to apply rigidities, 
‘viscosities’, etc. calculated from short- and intermediate-duration data to 
global considerations concerned with the distribution of orogenetic elements. 
Much unsatisfactory work has had its origin in this fundamental mistake. 

On the other hand, the /oca/l tectonics is governed by the laws of the inter- 
mediate time-range: During the over-all orogenetic process, the local stresses 
may be assumed to reach large values in periods substantially less than 
15,000 years. If a certain critical stress-difference is exceeded, a fracture-like 
process occurs which is felt as an earthquake and can possibly be recognized 
physiographically as a fault. 

The laws of the short time-range seem to have no application to geodynamics. 
They appear to be of importance solely in the propagation of elastic waves 
through the Earth. 
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VARIATIONS IN THE AMOUNTS OF POSITRONIUM FORMED 
IN LIQUIDS AND AMORPHOUS SOLIDS! 


R. E. GREEN? AND R. E. BELL? 


ABSTRACT 


Measurements are reported of the intensity, J2, of the long-lived (72) com- 
ponent of positron annihilations in a number of liquids and amorphous solids. 
Large variations in Jz have been found with different samples of similar sub- 
stances, It is shown, for example, that J, for fused quartz of good quality is 
about 51% of all the positron decays, and that an earlier published value of 
29% probably referred to a 96% silica glass rather than to pure fused quartz. 
The new value, taken with the results of angular correlation measurements, 
gives convincing evidence that the annihilation of triplet positronium in con- 
densed materials proceeds mainly by the ‘‘pickoff’’ mechanism. It is shown that 
the addition of NO;~ and NO2- ions to water in increasing concentration de- 
presses J, from its normal value of (21+5)% to a value at least as low as 2% 
without any marked change in rz. This behavior is explained in terms of the 
formation of the positron compound e*NO;~ in the water. This hypothesis 
explains the shape of the curve of Jz as a function of NO;~ ion concentration 
and also gives a cross section for the capture of positrons by NO;~ in rough 
agreement with expectation. 


INTRODUCTION 


Fhe measurements described here show that the amount of positronium 
formed when positrons are stopped in certain liquids and amorphous solids 
is much less constant than has hitherto been supposed. Previously Bell and 
Graham (1953) had detected a long-lived (72) component in the decay of 
positrons in a number of amorphous substances, and attributed it to the 
formation of triplet positronium; this conclusion is now widely accepted. 
They measured the intensity, J2, of this component with fair accuracy in 
only four substances, fused quartz (Jz = 29%), polystyrene (36%), poly- 
ethylene (29%), and fused borax (32%). The error attached to each of these 
figures was +5%. For a large number of other samples, Bell and Graham 
refrained from quoting definite figures for 72, remarking only that the value 
of I, for these substances seemed ‘‘consistent with the figure of ~30%”’ 
found for the four samples referred to above. They also mentioned that J» 
for teflon was about 33%. On these grounds they assumed that all substances 
showing the rz component would also show an intensity J: of around 30%. 
The impression that this is true has been strengthened by subsequent measure- 
ments by other observers as follows: polystyrene (27+3)% (Ferguson and 
Lewis 1953), benzene (31+43)% (Berko and Zuchelli 1956), teflon (814+5)% 
(Gerholm 1956). It will be noticed, however, that these measurements have 
added only one new substance, benzene, to the list. Graham, Paul, and 
Henshaw (1956) have found a component in liquid helium with a mean life 
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) of 2.7X10-* sec. and an intensity greater than 30%, but their more recent 
work (Graham, unpublished) seems to show that this is not the usual 72 
component. 

This accepted constancy of J. led to apparent difficulties in the inter- 
pretation of the experiments, but the additional measurements to be re- 
ported here show that this constancy was largely due to the chance selection 
of samples. 

In a very recent publication, Landes, Berko, and Zuchelli (1956) have 
shown that J, is quite different in molten naphthalene (29%) from what it is 
in crystalline naphthalene (9%). This work appears to be the first recognition 
of large changes in J». The present results, however, are concerned only with 
amorphous substances, including liquids. Berko and Zuchelli’s value of 
29% for molten naphthalene adds one more substance to the list of amorphous 
samples having J» near 30%. 

EXPERIMENTAL PROCEDURE 

For the measurements, sources of Na®* were dissolved or embedded in the 
sample materials and viewed by two diphenylacetylene 6342. scintillation 
counters. One counter was biased to accept only pulses from the 1.28 Mev. 
gamma rays emitted practically simultaneously with the positrons from the 
source; the other was differentially biased to accept only the upper part of the 
Compton electron pulse distribution from the 511 kev. annihilation quanta. 
The distribution in time of the annihilation quanta relative to the 1.28 Mev. 
gamma rays was measured with a time-to-amplitude converter (Green and 
Bell 1956) and a 28-channel kicksorter (Skarsgard and Bell 1954). The time 
calibration of the apparatus was performed by inserting known delays by 
means of cables, and observing the shifts of the centroid of a prompt coinci- 
dence distribution. The use of this converter, which we shall describe in another 
publication, enabled the whole delay distribution to be measured during a 
single counting interval. With a source of 10 uc. of Na®?, a time spectrum 
containing a total of over 200,000 counts could be recorded in about 2 hours. 
We call such a measurement a long run, and one having about one-tenth as 
many counts a short run. 

The Na” source material was prepared by a 45 Mev. proton bombardment 
of “Superpure”’ aluminum in the McGill synchrocyclotron, with subsequent 
chemical separation of Na®? by ion exchange. The specific activity after 
separation was about 0.1 uc./ugm. About 10 uc. of the dissolved material 
was deposited over an area of 1 cm.” on a mica sheet of thickness 1.8 to 2.0 mg./ 
cm.? and diameter 4 cm. and evaporated to dryness. A similar sheet was 
placed over the source and the two sheets were waxed together at the edges 
with black sealing wax. This source will be referred to as the sandwich source. 
Mica was chosen for the sandwiching material because of its chemical inert- 
ness and because it shows no long-lived component of positron decay. A 
certain fraction of the positrons from the source will annihilate in the mica, 
reducing the measured intensity 7, proportionately. This fraction can be 
estimated from a curve published by Rosenfeld, Swanson, and Warshaw 
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(1956) as something below 10% for the 1 Mev. positron spectrum of C'. 
Owing to the effects of backscattering, it will vary somewhat with the effective 
atomic number of the surrounding sample. Taking account of the lower 
energy (0.57 Mev.) of the Na*® positron spectrum, we correct J» for this 
effect by multiplying by 1.10 for water and organic substances, and by 1.15 
for quartz and similar substances. The effect of these corrections on the con- 
clusions of the experiments is small. 

In one measurement on water, the Na*? source was dissolved directly in 
distilled water contained in a lead cell. The only correction to J. in this case 
is for the fraction of positrons escaping from the water and annihilating in 
the lead walls. The correction factor to Jz was estimated from the dimensions 
of the cell as 1.04. In another case the active water was in a Pyrex cell; since 
Pyrex glass is fairly similar to water in behavior (see below), no correction 
was applied for annihilation in the glass walls. 

One sample of fused quartz, made from ‘‘Vitreosil’’ (stated to be >99.8% 
SiO:), was directly activated by a 70 Mev. proton irradiation in the cyclotron, 
the reaction being 

Si? (p, 4p3n) Na”. 


After 2 weeks’ decay, the sample was found to contain no activity other than 
Na®* in a detectable amount. The Na? in this case was distributed throughout 
the body of the fused quartz sample. The sample was surrounded by further 
fused quartz so that no positron could reach any other substance; this source 
therefore does not require any correction to the measured value of /2. The 
radiation damage to this piece of quartz was undoubtedly appreciable, since 
about 10'* protons per cm.? passed through it at 70 Mev. The discoloration 
of the quartz was not very marked, however, and it is likely that the heating 
effect of the proton beam on the sample caused at least a partial annealing 
of the radiation damage while the bombardment was in progress. This view 
is supported by the fact that the part of the quartz sample nearer to the 
water-cooled target clamp, and just outside the direct proton beam, showed 
more discoloration than the part that was more heavily bombarded and less 
well cooled. We refer to this specimen as the activated quartz sample. 
RESULTS 

Fig. | gives the results for three different high-silica glasses. Curve (a) 
shows the time delay spectrum for the activated sample of fused quartz, 
measured in a long run. The 72 component in this curve is extraordinarily 
prominent, and its decay is observed over a factor of 1000 with no sign of 
departure from the exponential shape. The mean life, 72, is (1.53+0.05) X 
10-° sec. Fig. 1 also shows the method of evaluating the intensity J. of the 
tT. component. The straight line representing the rz component is extrapolated 
back to the point where the center of a prompt resolution curve would lie 
(point P in Fig. 1); for the samples examined in this paper, this point is 
3.5 X 107! sec. to the left of the center of the 7; peak of the quartz curve, since 
the 7; mean life involved is 3.5X10-!'° sec. with an uncertainty of perhaps 
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Fic. 1. Time spectra of the decay of positrons in three different high-silica glasses: (a) 
Vitreosil, >99.8% pure fused quartz; (b) Vycor, 96% fused quartz; (c) Pyrex, ~80% fused 
quartz. Curve (d) shows the shape of the curve observed when the Na”? source is replaced 
by Co®, whose gamma-gamma coincidences are essentially prompt. The curves were taken 
at different times and their centers do not coincide exactly. The point ? of curve (a) is used in 
evaluating J2; see the text. 


+1107!" sec. (Bell and Graham 1953). For a value of r. of 1.5 107° sec., 
this uncertainty will cause an error of only +7% of the value of J». The point 
P having been located, the value of J. is given almost exactly by the area under 
the 7. line to the right of P, divided by the area of the entire curve (i.e. the 
sum of all the counts in the individual channels). A small further correction is 
required because of the finite width of the prompt response curve compared 
with the value of zr». This correction is analogous to the correction for finite 
duration of the counting interval when a rapidly decaying source is being 
counted. It is easily shown that the correction factor is given by (r0/rT2) X 
[sinh (r9/r2)]~', where 7, is the half width at half height of the prompt resolution 
curve. For the times involved here, this factor is usually about 0.95, and 
varies little from sample to sample. The value of 7» can be judged well enough 
from the width of the 7; peak in the curve. The same procedure for evaluating 
ZT. was followed in all the curves. The statistical errors in evaluating J» are 
seldom large, but the systematic errors are appreciable, and we have purposely 
assigned rather generous errors to J». The conclusions of the experiment are 
not affected by the errors assigned. 

The value of J. for the activated quartz sample, evaluated from Fig. 1, 
is (53+8)%. It is noticed at once that both 72 and J, for this sample are 
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quite different from those reported by Bell and Graham for fused quartz 
(v2 = 1.8X10-9 sec., Zp = (2945)%). It becomes clear on enquiry that the 
samples of ‘‘fused quartz’’ that one finds in the Jaboratory are not very well 
defined in composition. Taken with the variability in structure owing to the 
method of preparation of the quartz, it is perhaps not surprising that both 
72 and J» should be found to vary with different samples. 

Fig. 1 also shows the results of short runs made with the sandwich source 
on a pair of “Vycor’ quartz windows (96% SiOz) in curve (6), and on a 
sample of Pyrex glass (~80% SiO») in curve (c). The Vycor has 72= 1.75 X 10~® 
sec., Ig = (82+6)%, and the Pyrex 72 = 1.22X10-° sec., [2 = (8547)%. 
Thus the three samples of Fig. 1, all quite similar substances, show among 
themselves both a variable 72 and a variable J». It will be noticed that the 
Vycor result is close to that given by Bell and Graham for their ‘‘piece of 
fused quartz laboratory ware’’, which in fact is often made of Vycor. Curve 
(d) of Fig. 1 shows the shape of the curve for prompt coincidences from a 
Co® source, using the same biases on the scintillation counters as in the 
positron experiments. 

Samples of “‘Homosil”’ and “‘Ultrasil’’ optical-grade fused quartz, the latter 
stated to be the best available for ultraviolet transmission, showed J» = 
(49+8)%, essentially the same result as the activated quartz sample. We 
conclude that fused quartz of good quality has J; = (5146)%, and that 
even a 4% addition of other oxides to silica in a glass, as in Vycor, seriously 
depresses the value of Jo. 

Fig. 2 shows the results of three short runs for distilled water. Curve (a) 
is for water with Na”? dissolved in it, held in a Pyrex cell, curve (0) is for the 
same kind of source in a lead cell, and curve (c) is for the sandwich source 
immersed in water. The three curves are all very similar, and after correction 
for wall loss for (6), and mica loss for (c), they yield the following values of 
To: 


(a) 22.3%, — (b) 21.9%, () 20.6%. 


The value of J. for water may therefore be given as (21+5)%. Earlier runs 
on the sandwich source in water using the conventional single-channel delayed 
coincidence apparatus with 27) = 2X10~° sec. also gave values of 7, around 
21%, in agreement with the above figures. In all cases the indicated value 
of rz was about 1.80X10-° sec., in reasonable agreement with Bell and 
Graham's value of 1.710~-° sec. Analysis of Bell and Graham's published 
diagram for water also yields J, about 20%. The agreement of J. from curves 
(a) and (6) with that from curve (c) is evidence that the small correction 
factor for mica loss adopted here is sufficiently accurate. 

A measurement with the sandwich source in p-dioxane gave 72 = 2.671079 
sec., J2 = (8647)%, in almost exact agreement with Berko and Zuchelli’s 
results (1956) on benzene. The identity of these results is interesting because 
the p-dioxane molecule has a benzene ring structure. 

A very interesting effect is observed with water solutions of sodium nitrate 
in various concentrations. Fig. 3 shows a set of results for pure water, and 
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POSITRONS IN WATER 


COINCIDENCES 





° 2 4 6 8 10 (2 14 
TIME, 10°? sec. UNITS, ARBITRARY ZERO 


Fic. 2. Time spectra of the decay of positrons in water: (a2) water in a Pyrex cell, (6) water 
in a lead cell, (c) water surrounding the mica sandwich source. The statistical errors are shown 
only for curve (a). Those for curve (c) are similar, and those for curve (b) are somewhat smaller. 
The curves were taken at different times and their centers do not coincide exactly. 


for NaNO; concentrations of 0.25 M, 0.5 M, 2.5 M, and 5 M.* It will be 
seen that the effect of adding NaNO; to water is to depress the value of J2 
without substantially changing 72, At the highest concentrations, 2.5 M 
and 5 M, the results become uncertain statistically, and the values for J 
are only rough estimates. Any defects in the experimental techniques at 
these low levels will give values of J, that are too high. Fig. 4 summarizes 
the results of Fig. 3 in a curve of the experimental value of J. as a function 
of concentration of NaNO;. Two points in Fig. 4 have been taken from runs 
that were omitted from Fig. 3 in the interests of clarity. Fig. 4 shows that 
I, falls very rapidly as NaNO; is first added to the water, but seems to de- 
crease more slowly as the concentration is made large. The curve shown in 
Fig. 4 is drawn according to equation (5) of the discussion section. 

The behavior of NaNO; is presumably to be attributed to the NO;~ ion, 
because sodium hydroxide, for example, shows no effect on J, and only a 


*M = moles per liter of solution. 
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E1G. 3. Time spectra of the decay of positrons in water solutions of sodium nitrate of con- 
centration 0 MW, 0.25 M, 0.5 M, 2.5 M, and 5.0 M, showing the decrease of J» as the concen- 
tration is increased. A change of range of the time-to-amplitude converter occurs between the 
0.5 MV curve and 2.5 M curve. 
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Fic. 4. Points: experimental values of 2, observed from runs like those of Fig. 3, as a function 
of sodium nitrate concentration. Curve: theoretical shape computed from equation (5), fitted 
to I, = 21% for pure water (arrow) and to J, = 10.56% for 0.25 M NaNO. 


small effect on 7. at concentrations as high as 10 M. On the other hand, 
nitric acid (HNOs) shows the same effect as an equal concentration of sodium 
nitrate. It has also been observed that sodium nitrite (NaNO,) shows almost 
the same effect as sodium nitrate, differing only in requiring a slightly higher 
concentration for the same degree of suppression of /2. 
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Several other salts were found to depress the value of J. slightly in con- 
centrated water solution, but in every case the concentration was high enough 
(2-10 M) to affect the bulk properties of the water (mass density, electron 
density, etc.). Some salts leave J: relatively unaffected and modify 72 to a 
small extent. Table | is a list of the salts tried, their concentrations, and 
their measured values of rz and Jz. None of them possesses the apparently 
highly specific effect of suppressing positronium formation in dilute solution 
shown by NaNO;:, HNQO;, and NaNOs, although a wider search would 
presumably reveal other substances with a similar behavior. 


TABLE I 


EXPERIMENTAL VALUES OF 72 AND J. FOR VARIOUS SALTS DISSOLVED IN WATER. ALL MEASURE- 
MENTS WERE MADE WITH THE MICA SANDWICH SOURCE. THE ASSIGNED STANDARD ERROR IN 72 
Is +0.10, AND IN J, IS +5 IN ABSOLUTE VALUE, OR +3 FOR COMPARISON OF DIFFERENT 

ENTRIES IN THIS TABLE AMONG THEMSELVES 























Entry No. Sample Concentration t2 (107° sec.) Ie (%) 
1 H,0 Pure 1.80 21 
2 HCl 5M 1.83 16 
3 NaCl 5.2 M 1.76 15 
4 KCl 3.9 M 1.68 17 
5 NH,Cl 5.2 M 1.72 14 
6 ZnCle 6.3 M 1.82 16 
7 BaCle 1.7 M 1.72 17 
8 AICI; 3.9 M” 1.71 20 
9 LaCl; Sat. 1.80 18 

10 NaOH 10 M 1.49 21 

11 NaC.H;02 7.6 M 1.91 19 

12 K.CO; 5.8 M 1.57 17 

13 Na2CO; 1.3M Live 20 

14 Na2SQ, 3.3 M 1.75 16 
DISCUSSION 


The preceding results have shown that it is possible to find substances all 
having roughly the same value of 72 but with J: varying from ~2% to 53%. 
In terms of positronium formation by the incident positrons, these figures 
represent only the triplet positronium formation, i.e. three-quarters of the 
total; the total formation therefore ranges from ~3% to 71% of all the in- 
cident positrons. The high value for fused quartz has an immediate application 
in the interpretation of angular correlation experiments on the annihilation 
radiation. Page et a/. (1955, 1956) have found a narrow component in the 
two-photon angular correlation curve for fused quartz; a similar effect has 
been noted by Stewart (1955) for teflon. Page et al. evaluated the fraction 
of all their coincidences appearing in the narrow component, and found 
(17.641.5)%, later given simply as 18%. A narrow component can arise 
only when a positron in a slowly moving positronium atom annihilates with 
the electron to which it is bound. If the 7. component in fused quartz decayed 
mainly by conversion to singlet positronium, as suggested by Bell and Graham 
(1953), then essentially all the positronium formed would decay from the 
singlet state, and the fraction appearing in Page’s narrow component would 
be equal to the total positronium formation. On Bell and Graham’s original 
result for what they called fused quartz, this would be 4(29%) = 39%. If, 
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on the other hand, the rz component represents annihilations of triplet posi- 
tronium mainly by ‘‘pickoff’’ annihilation with electrons bound to neighbor- 
ing atoms or molecules (Garwin 1953), only the singlet positronium would 
yield a narrow component, and again on Bell and Graham’s original result, 
the fraction in the narrow component would be 4(29%) = 10%. Page et al. 
were puzzled by the failure of their narrow-component fraction, 18%, to 
agree with either value predicted by Bell and Graham’s result. Assuming 
pickoff, the discovery in this work that good fused quartz has J, ~ 51% 
allows us to predict a narrow-component fraction of 4(51%) = 17%, in agree- 
ment with that found by Page et a/. This agreement is strong evidence that 
the pickoff mechanism is in fact the main process determining the lifetime ro. 

The effect of NaNO; in suppressing the formation of positronium is difficult 
to discuss conclusively. A concentration of 0.1 M, corresponding to less than 
two molecules of NaNO; per thousand water molecules, produces a distinct 
reduction in J, (see Fig. 4) without appreciably affecting 72. It must be em- 
phasized that this strong effect is a quite different effect from the ortho-para- 
positronium conversion found by Deutsch (1951) for NO and NOz in gases. 
The present effect is not one of conversion, which would be reflected mainly 
as a variation of r. rather than Js, but one of inhibition of the formation of 
positronium. Dulit, Gittelman, and Deutsch (1956) have observed that 
the amount of positronium formed in argon can be influenced by adding small 
amounts of a polyatomic gas; their results, however, show an increase in 
positronium formation as the concentration of the polyatomic gas is in- 
creased, followed by a decrease to the value characteristic of the polyatomic 
gas itself. They do not give quantitative figures. This behavior can be under- 
stood qualitatively (Deutsch 1955) on Ore’s picture (1949) of the energetics 
of positronium formation. In a similar way, Deutsch and Brown (1952) and 
Deutsch and Dulit (quoted by Deutsch 1955) have shown that the enhance- 
ment of positronium formation in argon caused by an applied electric field 
can be nullified by the presence of a small amount of a polyatomic gas. 

No similar simple explanation of the present results for NO3;- and NO.- 
in water seems to be available. In the first place, it would require great bold- 
ness to apply Ore’s picture to the complicated liquid medium and it could 
only be expected to be a very general guide. In the second place, water itself 
is polyatomic, as well as the additives NO;~ and NO.~, and it is difficult to 
see how the additives could reduce the amount of positronium formed so 
drastically—according to Fig. 4, by a factor of 10 or more. It would also be 
difficult to explain why the other polyatomic ions from the salts appearing 
in Table I do not have a comparable effect; these other ions are the acetate, 
sulphate, carbonate, hydroxyl, and ammonium ions. Of course none of the 
above arguments can rule out the energetics of positronium formation as an 
influencing factor on J2, but they do show that an additional mechanism is 
necessary to explain the results. 

An alternative hypothesis is the following. Since HNO; is a stable com- 
pound, we may suppose that a positron may replace the proton and still 
form a stable compound. This compound would be written e+NO;7-; the paral- 
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lel notation for nitric acid would be H+NO;-. The work of Ore (1948, 1952) 
is not conclusive on the subject of the stability of such positron compounds, 
but he has demonstrated that e+H™ is stable in free space, and Simons (1953) 
has shown the same for e+Cl~-. We therefore assume provisionally that e*>NO;~ 
is likely to be stable, or at least that e+NO;7 is not dissociated by water within 
a few times 107-!" sec. after its formation. Given these conditions, NO; ions 
in water would capture positrons to form e+NO;-. The lifetime of this com- 
pound against annihilation would be of the same order as that of singlet posi- 
tronium, for the positron would be situated in a comparable density of electrons 
of favorable spin orientation. Thus the positrons captured by the NO;7 ions 
would be lost from the 7, component in the decay curves and would appear 
instead with a short lifetime of the order of 7;. We now show that this hypo- 
thesis explains the shape of the curves of Fig. 4. 

The fate of a positron in pure water—whether it will form positronium or 
not—is presumably decided within some average time of decision short com- 
pared with the 7; lifetime of about 3X107-'° sec. We call this short time Tp. 
During this time interval there is repeated capture and loss of electrons by the 
positron, and speaking in terms of averages, there is a probability IT, per 
unit time of its forming positronium and a probability I, of its remaining 
free or annihilating in a free collision. The total probability of forming posi- 
tronium is I’,7’p and the total probability, (T+ 1,)Tp, is unity. The fractional 
formation of positronium is T',7p, and the 7. component represents three- 
quarters of this, so that we have, from the water experiments, 


(1) I, = 3(T,Tp) = 0.21. 
Upon adding NO ;~ ions to the water, we introduce a positron capture process 
with a rate Nod where N is the density of NO;~ ions, ¢ is the mean capture 
cross section, and @ is the mean relative velocity of the positron and the NO;— 
ions. We then have, instead of (1), 


(2) I, = 2 T,Tp/(T,+T s+ Nod)T p. 
With the use of (1) and the fact that (T,+T,)Tp = 1, we get 
(3) T, = 0.21/(1+ Nei T)). 


If the concentration of NO; is varied, we vary only N in equation (3), so 
that we have 


(4) I, = 0.21/(1+RM), 


where / is the molarity of the NaNO; solution, and & is a constant. The 
experiments show that J, is about half normal, i.e. 1+] is about 2, when 
M = 0.25 (Fig. 4), and we get finally 

(5) IT, = 0.21/(1+4M). 


The solid curve in Fig. 4 is a graph of expression (5); the agreement with 
the measured points is good except at the two highest concentrations where, 
as already remarked, the measurements of 72 are unreliable and the experi- 


mental values may be too high. 
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A further confirmation of the NO;~ capture mechanism consists of showing 
that the positron capture cross section ¢ of NO3~ is of a reasonable magnitude, 
viz. of the order of the geometrical cross section of the NO 37 ion. In. order 
to estimate ¢ we must estimate Tp. If positronium formation occurs during 
the time when the positrons are slowing down from, say, 1000 ev. to 10 ev. 
as suggested by the electron capture calculations of Jackson (1955), and if 
the inelastic processes by which this slowing down takes place have cross 
sections of the order of a2, then the time TJ p is of the order of 101% sec. The 
positron velocities over this interval are around 4X10* cm. sec.~! We have 
seen that 1+M is equal to 2, i.e. Nad Tp = 1, when M = 0.25, i.e. when 
N = 1.5X10?9 cm.~* We get from these values of N, 3, and Tp, 

(6) ¢ = [NoT>p}" = 1.7X10-'* cm.? 


This value of ¢ is about equal to 27a,’, in the expected range of magnitude. 
The above estimates are probably rough enough to allow a factor of 10 change 
in either direction, but the order of magnitude would still support the hypo- 
thesis of et NO;- formation. 

It is now necessary to consider the substances in Table I, which do not 
show the positron compound effect. We can divide the entries in Table | 
into. four classes: entries 2 to 9, chlorides; entry 10, hydroxide; entry 11, 
acetate; and entries 12 to 14, salts having divalent negative ions. We dispose 
of the last group at once by noting that it would be necessary to have two 
positrons present simultaneously in order to produce a positron analogue of, 
say, NasSO,y. The remaining three negative ions, Cl-, OH~, and C.H;0.-, 
must now be supposed not to form stable compounds such as etCI~ in the 
presence of water; these compounds may of course be stable in free space, 
as indeed Simons has shown for e+CI~ (1953). The behavior of the NO;~ and 
NO.- ions now seems much less remarkable than it did at first sight, because 
they represent two out of five monovalent negative ions tried in the experi- 
ments. 

We may note in passing that in an experiment like that of Pond (1954), 
where the defect in the two-quantum annihilation rate (i.e. the three-quantum 
rate) was measured in benzene as a function of the concentration of a dis- 
solved free radical, the quantity measured is proportional to the product 
tl. In Pond’s experiment, the three-quantum rate decreased with increasing 
concentration of the free radical because 72 was being reduced by the action 
of the free radical in causing ortho—para—positronium conversion (Berko and 
Zuchelli 1956). An experiment like Pond’s performed on water solutions of 
NaNO; should show a similar effect because of the decrease in J2, and the 
shape of the curve of three-quantum rate as a function of concentration 
should not be very different. The experiment could also be done by direct 
detection of the three-quantum annihilations in the manner of the early 
experiment by Rich (1951). 

Returning to the subject of fused quartz, we find that nearly pure fused 
quartz has J, about 51%, while both 96% quartz (Vycor) and 80% quartz 


/ 


(Pyrex) have 7. about two-thirds this amount. In both the latter cases the 
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added substance is mainly boric oxide. Here we have a case where J» varies 
much less violently than in the case of NaNO; in water, and which resembles, 
so far as it goes, the gas results of Dulit, Gittelman, and Deutsch (1956) 
mentioned earlier. It is therefore likely that the energetics of positron capture 
are responsible for these results. Similar remarks probably apply to the 
salts in Table I that show slight effects on J». 


We wish to thank Mrs. N. Hollbach and Dr. L. Yaffe for the chemical 
separation of the Na*® source material. Valuable suggestions have been 
received during the course of this work from Dr. J. S. Foster and Dr. J. D. 
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FULLY DEVELOPED TURBULENT FLOW IN A CURVED CHANNEL! 


r, 0,3 = cylindrical polar co-ordinates of reference point of channel, 
referred to center of curvature of channel as origin (see 


R,, Re = boundaries of channel (see Fig. 2). 

R = parameter such that Ri < R < Ro. 

U, V, W = components of mean flow velocity in 6, 7, 2 directions 
respectively (V = W = 0). 

u',v', w’ = components of fluctuating turbulent velocities in @, 7, 
directions respectively. 

ten’, = components of fluctuating turbulent vorticities in 6, r, 2 
directions respectively. 

P = mean pressure at reference point. 

p = density of fluid, assumed constant. 

Fy(r) = eddy diffusivity for vorticity. 

y loa s ; iE 

¢ =-—(rU) = z component of mean vorticity. 

ror 

CsriGe = parameters in vorticity distribution. 

v = kinematic viscosity of fluid. 

t = temperature at reference point. 

tri, tro, tr = temperature at r = R,,r = Ro, r = R. 

dy = k/sp = molecular thermal diffusivity of fluid. 

k = molecular thermal conductivity of fluid, assumed constant. 


STEADY STATE HEAT TRANSFER TO 


A. W. Marris 


ABSTRACT 


A vorticity transfer analogy theory of turbulent heat transfer is developed 
first for the case of fully developed turbulent flow under zero transverse pressure 
and temperature gradients such as that in the annulus between concentric 
cylinders rotating with different angular velocities or in a “free vortex”. The 
mean flow is assumed to be two-dimensional. The theory, which requires that 
the turbulence be statistically isotropic, yields a temperature distribution in agree- 
ment with experiment except in narrow regions immediately adjacent to the 
boundaries. An argument is given to show that the boundary layer thickness 
should be of the order of the reciprocal of the square root of the mean velocity, 
these boundaries are introduced, and Nusselt moduli are defined and their 
dependence on Reynolds and Prandtl numbers is investigated. 

The temperature distributions for the case of non-zero transverse temperature 
and pressure gradients, i.e. for the case of flow in a curved channel in which the 
fluid does not flow back into itself, are then obtained and the applicability of 
the simpler equations for zero transverse gradients to this case is investigated. 


SUMMARY OF NOMENCLATURE 


Fig. 2). 


‘Manuscript received November 5, 1956. 
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Ss = specific heat of fluid, assumed constant. 
Fy(r) = eddy diffusivity for heat transfer. 
1 oP 1 oP 
Pr, Ps ” | an Oe Ce ae 
2 
Pu Pe a an ee 


ee a. 
apjae0 Cip 90" apyaero Cop 90 
\ = Fy(r)/Fy(r) = analogy parameter. 


QO; = AP. 

QO» = AP». 

r = —rq,/Sp. 

qr = radial heat flow per unit area per second at r (see Fig. 2). 
ei ee = width of boundary layers at r = Rj, r = Ro. 


Nre = (Ro—R1) Un/v = Reynolds number for channel. 


1 0 = mean velocity in channel. 

Np, = v/dy = Su/k = Prandtl number. 

Q,’ = O1/ U5. 

Q.’ = 02/ Un. 

a1, Bi, Y1, 61, @2, Be, V2, 62, 62 = parameters in the radial temperature distribu- 
tion. 
EB 1 qr 

& ~ Un Ot/80 — SpUmdt/7.d0 

te = mean mixed fluid temperature. 


INTRODUCTION 


In this work heat transfer phenomena for turbulent flow in annuli and 
curved channels are investigated by postulating eddy diffusivities for heat 
analogous to eddy diffusivities for vorticity which latter have earlier been 
shown to be capable of yielding formulae for the mean velocity distribution 
which agree with experiment as long as the ratio of the radius of curvature of 
the channel to the channel width is relatively small, i.e. so that the flow does 
not approximate to straight channel flow. 

Whereas it has become apparent that fruitful insight into the ultimate 
micromechanism of turbulence, and particularly the effect of solid boundaries 
on this mechanism, is not likely to be obtained from such a phenomenological 
investigation, it is felt that such an approach must retain a value for the 
design physicist and engineer in that it represents the most tangible means of 
predicting with adequate accuracy the effect of the geometry and general 
flow conditions on these mean phenomena (Goldstein 1950; Klebanoff 1954; 
Laufer 1951). 

This then is the main justification for this work, the need for applicable 
formulae to yield temperature distribution and heat transfer data in the 
numerous problems of rotary motion that arise in the fields of engineering and 
in meteorology and oceanography. Moreover, while such semiempirical results 
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probably cannot be used as a starting point for any fundamental theory of 
turbulence (this because of the generally slow response of the macroscopic 
results of the turbulence to changes in the microscopic activity), it would seem 
that they must be at least an approximate end-point to such a theory, and as 
such possibly serve as a guide in its construction. 
BACKGROUND 

In cylindrical co-ordinates the Reynolds equations expressing the compo- 
nents of mean velocity in terms of the time average of the turbulent velocity 
and vorticity components take the form, for steady motion, 
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The equations of continuity for the mean and turbulent velocity components 
are respectively 
oU ,aV,aW ,V _ 
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0, 
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For two-dimensional mean flow in circles, 
(11) V=W =0 


and apart from the pressure gradient 0P/rd0 all derivatives with respect to 6, 
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and the mean derivatives with respect to z, vanish, so that equations (1) and 
(2) become 


1 oP an oz aU ] 0U U 0 1 0 
12 =. see ia —w { | —>;+- ——5 i fel ay! ’4 | camel” el seed | 
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(13) eta eg tee. 
where, by equation (10), 
(14) v't’—w'r’ = -4 Suv) 
and 
fqn Ree be 0 q” 0 72 vy” u” 
f oo" me =—_— *——_- Ye —_— 
(15) ai ii dr 2 or’ oe : 


It is an experimental result (Taylor 1937) that for flow in circles under 
zero transverse pressure gradient, for example for flow in the annulus between 
concentric rotating cylinders, the mean velocity distribution is given by 
(16) rU = const. 
over all the channel except for the narrow region immediately adjacent to the 
walls. Putting dP/rd@ = 0 in equation (12) above, one sees that this is only 


possible if 
= tut tae o 1 @ J 
(17) vo —wy = Fen] 2 “ wy], 


Fy(r) being any non-zero function of r. The argument for the vorticity transfer 
theory is that for statistically isotropic turbulence it predicts just this result 
(Goldstein 1938). Near the solid boundary the turbulence is known to be 
non-isotropic and the mean velocity distribution is no longer given by 
rU = const. 

While equation (17) may not represent the whole story of the effect of the 
turbulent micromechanism on the mean motion, its form, once Fy(r) is deter- 
mined, lends itself as a basis for producing a mean velocity distribution 
formula of use in design problems. Further one may develop from it, on an 
analogy basis, expressions for mean temperature distributions in such a 
curved channel. The general sluggishness of the response of such mean 
phenomena to the microscopic impulses, which tends to prohibit an investiga- 
tion of the latter from an Fy(r) determined from mean distribution, has 
perhaps the compensating advantage of giving stability and reproducibility 
of these mean distributions. 

THE EDDY DIFFUSIVITY FOR CURVED CHANNEL FLOW 


Assuming that the relation 


tk ee . 0 1 0 r 
(17) voi—wy = Fe(| 2 ae ww) | 
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holds, at least in the region not immediately adjacent to the solid boundaries 
for the case when 0P/rd0 is small but not zero, i.e. for flow in a curved channel 
with a depth to breadth ratio sufficiently large for the flow to be two- 
dimensional, one would have for the region of statistically isotropic turbulence 
lL @P 010 og 
18 - —— = (y+ Fy(r))| — =~ — (UV) | = t+ Fe(r)) =, 
(18) p 00 erry cca ) ( vir) 3, 
U being the mean transverse velocity at radius r (see Fig. 2) and 0P/0@ the 
mean pressure gradient in the flow direction. 
Differentiating equation (13) with respect to 6 one has 


(19) 0°P/drdd = 0, 


so that dP /06@ is independent of r. It follows from equation (18) that, assuming 
dP/00 to be also independent of @, 
i NOG A L<aP 
(20) (v+Fy(r)) oo , where.4 = — — = ‘constant. 
or r p 00 

Wattendorf (1935) obtained for such a flow a radial vorticity distribution 
for which the vorticity is not constant as in the case of Taylor’s Couette tvpe 
flow (equation (16)) but which is representable by the equations (Marris 
1956) 


(21) C=) = Cy/(r— Ry) for: Riv< << R 
and 
(22) 6-60 = 62-0 = Cy (Ro—r) for R Se ae Rs, 


yr = R, and R» being the boundaries of the region of turbulence under con- 
sideration, and R a parameter such that Ry < R < Ry. The parameter C; is 
positive and Cy, is negative. 

One obtains immediately from equation (20) 


23) Fy(r) tv = Fy,(r)+v = Pi(r—R)?/r_ for Ri <r < R 
and 
(24) Fy(r)+v = Fy,.(r)+v = Po(Re—1r)?/r for R <r < Ro, 
where 
1 oP 

D5 2 —— ete 
(< >) Py Cip 00 ’ 

1 oP 
% P = a cae 
26) Ps Cp 00” 


and since 0P/0@ is negative, P; and P» are both positive. 
By equations (14) and (17) 
97 0 


] > : 
(27) vg’ —w'y’ = aah = (r°u'v’) = Fy 


at(«dils PCa 
ah. eee ee 
Or pr 06 or 
Equation (27) integrates to give for the total shear stress (Eskinazi and 
Yeh 1956) 
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a1 ,9 (U)_ _1 oP fn 
(28) Pp =O (2) a i s-1), 
where 7, is the radius at which 7 = 0. 

Eskinazi and Yeh (1956) measured the distribution of both u’9’/U2,5 
and +/U%,x for flow in a curved channel of the type here considered, for which 
R, = 27in., Re = 30in., and Umax = 98.36 ft./sec. They find that 7, = 27.96 in. 
Their results are shown in Fig. 1. The discrepancy between the value of the 
































Fic. 1. Radial shear stress distribution for Eskinazi’s channel. Full line is total shear 
stress 7/U2 ax. Points are values of u’v’/U}nax. 
total stress 7 and the value of w’v’ agrees with that found by Laufer (1951) for 
the case of straight channel flow. 

From equations (21), (22), and (27), one has 


" ty! )P % 
(29) me eee (%-1)+2 


r 


i 2 00 pU ax 
where 
- a £ r'dr 
(30) A — i r- . n(?—Ri)” for R < r < R 
and 
_ _ C2 \ rdr 
(31) A — Fil - 5 (Ro—r)* for R < r < Ro. 


From Eskinazi’s mean velocity distribution (see below, and also Fig. 5) 


Ci/ Umax = 0.0357, Co/Umax = —0.044, and R = 28.75 in. 


one obtains 
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(approx.). The distribution of the correction, positive for r < r,, and negative 
for r > r,, was found to be such as to decrease the error; however with the 
values of C;/Umax and C2/Umax chosen, the magnitude of A is insufficient to 
account for more than a small fraction of the discrepancy. 


HEAT TRANSFER IN CURVED CHANNELS UNDER ZERO TRANSVERSE 
PRESSURE AND TEMPERATURE GRADIENTS 


When 0P/06 = 0 equation (18) gives 


(32) o+ Fri] 212 eu) |= o+rim X=0 


so that ¢ = const. irrespective of Fy(r). 

It is verified experimentally (Taylor 1937) that for such flow in the channel 
between concentric rotating cylinders, ¢ = 0, and that in a so-called free 
vortex, ¢ = 0, except in regions near the center of the vortex. One obtains 
(see equations (21) and (22)) C, = C. = 0 and ¢» = 0 so that the constants 
P; and P» of equations (25) and (26) are indeterminate, taking the form 0/0. 

To overcome this difficulty one will regard this case as the limiting case of 
channel flow and set 


: | .ef 
me) oe ke we 
and 
(34) — 1 <~@F 


it ee 
apjae0 Cop 06’ 
where P; and P» are positive and finite. The eddy diffusivities for vorticity 
will then be taken as given as before by equations (23) and (24). 

Referring to the polar element in Fig. 2 the equation analogous to equation 
(18) representing the steady state heat condition in a curved channel is 


(35) U - = II (ont Ful) 4 Fe - I (ont Fatr)) 2 4 
- 2 Gout Falr))? at, uct Fale) S. 


where Fy(r) is the radius-dependent eddy diffusivity for heat analogous to 
Fy(r) of equation (18). 
With 02/06? small, equation (35) becomes 


me a8 ’ at 
(36) U 06 aie - (Qg+ Fu(r)) r _ ° 


For rotary annular flow (i.e. such that the liquid flows back into itself), 
no heat is transferred in the flow direction but only radially through the fluid, 
and one has 


(37) 0t/d0 = 0 


and from equation (36) 


0 ; Ot 
(38) . (Qy+Fyu(r)) r “| = 0. 








on 
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¢.0d0 + $(9,de)dr 





Fic. 2, Channel parameters and polar element. 


It is seen that equation (38) does not involve the mean radial velocity distribu- 
tion. 

Integrating equation (38) one has 
(39) (aQgt+Fy(r)) dt/dr = T/r, 


where 7 is independent of r. Neglecting @g in comparison with Fy(r) one 
has 


(40) Ot/dr = T/Fy(r)r, 
oT ir 

4] t—tr = rf —., 

GH) n= T J, Far 


Asserting analogy between the eddy diffusivities for heat and vorticity one 
has from equations (23) and (24), neglecting vy in comparison with Fy(r), 
(42) Fu(r) = \Fy(r) = A\Pi(r—R,)? y= Qi(r—R)? 'y for R; CFS R 


and 
(43) Fu(r) = Qo(Re—r)?/r for R<r < Ro. 


\ being taken as a positive constant, Q; and Q» are positive constants. 


Continuity of Fy(r) at r = R yields 


1 Qi _ [aE 
nt Q:~ LR-Ri J’ 
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By equations (40), (42), and (43) one has 

a ae 7 | r—R | 

SoD les ——| 2 | forRi<r<R 
Qider (r—Ri) Q:(R—R,1) oe 





and 


r—Ri 
Tus _32ar r r—R 
(46) t—tr = a J (7) = 0.(R2—r) =) forR <r < Ro. 


For the particular case of constant eddy diffusivity for heat one has 








a z 
(47) t-te = 7 log eera 


where 7’ and Fy are constant. 
Using equation (44) to eliminate Q» from equation (46) one has 


(45a) -t( 9) = aR Fea for Ri<r<R 
ns 1 oo 1 
and 
_ 2) _ 1/(R=R,)’ ! af 
(46a ) (t ta) T = 2-2 R—r for R<r< Ro. 


T may be identified with the radial heat flux across the annulus, see equation 
(58) below, and there is evidence that Q; and Qs» depend linearly on U’, 
(Marris 1956), see below, section ‘‘Nusselt modulus’. It follows that there 
is only the one radius-dependent quantity, Q:i/U,, to be determined by 
experiment. 

Taylor's (1937) experimental results for flow under zero transverse tempera- 
ture and pressure gradients in the annulus between concentric rotating cylinders 
of radii R; = 1.865 cm. and Ry = 3.745 cm. are shown as curve 1 of Fig. 3. 
The mean velocity distribution is given by Ur = const. over the region 





38 
Phy 3.795 





Fic. 3. Radial temperature distribution for flow under zero transverse temperature and 
pressure gradients. 

Curve J}. Taylor's experimental distribution, 

Curve 2. Distribution calculated from equations (50) and (51). 

Curve 4. Distribution calculated from equation (52), case constant My 


n 
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(48) 1.995 <r < 3.595. 
Taking 
(49) R = 3.4 cm., tr = 30.05° C., T/Q: = —0.617, 


one has by equation (45a) 
3 3 0.617 r—3.400 3.400—r 
: == 26h Si i ea ee ee eo. Pie 
GS). (= = —5 ae | | ae E44 
for 1.865 <r < 3.400 


and by equation (46a) 
(51) t—30.05 = —0.103 2] for 3.40 < r < 3.795. 

3.795—r 
Curve 2 of Fig. 3 is calculated from equations (50) and (51). The calculated 
curve is seen to fit the experimental curve well except in a region immediately 
adjacent to the outside wall of the annulus. 

The constancy of Ur would seem to require that the turbulence be statisti- 
cally isotropic (a necessary condition for equation (18) to hold). One would 
also think statistical isotropy necessary for the analogy relative Fy/Fy = 
constant. It would therefore seem reasonable to expect equations (50) and 
(51) to fail outside the region Ur = const. 

Also on Fig. 3 is shown for reference the temperature distribution obtained 
from 


(52) t—tr = 2.44 log(r/3.4), 1.865 < r < 3.795 cm., 


i.e. from equation (47) with 7’, Fy = 2.44, for constant diffusivity. The dis- 
tribution is seen to be generally too flat and to fall well below the experimental 
distribution for almost a third of the channel width. 
NUSSELT MODULUS FOR FLOW IN CIRCULAR CHANNEL 
UNDER ZERO PRESSURE GRADIENT 
The Nusselt moduli for heat transfer for the case of the Couette flow be- 
tween concentric cylinders (or for ‘‘free’’ vortex) will be detined by 
ge,(Re—Ri) 
R(te, — try) 


Il 


53) Nu, ; 

», (Re— Ri) 
5. iii Qr, Ute 1 
(O4) Us ee he 


where gg, and qx, are the quantities of heat per second per unit area being 
transferred radially across the solid boundaries r = Ry, 7 = Re, respectively, 
tr, and fr, are the temperatures of these boundaries, and & ts the molecular 
thermal conductivity of the fluid. 

It must be noted that equations (53) and (54) are not the conventional 
definitions (Kreith 1955) of Nusselt moduli in that total temperature ditterence 
between the channel walls is employed rather than the diflerence between a 
wall temperature and the mean flow temperature. The above definitions are 
employed in order to maintain symmetry and simplicity in the analysis and 


at the same time produce expressions for experimentally measurable quantities 
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One has 
(55) -/Sp = —(dy+Fy(r))dt/ar, dy = k/sp, 
and equation (55) will hold for all r such that Ri < r < Re, when a suitable 
form is given to Fy(r). 


Also (by considering no heat accumulation in the polar element of Fig. 2) 


(56) pen ee 


00 sp or 


(Equation (56) of course follows from equations (36) and (55).) 
For the case being considered, d¢/06 = 0, so that from equation (56) 


(57) rg, = const. = Rigr, = Roger, = —spT (see equation (39)) 
and also 
(58) RiNuy = RoNusz. 


It will be assumed that the inner cylinder (of radius R) rotates with peri- 
pheral velocity U, (i.e. angular velocity U,/R,) and the outer cylinder (radius 
R2) is at rest. 

Except in the case of liquid metals the value of the Nusselt modulus is 
governed largely by the thickness of the layers of laminar flow adjacent to 
the boundaries. Buffer layers will not be considered in the present analysis. 
The flow channel is considered to be divided into four regions as follows 
(see Fig. 2) with the velocity and temperature distributions as given. 


(1) Inner Boundary Ri <r < Rite 


(59) U = U;R,/r 

(60) t = tri—qr,(r—R1)/k 
(2) Turbulent Region, Inner Section Rite. <r <R 

(61) U= D/r 

me t= nt GRR | a | 
(3) Turbulent Region, Outer Section R <r < Ro—€2 

(63) U=D/r 

(64) = et gaa LES | 
(4) Outer Boundary R2-e. <r < Ro 

(65) U = 0 (approx.) 

(66) t = tret+qr,(R2—1)/k 


One then has by equations (60), (62), (64), and (66) 


(67) try —trs = (try —tr, a+ (tr, +€1 —tr)+ (tr—trs_es) + (tr. —tro) 


- ee ao PE ane Beg qr: 
k at Oi (R—-R) €1 Q2(R.—R) €2 + k * 
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where (Rit+e:)—R and (R2—e€:)—R are written as Ri—R and R:—RK re- 
spectively. Using equation (57) one has 


oo a] qr * z Rages 
(68) tr, tre = [at R; k + Gn” Ce s 


5 Rea | qrs a! ae Rega 
a Ee R, k * 7 a Ss F 


By equations (59), (55), and (68) one obtains for the Nusselt moduli 








(69) ha = ls ines aH 
Ry | # 1 | 
| +B a2 Ga,” Un ay 
(70) Nu; = wesc ca? ; 
Rs | 2. J] 
=. €1 |+ Gas * Ces ays 
and 
(58) RiNu, = RoNuo. 
Equations (69) and (70) show the very great dependence of Nu; and Nuz 
on the boundary layer thicknesses €; and e. These must now be investigated. 
For case of flow through a channel with dP /0@ # 0, one would have 
(23) Fy, = —v+P,(r—R;)? rf, Rite < F< R, 
(24) Fy, = —v+P.2(R2.—r)? T, R Ge T< Ro— €, 
and 
lh. oF 
25) Ye ae 
( : Ps Cip 00 
1 oP 
oR i = — 
(26) P» Cp 00” 


where P,; and P» are both positive. One has 

. 2 

Fy, _ Pilr—Ri)” 
Vv 


(71) : : —1 forRkite <r<R, 
Pe D.(Re—r)* 
(72) Pry ” PAR | forR <1 < Ri—-e, 
or, near the boundaries, 
(73) Fv _ Pik ¢ ey 
v v r , 
ii Fy, _ PRs ( =) 
(74) oe l R: 5 


and it is seen that for regions of negative Fy, and Fy, to be excluded then 


> 14Vr/PiRi, ier > RitWV/Rw/P1, 
1 


r 


R 
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(76) =x 1—~/v/P2R2 , i.e. r<Ro- VRov/P2 . 


These boundary regions of negative Fy would seem to be incompatible with 
any physical conception of the turbulent mixing process, so that equations 
(73) and (74) are to be considered valid only outside this range. 

One therefore postulates boundary layers at least as thick as these regions 
of negative Fy and takes 


(77) a = VRw/Pi, 
(78) é = VRw/P2 3 


In further justification of equations (77) and (78) it may be remarked that 
in the case of rectilinear turbulent flow through a circular pipe, the relegation 
of the region for which the eddy diffusivity for momentum is negative into 
the boundary layer yields a boundary layer whose thickness seems to depend 
on Reynolds number in the experimentally required way (Marris 1954). 

Taking P; = \Qi, P2 = AQe, and assuming that the foregoing argument 
holds in the limiting case of zero pressure gradient, one writes 





(Tia) a= VRoPQ 
(778) = VEoTOs 
and so 
= = 
(79) ii seats ceil . 


(Gt ome) al 


or if \ is taken as unity 
(Ro—Ri) 


ra tt] d+) Ce 


where Np, = v/dy = Su/k is the Prandtl modulus. Also 
(58) Nue = (R;/R2)Nuy. 


(80) Nu; = — 





Equations (79) and (80) may be tentatively developed a little further. For 
flow past a rectilinear boundary it can be shown that the adjacent-to-boundary 
region in which the eddy diffusivity for vorticity becomes negative is 
OW N/(h/v)V(t0/p), Where N is constant in von KarmAn’s velocity distribu- 
tion law, / is the channel depth, 79 is the shear at the boundary, and for fully 
developed turbulence the shear velocity ~/70/p is proportional to U,. For 


this case then 
= 
_ on™ 7 


A similar result holds for the necessary boundary layer thickness when the 
vorticity transfer theory is applied to flow in a pipe of circular section (Marris 
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1954). These results can be identified with equations (77a) and (77) if Qi 
and Q» contain the mean flow velocity U,, as a factor. 

For the case of pressure flow in a channel 0P/06 is proportional to U,,? 
provided the Reynolds number is high enough. The fact that the velocity 
distribution obtained by integrating the vorticity distribution (from which 
C, and C; were derived) is capable of being fitted to experimental distributions 
for varying U,, when the velocity is expressed by the quotient U/U,, would 
indicate that C,; and C, contain U,, as a factor (Marris 1956). Thus it again 
follows that Q; and Q» are each proportional to Un. 

The final justification for believing Q; and Q2 contain U,» as a linear factor 
is that the kinematic viscosity »v would be expected to occur in equations (77) 
and (79) in the Reynolds number grouping. 


Writing 
(81) QO; = OF U a, 
(82) Q2 _ Q.’ Uns 


where Q,’ and Q,’ are functions of R; and R, alone, equation (80) now becomes 


4 
(#-1) Nne’ 
(83) Nu, = ———+ 


oD eGl al 


where Q,’ and Q,’ are independent of U,, and functions only of the radii. 
It follows that for a given channel 


+ Zz Nr ; Np ry Re 
84 Nu = =~, Nuz = = Nui, 
( ) ul (1+N m,* 2 R; 1 
where Z is a function of the radii, and that for low Prandtl number (liquid 
metals) Nui, Nus vary as Npe!N pr. 

TEMPERATURE DISTRIBUTION FOR TURBULENT HEAT TRANSFER FOR 
TWO-DIMENSIONAL PRESSURE FLOW IN A CHANNEL OF CONSTANT 
RADIUS OF CURVATURE 

From equations (21) and (22) for the radial vorticity distribution for two- 
dimensional curved channel flow, one has 





‘ er eee: Te 
(85) eo (rU) = wey for Ry <Gre<iR 
and 

. 1 0 ; 7 eet Cy 
(86) a (rU)—fo = ye" forR<r< R:. 


Integrating equations (85) and (86) one obtains for the mean radial velocity 


distribution across the channel 


(87) is eae forRi<r<R 


and 











424 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1957 


(88) U= re la iat ik ie 


Like the eddy diffusivities (equations (23) and (24)) these velocity distribu- 
tions have their origin in the experimentally observed vorticity distribution 
and must be regarded as empirical on this account. They were derived pre- 
viously (Marris 1956) on a mixing length basis analogous to that used by von 
Karman and Goldstein for flow in straight pipes and rectilinear channels, 
a way which showed that a reasonable assumption concerning the dependence 
of mixing length on channel curvature would give rise on Prandtl’s hypothesis 
to the eddy diffusivities and consequent vorticity distribution. One can say 
perhaps that the status of the mean velocity distributions is thus raised from 
empirical to semiempirical, and hope that sometime the mixing lengths 
referred to may be approximately deducible from fundamental considerations. 

The velocity distributions of equations (87) and (88), while they must fail 
in the immediate vicinity of the channel walls, giving velocities of minus 
infinity instead of zero, seem however capable of being fitted to experimentally 
observed velocity distributions for channels of relatively small radius of 
curvature, the fit breaking down at the inner wall of a 2 inch wide channel 
when the radius of curvature exceeds about 2 feet (Marris 1956). 

That the functional form of Fy allows such a fit for small radius of curvature, 
and the same functional form Fy gives a radial temperature distribution in 
agreement with experiment for the Couette channel of small radius, would 
seem to be justification for taking the analogy parameter \ = Fy/Fy as 
independent of r. 

Considering now the temperature distribution in the curved channel, one 
has by equation (36), neglecting the molecular diffusivity ay, 


ae ® at 
(89) U0 Trait) , _ 


and from equations (87) and (88) for the heat transfer case 


Ra log(r— Ra) Di) at _ - ; } 
(90) | co 14 F + og ee Fy(r) ts forRi <r<R, 


R2 log(R2—r) | , for , De at 0 ot 
(91) | — cf 1-RelonRamn | te ae = AT ryt) 7a] for R<r<Ro. 


r 


There seems to be experimental evidence that Fy, C;, D:, C2, Ds, &o« depend 
linearly on the mean flow velocity U, provided U,, is great enough to ensure 
fully developed turbulence (Marris 1956; Wattendorf 1935). Thus it is possible 
to fit equations (87) and (88) with experiment when they are divided through 
by U, to give similar expressions for U/U, (Wattendorf 1935). The fact 
that Fy and ¢ depend linearly on U,, is further indicated by equations (27) 
and (28) and from the experimental fact that a continuous experimental 
plot of u’v’/U,,? for varying U,, is obtainable (see equation (28)). Equations 
(90) and (91) can thus be written 


y 


~ W 


wa = 


aa ES S——<C 
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Ri leet) =) oe a = 
(92) | a 14 - + - Un a= Fy(r) — for Ri <r<R 


and 


Rs lo is) By ] a = ~ 
(93) | ad 1+ : + : +Eor | Un as Fu(r) r . 
for R <r < Ro, 


where 
4) 4-7, Bap, A--p, Bap, Baa. 


and A, Bi, As, Bs, E2 are all positive functions of the channel radii. 
From equations (92) and (93) one obtains 


oe a ae flat 2) | 421] 
(95) Fu(r) 13, = Um ae "4 be ‘ale dr forRi<r<R, 


(96) Falr) rS = Un an fp tetteemn D4 Bs |e 
for R<r< R:, 


taking 0t/dr as zero atr = R. 
Equations (95) and (96), like equations (87) and (88), fail near the bound- 
ariesr = Rj,r = Ro, because the logarithmic terms become very large negative 


?. ot / or 
causing 
ot 





to become negative. 


The constants A; and Az, in the velocity distribution are small compared 
with B,; and By (£2 is also small), and the velocity in the center region of the 
stream is very nearly given by V = k/r as in the case of flow in an annulus 
under zero pressure gradient. The effect of the terms A,R,log(r—R)/r and 
A2Rlog(R.—r)/r is negligible in the center region of the stream but they 
become large negative in the vicinity of the boundaries of the region, can- 
celling out the positive terms in the expression and thereby causing the rapid 
decrease in U at the boundaries. 

Whereas these terms could be integrated in series, it is considered that the 
boundary velocity distribution is too uncertain to warrant this arithmetic 
complication. Instead, the expressions on the right-hand sides of the 


inequalities 
: R) Ri 
( denies Pee 
(97) loe( 1 ; > ne 
r 
( ' See 
(98) log(1 ie) as ins 


will be used instead of the logarithm terms to represent the velocity decay at 
the boundaries. 
For R; <r < R, by equation (95), 
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B= 08 {Lote a8) se) 
Fa(r) ra = Un yr} "| act : + “ log\ 1 ; +log r} | dr 
(99) =Ur A a, ae ALR5( tog = + (Bit AiRi+-AiR; log R) log % 


—A,R, log 5 hs 





and for R <r < Ro, by equation (96), 


a ee’ Aske = )] 
Fu(r) ra = Um 39 [4 4,422 + Ew+ : log\ 1 RB +log R2} | dr 


(100) = Up, sais (B.+A>2R> log Re) log ‘ 
bes | 
+AoRs log ee e 
Setting 
(42) Fy(r) = AFy(r) = Qi(r—R,)2/r = OY Un(r—R)?/r_ for Ri <r < R 
and 


(43) Fi(r) = Q2(r—R2)?/r = Qe’ Um(r—R2)?/r for R << <, Rs 


and neglecting the small and slowly varying term }A,R,(log r/R)? in equation 
(99), one obtains from equations (99) and (100) the expressions for radial 
mean temperature gradient 


(101) Biel 5 ee 
Af 


j* i(r—R)+(Bit+AiRitAiRilogR)log (r/R) — Ai Rilog[ (r— Ri) /(R—- Ri] 


R (r—R,)° 
for: Ry. <4 < R 


and 


Lie» .. 1.(r—R)+3E2 u(r" — R’) 
(102) ol oo 1? S( (Ri—ry 


(B.+A,RelogR»)log(r/R)+A,Relog|(Rs—r) / (Rs 2) 

dr 
(R,—r)? 

for: R <* < Re: 


Evaluating the integrals /; and /., one obtains 


1 ot 
iy call 
(1 )3 ) t lr 0,’ 00 


By , R, R e r 2 Seth Sis aa 
(« ry 1) oe R RTT 1 Ry +8 A re Eee 


where 


al 
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a= A,(2+log R) + ’ 
} Bi = AiR, 
¥1 wes ‘aw S 
6, = A,(1+log Rr)+2 = a,—A}, 
Ri 
and 
ot 
(105) t—tp a 30 
4 aft) log R.- -R,, ry r,, (rR) 
7 (n+ gH ER PRT R Rr 7 Ee Mn 
forR<r< R:, 
where 
Bz 
ag = E; oRo+. 1, 2(1+log Ra) +p» 
Bz = A2Ro, 
‘ J Ax(2R.—R), 

10¢ eo SS ees Vo 9 
—_ \" 

do = Agel ‘ms = —A.—E2R 

eS eee ae eos w2it2, 


CORRELATION WITH EXPERIMENTAL DATA 

The mean radial velocity distributions given by equations (87) and (88) 
are fitted to the experimental distribution of Wattendorf (1935) for water 
flow in a channel for which R; = 20 cm., Rs = 25 cm., and R = 23 cm. 
(approximately), and to the recently obtained distribution of Eskinazi and 
Yeh (1956) for air flow in a channel with R; = 27 in., R. = 30 in., and R = 
28.75 in. (approximately). 

For Wattendorf’s channel one has for the ratio of the mean velocity at r 
to the mean flow velocity (Marris 1956) 


3.527 07 . 
| Lu OE logio(r — 20) + - for 20 <r < 23, 


(107) Pe 

U cs 34 9.20 ; 

U = ().0596+ ; Diet) +0.02373r+ oe for 23 <7 < 2, 
sO that A ; = Ci l = 0.0766, By = dD, v = 20.7, As = Cs l fo = 0.0596, 
B, {o/2U,, = 0.02373, and E: D,/U,, = 9.2. 


The plot of equation (107) together with Wattendort’s experimental results 


is shown in Fig. 4. 
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r (cm) r (inches) 
Fic. 4. Radial velocity distribution for Fic. 5. Radial velocity distribution for 
Wattendorf’s channel. Eskinazi’s channel. 
Curves 1 and 2. Plot of equation (107). Curves 1 and 2. Plot of equation (108). 
Curve 3. U/Um = 25/r. Points are Eskinazi’s experimental readings. 
Points are Wattendorf’s experimental 


readings. 


For Eskinazi’s channel one has for the ratio of the mean velocity at r to the 
maximum velocity 


U 7 2.220 ak LD 6 ae 
— 0.0357 + a logio(r —27) + for 27 <7 < 28:75, 


(108) ‘ cia 


3. 479 oe 
# aa 0.01404 9.088 logio(30—r) +0.032027 +92 for 28.75 <r<30, 
SO that Ci U mex = 0.0357, dD, Umax = 27, Cy Unex = — 0.044, Co 2U was = 


0.03202, and E2/Umax = 0.479. 

The plot of equation (108) together with Eskinazi’s experimental results is 
shown in Fig. 5. 

Comparing Figs. 4 and 5 it is seen that the present theory fails near the 
linear boundary when the radius of curvature becomes large. This was remarked 
on previously (Marris 1956), the same having been found to be true when the 
theory was correlated to Wattendorf’s experimental data for a channel of 
larger radius of curvature. 

Eskinazi’s distribution can be fitted in the region adjacent to the inner wall 
by taking a larger value of Ci/Umax, Say Ci/ Umax = 0.1416, but only at the 
expense of the fit in the center region of the channel. 

Owing to unavailability of data it is not possible to compare the temperature 
distribution given by equations (103) to (106) directly with experimental 
data. What will be done instead will be to employ the Wattendorf channel 
quoted above, for which the most satisfactory over-all fit is obtained for the 
mean velocity distribution, and determine the temperature distribution that 
would occur in it. 
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Using the above values of Ai, B,, Ax, Bo, and E, the radial temperature 
distribution given by equations (103) and (105) becomes, for this channel, 


(109) pink & 


Q,’ 30 
- (3. 290 + 3528) logio 7 4 sau (3) +30 1” r logio ~~ 2 for 20<r<23 
20 3 20 20 
and 
1 at 
(110) -te/ 35 ] 
25 
= (2 so24 3:18 logio =~ + 1.398 ra 4. 2920 5-"-- logie 53 
—0.0119 (28) for 23 <r < 25 
25—r 


The radial temperature distribution as calculated from equations (109) 
and (110) is shown in Fig. 6. The distribution fails near the outside boundary 

= R» owing to the velocity distribution (88) giving U = —@ atr= Rs». 
The error would seem to include all r such that 24 < r < 25. The extrapolated 
temperature distribution is shown by the dotted curve. 

















Fic. 6. Radial temperature distribution as calculated from equations (109) and (110), for 
non-zero transverse temperature and pressure gradients. Extrapolated curve near outer 
boundary is shown dotted. 
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Using equations (44), (45), and (46) for the case of flow under zero trans- 
verse temperature and pressure gradients, for this case also, one would write 





4 5 ee — | 
(111) t tn $3 eo” Rh Lik, forRi<r<R 
and 
i tae Me = | 
(112) t ie 07 0 R,—-RLR:—R for R <r < Ro, 
where 
sa sila apie a 
(113) g= U,, 31/86 ~ spU. at/r00 approx. (see equation (57)) 
and 
(81) Q1 = UnQ1’, 
(82) Qo = UnQz’. 
For the channel being considered one would then have 
aysiu — 
(114) t el 0, 06 ~ 3 Lr—20 for 20 < r < 23, 
‘ _ 4.3 ee 
(115) t te/ 3 a6 2 | 95-7 for 23 <r < 25. 


The temperature distributions given by equations (114) and (115) for g = 0.15, 
1.5, and 3.0 are shown by curves 2, 3, and 4 of Fig. 7. Curve 1 of Fig. 7 is the 
distribution as calculated from equations (109) and (110) and is given for 
reference. 

It is seen that in the Ri <r < R region t—tpz as given by equation (114) is 
increasingly less than t—tp as given by equation (109) as g decreases from 
about 2.25 and that in the outer region of the channel the values of t—tz as 
calculated from equation (115) are strikingly greater than the values as 
calculated from equation (110) for all g > 0.15. For g = 0.15, the values of 
t—tp as calculated from equation (115) and equation (110) are substantially 
the same until the point of failure of equation (110) is reached. 

The conclusion seems to be that the employment of equations (114) and 
(115) and the foregoing Nusselt modulus analysis to the case of non-zero 
transverse temperature and pressure gradients would only be possible if g 
was given two different values, a relatively high value in the inner region of 
the channel and a low value in the outer region. 


SUMMARY AND CONCLUSION 


A radial temperature distribution for Couette flow in the annulus between 
concentric rotating cylinders which agrees satisfactorily with experiment 
(Taylor 1937) except in the immediate neighborhood of the outer wall of the 
channel is obtained by using an eddy diffusivity for heat which is a constant 
times the eddy diffusivity for vorticity which latter yields the velocity distri- 
bution for flow in channels of relatively small radius of curvature. 
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Fic. 7. Radial temperature distribution for non-zero transverse temperature and pressure 
gradients using equations (114) and (115). ; ; 
Curve 1. Reference curve as calculated from equations (109) and (110). 


Curve 2. g = 0.15. 
Curve 3. g = 1.5. 
Curve 4. g = 3.0. 


The temperature distribution is given by the relatively simple formulae 


a Er r—R 
(45) t—trp = 0:(R—R:) [ =e) for RR: <r <R, 
(46) ~~. as (= forR <r < Rz, 


where ry = R is a circle dividing the annulus into two regions such that for 
R, <r < R the effect of the inner wall predominates and for R < r < R, the 
etfect of the outer wall predominates. For the case of flow in a curved channel 
(see below) for which 0¢/00 and 0P/06@ are not zero the point r = R is that 
at which dt/dr = 0. 
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The temperature distribution given by equations (45) and (46) is independent 
of the velocity distribution provided 0P/d6 = 0 but the eddy diffusivities 
for heat are taken as analogous to eddy diffusivities for vorticity under a 
condition of isotropic turbulence and therefore it is expected that the tempera- 
ture distribution could only be valid over the region of isotropic turbulence. 

The distributions (45) and (46) agree with experiment over a region some- 
what in excess of the region of Ur = constant. 

The annulus is then divided into three regions, laminar layer at inside wall, 
turbulent region, and laminar layer at outside wall, and Nusselt moduli for the 
channel of the forces, 


5 Ny, a Gar(Re— Ri) 
sii mn R(tr, —tr.) : 


: vay = $02(Re—Ri) 
” ae R(tr, —tr:) ; 


are defined and constructed on this basis that the laminar boundary layers 
are of width O(1/U,,)*. As in the case when the vorticity transfer theory is 
applied to flow in a pipe (Marris 1954) the laminar layers are shown to be of 
necessity at least this width in order that the eddy diffusivities do not become 
negative. 

For a given channel the Nusselt modulus of the form given above seems to 
vary simply as Nae?/(1+1/Np,) where Nae and Np», are respectively the 
Reynolds and Prandtl numbers. 

The case of turbulent flow through a channel for which 0P/0@ and dt/ 00 
are not zero is then considered. It is necessary however to assume that 0¢/ 06 
is independent of r as well as of 6; however there is evidence that this approxi- 
mation causes but little error (1-2%) away from the boundaries (Marris 
1955; Seban and Shimazaki 1951). 

The temperature distribution is calculated on the basis of the velocity 
distribution for the 0P/d0 = 0, dt/d0 = 0 case. It appears that such an 
adaptation would require different values for the parameter 


ata 8 spUm at/rae 
for the inner and outer regions of the channel, and on this account would seem 
to be unsatisfactory. 

One closes with a remark concerning the differences in heat transfer 
coefficients at boundaries convex and concave to the flow. 

It is known that for Nusselt moduli defined by 


» — dr(R2— Ri) 


117 N 
" ~~ R(tr:—tm) 

for heat flow outwards from the inner (convex) surface and by 
(118) Nu,’ = 226%: Ry) 


R(tr,—tm) 
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for heat flow inwards from the outer (concave) surface, then (Kreith 1955): 
(119) Nu,’ = Nu’. 


Consider two points A and B in the turbulent region at distance x from 
the inner and outer walls respectively, where 


(120) €1, €2 <NX < R-R,, R.—R, 
i.e. so that 

(121) an Rit+x, 

(122) rp = Ro-x. 


By equations (42) and (43) 


 . On" 
(123) Fy, = Fu(Ri+x) ere 
Qh” 
D4) ; ale bes ee 
(124) Fu, Fy (R2 x) Re~x’ 
so that 


On Fay baa Q2(Ritx) 
ia Fu, Qi(R2—*) 


Ro-x 


i) Ritx 
R:—R 


by equation (44). 
It follows that 


Fag. Ri ar 


fim ee i S83 
oo Fo,” Me tLm—n 


for all x satisfying the inequality (120). 
If 


" R~Ri — 
9 nee c 
(127) Rick > VR:/Ri, 


as seems to be the case for a channel of relatively small radius of curvature, 
(128) Fy, /Fa, > 1 
and Nu.’ would be expected to be greater than Nuy’. 
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DIFFUSION IN MULTICOMPONENT METALLIC SYSTEMS! 


J. S. KiRKALDY 


ABSTRACT 


Onsager’s phenomenological scheme for diffusion in multicomponent liquid 
systems is examined for suitability as a description of metallic interdiffusion. 
Subject to certain restrictions and approximations a set of non-linear differential 
equations is obtained which can be simply applied to important boundary 
conditions. A solution of the set for the system iron, carbon, silicon is shown to 
provide a good fit to Darken’s experimental results. 


1, INTRODUCTION 


The present purpose is to find a convenient and concise mathematical 
framework for the description of interdiffusion in multicomponent metallic 
solid solutions, for the correlation of available experimental results, and for the 
calculation of concentration distributions corresponding to simple boundary 
conditions. The subject of diffusion mechanisms will be introduced only 
insofar as it might assist in the simplification of the phenomenological equations. 


2. ONSAGER’S PHENOMENOLOGICAL SCHEME 


Onsager (1945-46) in his consideration of interdiffusion of liquid electrolytes 
assumed as a generalization of Fick’s first law (Fick 1855) that the diffusion 
current of each component is a linear function of all the concentration grad- 
ients. Accordingly, in one dimension, 


(1) = -D Dy Xe, (¢ = 1,2,...), 


where the diffusion coefficients D% are assumed to be functions of the C,. 
Some authors (Darken 1951; Bardeen and Herring 1951) have expressed the 
view that the expression 


(2) J,=->Ma cI (Gi = 1,2,...), 
i Ox 


where the yu, are the partial molal free energies, is more fundamental since 
the virtual forces for diffusion are represented as the free energy gradients. 
Equation (2) has been used successfully by Bardeen and Herring (1951) in 
the description of radioactive and chemical diffusion in pure metals and binary 
alloys. However, this form is not very useful for the solution of practical 
alloy problems since few of the thermodynamic data are available, nor are 
they easy to obtain. On the other hand, concentration gradients are relatively 
easy to measure, and it is usually the variation of concentration in diffusion- 
dependent alloy transformations which we desire. Accordingly, we will focus 
our attention on equations (1) and investigate methods for their simplification 


and manipulation. 
1Manuscript received December 20, 1956. 
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In all the discussions it is specified that distance is measured in numbers of 
lattice planes, that concentration is measured in atoms per unit cell (Birchenall 
et al. 1948), and that the total number of lattice sites is very closely conserved. 
This specification of concentration units is fundamentally different from the 
usual ones (mole fraction, g./cm.*, or weight per cent) since it implies that in a 
binary interstitial solution the solvent atoms have no concentration gradient 
even though the interstitial solute atoms are distributed non-uniformly. 
For practical purposes, at small solute concentrations, the four specifications 
are approximately equivalent. 

For simplicity it will be assumed that interstitial components diffuse 
interstitially and that substitutional components diffuse by a simple exchange 
mechanism, ignoring thereby the complications involved in a vacancy mech- 
anism and the associated Kirkendall effect. 

Consider now an n-component system with 7 components diffusing inter- 
stitially (numbered 1 to 7) and n—r components diffusing by an exchange 
mechanism. In this case equation (1) applies subject to the conditions that 


n 
0 and >> C; = constant. 


1 i=r+l 


bed, 
= 
I 


(3) 


i=r+ 
With these equations, (1) reduces to 
n—1 
OC, 
(4) Jpz-2 Des, 
k=1 x 


where Dy = Dy, 1<k <r, and Dy = Dy —Dy, r+1<k < n-1. One 
of the J; is of course dependent. In agreement with Darken (1951) we see that 
with a suitable choice of coordinate system we can always reduce the 
phenomenological equations for an n-component system involving n? diffusion 
coefficients to m—1 equations involving (m—1)? coefficients. 

Onsager (1945-46) demonstrated by invoking his hypothesis of ‘microscopic 
reversibility’ that for a liquid diffusing system subject to no bulk flow and no 
pressure gradients, the D% can be related to thermodynamic quantities by 
Ym Ou: ry Ou , 

(5) X ac, O* = x ac, 2" 

These subsidiary relations allow the further reduction of the number of 
independent coefficients from (n—1)? to n(m—1)/2. Darken (1951) has 
recommended the application of the same relationship to metallic systems. 
This is probably a valid procedure since the usual assumption of conservation 
of lattice sites is approximately equivalent to Onsager’s conditions on liquid 
diffusion. However, in systems which show porosity and lateral distortion 
associated with the Kirkendall effect, large pressure gradients have necessarily 
been present. Consequently, in such cases, some doubt is cast on the applica- 
bility of equation (5). 

In any case, for the present these equations are of little practical use since 
none of the required thermodynamic data are available and they are generally 
quite difficult to measure. We have therefore to treat a system of m components 
with »—1 differential equations involving (m—1)? diffusion coefficients. 
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While this might appear quite unmanageable for three components and out 
of reason for four, we shall see that in special cases the values of the D’s are 
such that major simplifications are possible and accordingly the equations 
can be manipulated and integrated with reasonable facility. 


3. MANIPULATION OF THE DIFFERENTIAL EQUATIONS 
If we subject equations (4) for an m-component system to the corresponding 
continuity equations 


(6) 14S = 0, 


we obtain the »—1 simultaneous diffusion equations in one dimension: 


5 7s) aC, 
@) at eax LO* ax 1° 
Following a procedure often used in binary diffusion (the Boltzmann—Matano 


analysis (Boltzmann 1894; Matano 1933)) we seek parametric solutions of 
the form 


(8) C,= CA), 3¥=x/vy1, 

which must satisfy the non-linear ordinary differential equations 
Adc, _ a aC.) 

(9) 3 dn = Ua Da I 


Solutions of these latter exist such that all dC;/d\ vanish as \ — +o. This 
implies, in view of \ = x//t, that appropriate initial conditions C;(x, 0) 
are single step functions at the origin. Solutions of (9) therefore apply to the 
semi-infinite diffusion couple for which the initial conditions are 


C,= Cf torx >OCandt=0 Ge. X—«), 


(10) : 
C,=C? forx <Oandt=0 (ie.\ = —@&), 


Equations (9) can be integrated once to obtain 

Tan n—1 
(11) J MC; = —-2>> pa. 

ci? k=1 dx 

Since these equations relate the diffusion coefficients to the properties of 

the measured experimental curves, we can use them for the experimental 
determination of the coefficients. In binary alloys the procedure is a standard 
one and is relatively simple to carry out since there is only one coefficient to be 
determined from a single curve. For higher order alloys it is apparent that 
the extent of the experimentation and computation would be prohibitive 
if no approximations were made. In most cases suitable approximations 
should be made apparent by the experimental results. 


4, ANALYSIS OF EXPERIMENTS ON THE IRON, SILICON, CARBON SYSTEM 


Darken (1949) has performed an interesting series of experiments on 
carbon diffusion in y-iron in the presence of other soluble alloying elements. 
In one example he prepared a diffusion couple in which the carbon (component 
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1) was in approximately equal concentrations on both sides of the weld while 
the silicon (component 2) was present essentially on one side of the couple 
only. Fig. 1 indicates the initial concentration distributions as well as the 
observed carbon distribution (points) after a 13-day anneal at 1050° C. 


0.6 








SILICON 3.89% SILICON 0.05% 


% CARBON 





0.8 0.6 0.4 0.2 O 0.2 0.4 0.6 0.8 
DISTANCE FROM WELD (INCHES) 


Fic. 1. Interdiffusion in the iron, silicon, carbon system. 


To facilitate solution of the condensation of equations (9) applicable to 
this couple, the following assumptions have been made: 
(a) Dy, Dox, and Dy. are constant, 


dC, TC: 
(b) Do dn. K Doe dd? ° 


The final justification of these is that the solution of the differential equation 
based thereon can be fit quite precisely to the experimental curve. The validity 
of relation (0) is least certain, since for these particular initial conditions the 
calculated carbon distribution is rather insensitive to assumptions about the 
relative magnitude of the terms. In the absence of a measurement of the 
silicon distribution, the relation cannot be confirmed so it remains essentially 
a mathematical convenience. Should it be invalid, and the terms more nearly 
equal, the value of Dy» determined by fitting the curve will be incorrect but 
D,, and Dy, will be relatively unaffected. 
Applying (a) and (0) to (9), the equations become 


; A dC; Fe, rc. 
2 eS Se | ae Bare em, oe 
(12) 2 dX Du dn” +Di2 dn 
and 
3 A dC2 _ dC: 
(13) 3 an = OF De De 
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Equation (13) has the independent solution 


co —C2 
14 = C+ il ek 
(14) C, = C2 + i~eri = -/Ba 
which can be substituted in (12). The solution of equation (12) subject to the 
approximate conditions stated is 
A. 0 1 
(15) Cc, = CY. +Es -+-4 [ A 24 —-Ci +O, orf — A f 


erf =--"-— 


2 V Doe 2A ” 2 VD 


where 
(16) A = - aes eS ee : 


The coefficient A (and therefore Dj.) can in general be obtained by fitting 
equation (15) to the results. In the present case, since Do. K Dy, it is sufficient 
to read A directly from the curves by extrapolation to the origin. For the 
calculation of the theoretical curve shown in Fig. 1 the following values have 
been used in equation (15): C,° = 0.441%, Ci! = 0.478%, C.'—C.® = 3.84%, 
A = 0.135%, Dir = 4.8X 10-7 cm.?/sec. (deduced by trial fits to the experi- 
mental curve at points well removed eon the origin and in agreement with 
Darken’s values from the same data (Darken 1949)), Do. = 2.31079 
cm.*/sec. (a very rough value obtained by trial fits near the origin), and 
D2 = 0.34 10- 7 cm.*/sec. (calculated from equation (16)).* As we see the 
final fit is quite good and information has been obtained from this one couple 
about all three coefficients. 

The value of Ds. used in the calculations shows only order of magnitude 
agreement with data given by Smithells (1955) for the diffusion of silicon in 
iron but this is not surprising in view of the meager data for the critical 
region near the origin and the rather doubtful assumption (0). 

The value obtained for the off-diagonal coefficient D,,. indicates a surprisingly 
large interaction between carbon and a silicon gradient. The significance 
of such a result can best be made apparent by establishing the relationships 
between the D;, and the My, through equation (2). 


5. SUMMARY 
The very satisfactory description of interdiffusion in the iron, silicon, 
carbon system as evidenced by Fig. 1 suggests that equations (1) or (9) may 
be applied with equal directness and facility to other systems and to the solu- 
tion of various metallurgical problems with simple boundary conditions. 
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SOME SCATTERING PROBLEMS IN CONDUCTION THEORY! 


P. G. KLEMENS? 


ABSTRACT 


An expression is derived for the scattering of electrons by inhomogeneous strain 
fields in terms of the scattering of electrons by low frequency lattice vibrations, 
as deduced from the low temperature intrinsic thermal resistivity and the 
lattice thermal conductivity. This is applied to the scattering of electrons by 
dislocations and stacking faults. An expression for the scattering of lattice waves 
by stacking faults is also derived. The validity of these results is discussed. 


INTRODUCTION 


In order to interpret changes in the conduction properties of metals at low 
temperatures due to cold work or irradiation, it is of interest to calculate the 
scattering by various model imperfections both of electrons and of phonons. 
The former is related to changes in the residual electrical resistance, the latter 
to changes in the lattice thermal conductivity. 

While the theory of the scattering of electrons by vacancies and inter- 
stitials (Jongenburger 1953, 1955; Blatt 1955; Overhauser and Gorman 1956) 
seems to yield fairly reliable results, the scattering of electrons by extended 
imperfections seems not as well understood. Hunter and Nabarro (1953) 
have obtained an expression for the scattering of electrons by static strain 
fields, and applied this to the scattering by the strain field about a dislocation, 
deriving the magnitude of the interaction between electrons and the Fourier 
components of the strain field from first principles. An alternative treatment 
is developed here, which is formally similar to Hunter and Nabarro’s, but in 
which the interaction is deduced from the magnitude of the electron-phonon 
coupling as derived from the observed conduction properties. 

Care must be taken in interpreting the observed conduction properties, and 
the most reliable procedure seems to be that of taking the low temperature 
lattice thermal conductivity, limited by the scattering of phonons by electrons, 
as a measure of the electron—phonon interaction, since this quantity is not 
sensitive to details of the electronic band structure. Another reliable measure 
is the intrinsic electronic thermal conductivity at low temperatures. 

The expression for the scattering of electrons by strain fields thus derived is 
similar to that of Hunter and Nabarro, and is applied to the scattering of 
electrons by dislocations and stacking faults. 

An expression for the scattering of phonons by a static strain field was 
derived earlier (Klemens 1955), and applied to the cases of dislocations and 
grain boundaries. An expression for the scattering of phonons by stacking 
faults is now similarly obtained. 


‘Manuscript received December 11, 1956. 
Contribution from the Division of Pure Physics, National Research Council, Ottawa, 
Canada. 

Issued as N.R.C. No. 4229. 

2On leave from Division of Physics, C.S.I.R.O., National Standards Laboratory, Sydney, 
Australia. 
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The present treatment is more suited to the case of the scattering of phonons 
than to the scattering of electrons, because the wave-lengths of the latter 
are rather short, and most of the scattering arises from regions in the im- 
mediate vicinity of the imperfection core, when the elastic continuum theory 
cannot be expected to apply. 

SCATTERING OF ELECTRONS BY STRAIN FIELDS 

As pointed out by Hunter and Nabarro (1953), electrons are scattered 
only in an inhomogeneous strain field, which destroys the periodicity of the 
lattice, and it is the strain, rather than the displacement, which provides an 
appropriate description of the deviations from periodicity. Lattice vibrations 
also produce an inhomogeneous strain field, so that the observed scattering 
of electrons by lattice waves gives a measure of the scattering of electrons by 
the various Fourier components of the static strain field, provided the fre- 
quencies of the lattice vibrations are sufficiently low to allow the conduction 
electrons to adjust themselves fully to the ionic displacement, as they must 

the case of static strain. One can therefore use the interaction between 
electrons and low frequency phonons as a measure of the interaction between 
electrons and the components of a static strain field. 

The displacement u(x) of the lattice point x due to lattice vibrations is 


given b 
] 
} u(x) = - >, a(q) e, e* 
VG qj 
here e, is the polarization vector of the polarization branch j, q is the wave- 
vector, w is the frequency, G is the number of unit cells in the erystal, each of 
nass ‘J, and a(q) is the harmonic oscillator matrix, whose elements are 
‘q | h | 
a Ve N—I , 
IM 
, L2Mw - 
“yh 
h F 
¥ /( re ( 
o'er Ee nea 
Phe rate of chanye of the electron distribution function { due to interaction 
, viven b 
df | 2 - Ly 
3 f Ci? QAE/)(fO—fy(N+1 (1—f’) fA 
: LC Ass )| ) ) LIN], 
‘ at 4 WGh 4 , Ww | 
being due to all processes in which an electron kK and a phonon q combine to 
form an electron k’ or vice versa, Conserving Wave-vector 
k+q =k’, 
tyr eol the re SMOhialicd fine tor 
WAI 9 J cos(Alu ‘f) 
(Alt /h) 
here LJ Int ha It’, conserving energy. The constant C, is the electron 


ONO eraction parameter us defined, for example, by Sommerteld and 
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Bethe (1933), but generalized to describe the interaction between electrons 
and phonons of polarization 7. 
Longitudinal waves cause an inhomogeneous dilatation 


(4) div u(x) = - __ > ga(q) e** 
Vv G q 


whose effect is comparable with that of an inhomogeneous static dilatation 
A(x), resolved into Fourier components 


s ‘ 1 
6) At) = 3 Dig) e™* 
VG‘ 

so that one can replace ‘a(q)q in (3) by D(q) to obtain the rate of change of f 
due to the scattering of electrons by the strain field A(x). Noting the relation 
(2) between N and a’, one obtains 

df 4 2 72 ewes as 
(6) ‘if. 3D, Cy D*(q) Q(AE)[f’ —f] 

dt 9Gh? 2 ae ae enn oes 

where now AE = E—F’, and C, is the electron—phonon coupling constant 


for longitudinal waves. 
Transverse lattice waves similarly cause an inhomogeneous rotation 


= t 
(7) curl u(x) = —3 > (qXe) a(q) e™* 
VG“‘G 


which can be compared with the components of rotation of an inhomogeneous 
strain field 


. l 
(S) curl u(x) I(x) = > P(q) e“* 
qa 


VG 
so that in (3) one can replace ra(q)g for a transverse wave by P(q) orientated 
in the direction of (qXe). An expression, similar to (6), can be derived for 
the scattering by an inhomogeneous rotation field in terms of C),, the coup! 
constant for transverse vibrations, and ??(q). If it is assumed that there is 
no systematic interference between scattering by the dilatation and rotation 


telds, one obtains tor the scattering by an inhomogeneous strain field 


dy t “ » » > > 
0) t —= » {C, D'(e)+C,, P'(e)} Q(AE 
dt 0GK <= *™* \q at . 
Since the interaction parameters Ce and Cy? depend sensitively on the 


departure ot the electron wave-tunctions trom tree-electron wave-functions 
the estimation of these parameters trom first principles is a tormidable task 
and we shall rather estimate them trom the observed conduction properties 
Che electronic conduction properties lead to an estimate of C CY+23c 
However, the high temperature resistance does not lend itse 
estimate ob OC?) because the ettects of Umklapp- processes and the dispersio 
ot the lattice waves are not easy to take into account. Lt these effects were 


disreparded, one would arrive at an overestimtate of GC?) Phe electrical resistances 


*See, for example, Hardeen (ia0) 
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at low temperatures is likewise unsuitable, because it is very sensitive to the 
shape of the Fermi surface, as discussed elsewhere (Klemens 1954). However, 
the thermal resistance at low temperatures is not sensitive to the shape of the 
Fermi surface, and can thus be used to reliably estimate C?. In Table I values 
of (C/¢)? are given, estimated from the high temperature resistance (electrical 
or thermal) and the low temperature thermal resistivity. The latter leads to 
lower values. 

The scattering of phonons by the conduction electrons limits the lattice 
component of thermal conduction at low temperature, which is thus an 
additional measure of the electron-phonon interaction. The lattice thermal 
conductivity varies inversely as either C;?? or C,;°, whichever is the bigger, 
because the two polarization branches are closely coupled to each other by 
phonon-phonon interactions (Klemens 1954a, 1956a). Comparing the mag- 
nitude of the electronic thermal resistivity—which yields C*—to the lattice 
thermal conductivity, it has been concluded that 


(10) Ch = Cy = $C* 


in the cases of silver (Kemp, Klemens, Sreedhar, and White 1954, 1956) and 
copper (White and Woods 1954, 1955). In other cases no information is 
available yet on the lattice thermal conductivity, though it appears likely 
from the effective ratio of the electrical to thermal resistivities at low tempera- 
tures that the electrons interact equally strongly with waves of all polarizations 
also in the cases of gold and sodium (Klemens 19540). 
TABLE I 
ESTIMATES OF C? 


Element (C/f) 2 (C/f) 2? (C1/f)? 
Cu 1.23 0.20 0.07 
Ag 1.46 0.28 0.09 
Au 1.44 0.25 — 
Na 1.08 0.20 a 


(C/¢):? is the estimate from high temperature resistance data, 
(C/¢)2? from the intrinsic thermal resistance at low tempera- 


tures, and 
(C1/¢)? from the lattice component of thermal conductivity. 


SCATTERING OF ELECTRONS BY DISLOCATIONS 


The components of the strain field about a screw dislocation, whose axis 
lies along the z-axis, are given by 


(11) Tl, = oe cos 8, Il, = 3 sin 6 


where } is the Burgers vector, and (7, 8) are cylindrical coordinates. The 
strain field of an edge dislocation is 
6 1-—2y. 


A= co . sin 6, 


(12) 


II, = 3 cos 6 


where v is the Poisson ratio (Cottrell 1953). 
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The expression (9) for the scattering of electrons by a static strain field is 
formally similar to an expression for the scattering of phonons by a static 
strain field (Klemens 1955), which had been applied to the scattering by 
dislocations, so that some of the results derived there can be used for the 
present purposes. Thus the Fourier transforms of the strain field are 


G B 
(13) F (q) = E 3 a‘g a 
for screw dislocations, while for edge eae they are 
; G Bb ; Z a G Bb 
Ea = 8 “a oe = aa. = ~ 9 
(14) f (q) Ten a‘g” ’ D (q) 2Gy° a‘g 


where q is perpendicular to the z-axis (which implies that cial by a 
dislocation cannot alter the z-component of the electron wave-vector), Go 
is the number of unit cells of the crystal in a plane perpendicular to the z-axis, 
and a is the lattice constant. Noting that, if » is the direction cosine of k 
relative to the field, 


(15) E 4a) 2 [1- a | 2 I O(AE) = 2G, po(e5E) 


t 


in analogy with a similar summation discussed previously (Klemens 1955), 
and that this should be multiplied by a factor 0.55 if the dislocations are 
arranged at random, one obtains for the inverse relaxation time 


i at) eit Cy 


7 ~ P—fatl~ "9 Goa! k(kdE/dk); 


in the case of a screw dislocation. The scattering due to an edge dislocation is 
of similar magnitude, with a factor 


Lore =e)er| 


in place of C,,;°. For free electrons (k dE/dk); = 2¢, and since (a?Go)“! = N, 
the number of dislocation lines per unit area, one obtains 


1 acne 2a hey rz2 


(16) 


(17) - = V.009 “9 he NO. 
For copper ¢ = 7.1 ev. and 6 = 2.5X10-8 cm., so that 
(18) 1/r = 0.027N sec.“! 


and the increase in electrical resistance at room temperatures due to disloca- 
tions is 


(19) Ap/p = 0.7X10-!N. 


The values for the inverse relaxation time and electrical resistance due 
to dislocations derived here are somewhat smaller* than those derived by 
Hunter and Nabarro (1953) on the basis of a similar calculation, but using the 


*The results of these calculations were presented at the Conference on Electron Transport 
in Metals and Solids, Ottawa, (Klemens 19565). A numerical error was incurred in going from 
(17) to (18) and it was erroneously stated that the present considerations lead to a somewhat 
higher scattering probability than those of Hunter and Nabarro. 
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deformation potential to estimate the interaction between electrons and the 
strain field. This is to be expected, since the deformation potential yields values 
of C of the order of ¢—see Jones (1956). 
SCATTERING OF ELECTRONS BY STACKING FAULTS 
The present formalism can be used to calculate the scattering of electrons 
by stacking faults. A face-centered structure can be regarded as composed 
of closely-packed spheres of radius 7, arranged in hexagonal layers, the distance 


between layers being 


20) D. = 2r(2/3)'* 
and 
(21) y= QM8(aq/2) 


where a is the volume per atom. The layers are stacked one above the other, 
each being displaced in its plane by a distance 
(22) S = 2r/V/3 = a 26 3-12, 


successive displacements forming an equilateral triangle. 

A stacking fault can be regarded as the displacement of one half of the 
crystal through a distance S in the layer plane in one of the three possible di- 
rections, and thus corresponds to a shear S/D localized within a distance D. 


Then 


a 


(23) P(q) = YF (a) 52) 


where G; is the number of layers in a crystal containing one stacking fault. 
Thus P(q) vanishes unless q lies perpendicular to the plane of the stacking 
fault, which implies that the component of the electron wave-vector k in 
the plane of the stacking fault is not changed. Thus an electron can only 
suffer specular reflection at a stacking fault. 


Also 
zy 
2 Q(AE) = G; DiS) — 
ef) d ae) ee (44 lcos 0} wm 


where @ is the angle of the incident electron with the normal to the stacking 


fault plane. Taking an average value for 


1 | 
) es = eee. pa 
acid lcos 6] E be 


one obtains, on substituting into (9), 


1 _ 16 Cy'1_S" 


Wit 





26) -=——}- 

- 7 27 vGD 

while in the case of perpendicular incidence, the reflection probability is 
4Cyr S 

27) i= == 

o 9% v 


where v is the velocity of an electron of Fermi energy. In the case of copper, 
R= 0.4. 
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SCATTERING OF PHONONS BY STACKING FAULTS 
It is possible to use the strain-field representation (23) of a stacking fault 
to calculate the scattering of phonons. It has been pointed out (Klemens 1955) 
that a static strain field scatters lattice waves because of the anharmonicity 
of the interatomic forces, and that, in the case of a non-dilatational strain 
field, 
1 a’y’ J 2 2 = a] 
2 == ———4 wT Fr *—7 1-—— 


T 





where y is the Grueneisen anharmonicity parameter, and the integration is 
taken over the three surfaces (one for each polarization j’) w’ = w in k-space. 
With the value (23) for P(q), and using (25), the effective relaxation time is 
given by 
(29) 1 os a yw ss. 

Tv 18 G3c 3 
where c is the velocity of sound. Here again, phonons are specularly reflected, 
but the reflection probability, in the case of long waves, varies as the square 
of the frequency, so that the corresponding thermal resistance varies inversely 
as the absolute temperature. Th»: reflection probability for perpendicular 
incidence is 
(30) R = a@ky2/18 
where & is the phonon wave-number. 

DISCUSSION 

The expression for the scattering of electrons by inhomogeneous strain 
fields developed here is only valid for slowly varying strain fields, where the 
important Fourier components are of low wave-number, comparable to the 
wave-numbers of the phonons important in determining the low temperature 
thermal resistivities, from which C? has been determined. In general the 
effective value of C? changes as the wave-number of phonons is increased 
(Jones 1956). This is largely due to the inability of the conduction eleetrons 
to adjust themselves as well to high frequency vibrations as to low frequency 
vibrations, an effect which does not arise in the case of static strain fields, so 
that the use of the present formalism for rapidly varying strain fields does 
probably not introduce a large error. More serious, however, is the objection 
that the electrons of Fermi energy have a wave-length of the order of the 
interatomic distance, so that much of the contribution to the scattering 
Hamiltonian takes place in regions where the imperfection cannot be ade- 
quately described in terms of the strain field alone. This is an objection which 
does not arise in the case of the scattering of long lattice waves. 

Thus in the case of dislocations, for example, the electrons are scattered in 
the immediate vicinity of the dislocation core, where the displacement of the 
atom is not adequately described by the continuum theory. This objection 
applies to the present calculations as well as to those of Hunter and Nabarro 
(1953). A proper theory of the scattering of electrons by dislocations would 
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have to take account of the detailed behavior of the conduction electrons at 
the dislocation core. An attempt along such lines has recently been made by 
Seeger and Stehle (1956). 

When comparing the electrical resistivity increase on plastic deformation 
with the electrical resistance increase to be expected from (19), or from the 
corresponding relation obtained by Hunter and Nabarro, the difficulty arises 
that one is not sure how much of the observed resistance is due to dislocations, 
and how much due to point imperfections (vacancies, interstitials, or clusters 
of these). This difficulty has been exhaustively discussed by van Bueren 
(1955, 1956) and also by Seeger (1956a). It appears that in some cases 
most of the electrical resistance can be ascribed to point imperfections, though 
this does not necessarily imply that the resistance due to extended imperfec- 
tions is negligible, as in the case of brasses, where lattice thermal conductivity 
and electrical resistance were studied simultaneously (Kemp, Klemens, 
Tainsh, and White 1957). In other cases a substantial fraction of the resistance 
seems to be due to dislocations. It would then appear that there is a sub- 
stantial discrepancy between the observed dislocation resistance and that 
calculated by Hunter and Nabarro, as discussed by van Bueren (1956), and 
an even greater discrepancy with the present calculations. One way out of 
this difficulty would be to suppose that a large fraction of the dislocations 
are split into partial dislocations, joined by stacking faults, and that the 
stacking faults contribute most of the electrical resistance of this arrangement. 
This has been suggested, among others, by Seeger and Stehle (1956). 

If stacking faults are to contribute an appreciable fraction of the resistance 
increase due to plastic deformation, the reflection probability for a conduction 
electron impinging on the stacking fault must be somewhere of the same order 
as the value calculated here (equation 27), as discussed by Seeger and Stehle 
(1956) and Seeger (1956). 

The electrical resistivity due to stacking faults was estimated by Klemens 
(1953) by means of a wave-fitting calculation, assuming the electron wave 
functions to be strongly modulated by the periodic lattice. A reflection 
probability of 0.5 was obtained in this way. A more general wave-fitting 
calculation was carried out by Blatt, Ham, and Koehler (1956). As would 
be expected, the reflection probability depends intimately on the deviation 
of the electronic band structure from the free electron model: in the case of 
perfectly free electrons a stacking fault cannot influence the electronic motion, 
because it has no dilatation. Using a Wigner-Seitz approximation, these 
authors get a reflection probability which is considerably smaller than the 
value obtained by Klemens for strong modulation. 

It would appear at first sight that the present method of calculating the 
reflection coefficient in terms of the shear strain field of a stacking fault is 
not related to the above considerations. However it should be remembered 
that for free electrons (unmodulated wave functions) C? = 0, so that C,,’, 
and hence the reflection probability, depends upon the deviation of the 
electronic band structure from the free electron model. In the absence of 
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modulation of the electronic wave function C;;2 = 0 and a stacking fault 
would not scatter electrons. 

The fact that the present calculations lead to a reflection probability of 
the same order of magnitude as that calculated by a wave-fitting method, 
assuming strongly modulated wave-functions, gives some confirmation that 
the present method yields scattering probabilities of the correct order of 
magnitude, and confirms that, at least in the noble metals, the conduction 
electrons are far from free, as is also apparent from the various intrinsic 
conductien properties (Klemens 1954), 1956a). On the other hand one should 
not rely on these estimates to more than order of magnitude: firstly one can 
question the replacement of a stacking fault by a shear over one atomic plane, 
and more seriously, the applicability of perturbation theory in the present 
case, when the reflection probability approaches unity, is doubtful. Neither 
of these effects is likely to change the order of magnitude of the result, but 
probably falsifies it to some extent. 

Another approach still to the problem of the scattering of stacking faults 
has been made by Seeger (1956b), who treated the diffraction of the electron 
waves due to the displacement of one half of the crystal relative to the other. 
His result also yields reflection of comparable magnitude. 

No reliable experimental results are yet to hand on the effect of stacking 
faults on the lattice thermal conductivity. The frequency dependence of (29) 
must be regarded as well established, as it can be derived from quite general 
arguments. The absolute magnitude of (29) is more questionable, because of 
the approximations involved, the most serious being the use of the cubic 
anharmonicities for the case of a very strong shear, where this approximation 
could break down. Again the present result is probably an estimate of the 
order of magnitude only. One uncertainty, which is also of importance in the 
problem of scattering of lattice waves by dislocations, is the correct choice 
of the Grueneisen constant y, as discussed by Barron (1956). 


The work described here was carried out while the author was on a visiting 
appointment to the Division of Pure Physics, N.R.C., Ottawa. The author 
wishes to thank the National Research Council for their hospitality, and 
Professor F. R. N. Nabarro for helpful discussions. 
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A STUDY OF NUCLEON FORCES WITH REPULSIVE CORES 
III. LOW ENERGY PROPERTIES OF THE LEVY POTENTIAL! 


M. A. PRESTON? AND J. SHAPIRO® 


ABSTRACT 


An attempt has been made to select the core radius and coupling constant of 
the Lévy potential for the interaction of two nucleons in order to fit the binding 
energy of the deuteron and the singlet state neutron—proton scattering length. 
It was found that these two quantities cannot be fitted simultaneously. For any 
given choice of coupling constant, a somewhat larger core radius is required to 
fit the deuteron binding energy than is required for the scattering length. This 
spin dependence of the core radius does not preclude the possibility of a fit to the 
low energy data with the Lévy potential. 


An approximate two-nucleon potential has been obtained by Lévy (1952) 
from symmetrical pseudoscalar meson theory with pseudoscalar coupling by 
means of the Tamm-Dancoff method. Lévy was able to show that strong 
repulsions would occur in the relativistic region (i.e. at small separation 
distances). These were taken into account by using a hard core cutoff. For 
even parity states (with which we shall be concerned here) the potential is 

Vir) = @ {on tr’ <%,, 


1) 
( Vir) = V(r) +Vr(r)Si2 for r > r-, 


where r, is the core radius, and 

9 ae 2 ts) e” «(s) 1 S26 
1 °\ 
m2 Kun | f |. 


2( m, \" 3 ee 
3) vote) = — ance (ay Bie Se. 
( r(r) Ml OM +t) on 
M is the reduced mass of the two-nucleon system, A, the Bessel function of 
imaginary argument, Sj. the usual tensor operator, and yu is connected with 
the -meson mass by m, = u/c. Terms up to the fourth order in the pseudo- 
scalar coupling constant G are included (A? = G?/4rhc). The potential has two 
parameters (r,, \2) to be determined from experiment. Lévy obtained these by 
fitting the \S NP scattering length a,, and the binding energy of the deuteron, 
e. Using rather approximate numerical methods he found the result r. = 0.38 
+0.03 (in units #/m,c, which will be used throughout this paper unless 
explicitly stated otherwise) and \? = 9.7+1.3. With these parameters, the 
‘Manuscript received December 6, 1956. oan f 
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potential gave good agreement with the low energy properties of the NP 
sysiem, except for the deuteron quadrupole moment, which was 20% too low. 
Lévy also obtained reasonable agreement with the NP differential cross section 
at 40 Mev. Martin and Verlet (1953) used the Lévy potential to compute the 
PP differential cross section at 18 and 32 Mev., obtaining good agreement with 
experiment, provided, however, that the coupling constant was restricted to 
” = 10.36+0.02. 

Klein (1953a, 6) has shown that the term in [K,(ur)]}* should not appear in 
the fourth order central potential. This term contributes about 10% to the 
central potential. Calculations have been performed by Blatt and Kalos 
(1953) and by Jastrow (1953) in the low energy region neglecting this term, 
and allowing different coupling constants , \’ for the second and fourth 
order terms. They found no potential that simultaneously fitted all the low 
energy parameters. Nevertheless, the Lévy potential appears to be of some 
value, even if it is considered only as a phenomenological potential giving 
reasonable agreement with the experimental data up to 40 Mev. The calcula- 
tion of G? and r, has been repeated, using the potential of equations (1), (2), 
and (3). The reason for using the original potential, rather than Klein’s 
modification, was that the computations reported here were nearly completed 
by the time Klein’s work appeared. The calculations were performed on 
FERUT, the digital computer at the University of Toronto. 

The method is the same as that of Lévy. For various choices of the coupling 
constant, the radial equation for the 'S state of the NP system was integrated 
inwards, starting from the asymptotic form 


(4) Up = 1—7r/d gy. 


The integration proceeded until « = 0, which determined the core radius r,. 
In this way r, was determined as a function of X?. 

This procedure was repeated for the coupled *S and *D radial equations for 
the deuteron, starting from the asymptotic radial wave functions 


u = Aexp(—yr), 


w= pA exp(—1")| 1+ 3 oo $ | ; 
yr (yr) 

where u and w are the 4S and *D radial wave functions, respectively, and 
y? = Me/h?. The parameter p was varied until both radial wave functions 
vanished at the same radius r,. This again determined 7, as a function of ?*. 

The numerical integrations were performed by the predictor—corrector 
method of Milne (1949). The error in the solution was estimated by performing 
runs with different interval sizes and different starting points for the integra- 
tions. The effect of variation of the meson mass and of the scattering length 
Asn Was investigated in the 'S case. The effect of the experimental uncertainty 
in Gs, Was negligible. 

The two curves obtained from fitting the deuteron binding energy and the 
singlet scattering length are given in Fig. 1, in which 2? is plotted against r,. 
The singlet case was done for meson masses of m, = 274 m,, 276 m,, and 
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278 m,, the deuteron for m, = 276 m,. The experimental parameters used were 
Asn = —23.69X10-" cm., € = 2.226 Mev. 
The following points should be noted regarding Fig. 1: 


' (a) (b) 
S scattering length deuteron 


3 


COUPLING CONSTANT 
wo 


@ 





30 32 34 36 38 40 42 44 46 
CORE RADIUS X, (units — fi/m,c) 


Fic. 1. Coupling constant plotted against core radius for the Lévy potential, as determined 
from (a) the 1S NP scattering length, (6) the binding energy of the deuteron. 


(a) The two curves are very nearly parallel and do not intersect, at least 
for the range of core radius used here. 

(6) For a given A’, the core radius obtained from fitting the deuteron binding 
energy exceeds that obtained from the 'S scattering length by approximately 
0.045. 

(c) The estimated error of 7, in the calculations is approximately 0.0001 for 
the |S state and 0.004 for the deuteron. 

(d) The effect on r, of changing the meson mass by two electron masses in 
the |S case is 0.0014. The effect is expected to be very nearly of the same size 
and in the same direction for the deuteron, since the potential will be changed 
in similar ways in both cases. 

(e) The point of intersection quoted by Lévy lies between the curves obtained 
here. 


The separation of the curves of Fig. | is approximately five times the total 
error of the computations. Thus it must be concluded that one cannot fit 
simultaneously the singlet scattering length and the deuteron binding energy. 

An independent check of our calculations is contained in the work of 
Jastrow (1953), mentioned above. For \? = X” = 10, he obtained a singlet 
core radius of 0.34 and a triplet core radius of 0.38. Both of these are slightly 
less than those shown in Fig. | (as is to be expected, since he used a slightly 
less attractive potential), but the difference is the same as obtained here. 

The possibility of a fit to experiment by the Lévy potential cannot, however, 
be excluded by the present work. This question has been discussed in Paper I] 
of this series (Shapiro and Preston 1956). 
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MEASUREMENT BY PHOTOELECTRIC FRINGE SCANNING OF 
THE PRESSURE SHIFT OF Hg-198 EMISSION LINES! 


Kk. M. Batrp AND D. S. SMITH 


ABSTRACT 


Photoelectric fringe scanning apparatus has been developed at the National 
Research Laboratories for the measurement of very small wavelength changes 
(of the order of 5X10-5 A). This paper describes its application to a direct 
measurement of the shift of some of the visible lines of Hg-198 caused by the 
presence of the various rare gases used as a carrier. The values found for the 
shift in wave number are all of the order of 2X10~4 cm™ per mm. (Torr) pressure 
and to the red, except for He, for which it is somewhat less than 1 X 1074 cm™ per 
mm. and to the blue. 


INTRODUCTION 

A spectral line that is to be used as the primary standard of length must 
be very narrow and symmetrical to enable precise determination of its wave- 
length. A second requirement is that any variations of the wavelength depen- 
dent on conditions of excitation or observation must be small and well known. 
One source of lines satisfying the first requirement is the electrodeless Hg-198 
lamp developed by Meggers and Kessler (1950), who proposed that its bright 
green line be used as the primary standard. 

The wavelengths emitted by this lamp are shifted toward the red by the 
presence of argon, which is normally used as a carrier gas. The amount of the 
shift has been inferred by Barrell (1951) from separate determinations of 
wavelengths of various lamps at the National Physical Laboratory and from 
wavelength measurements made at other laboratories. He obtained the 
value of about 410-4 cm~! per millimeter pressure of argon for the wave 
number shift. However, the accuracy of this is not great because of errors in 
separate wavelength determinations as well as in the uncertainty of the argon 
pressure in the lamp at the time of observation. 

Margenau and Watson (1936) give a detailed treatment of the theory 
of broadening of absorption lines and Barrell (loc. cit.) has summarized in- 
formation relevant to shift and broadening of the emission lines of Hg-198. 
Under the conditions of use of Hg-198 lamps with which the present paper 
is concerned, broadening is mainly due to the carrier gas and is expected to 
vary linearly with its density. 

This paper describes a direct measurement of wavelength shifts for an 
electrodeless Hg-198 lamp whose carrier gas pressure could be changed and 
measured at will. The shifts were observed as changes in the fractional order 
of interference in a Fabry Perot etalon by means of apparatus for scanning 
the interference fringes photoelectrically. The latter was an adaptation of 
apparatus previously described by Baird (1954) in application to length 
measurement. Its sensitivity in a single setting was about +210-5 A. 


1Manuscript received January 7, 1957. , 

Contribution from the Division of Applied Physics, National Research Council, 
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Pressure shifts in the range 0 to 10 mm. pressure of Xe, Kr, Ne, and He as 
well as A were measured for several of the visible lines of the Hg-198 spectrum. 


FRINGE SCANNING APPARATUS 


The apparatus is shown schematically in Fig. 1. Light from the electrodeless 
lamp was collimated and passed through the etalon, which was mounted in a 





Fic. |. Schematic diagram of fringe scanning apparatus. 


pressure-temperature controlled chamber, C. The circular fringe system was 
focused and then centered on the slit of the spectrograph at S where a very 
small aperture allowed only light from the center of the rings (corresponding 
to 1/50 order) to pass into the spectrograph. At the back of the latter a 1P21 
photomultiplier tube received light from the spectral line under study. 

The pressure in the chamber C, and therefore the optical separation of the 
etalon plates, was controlled by two pistons, P; and P:. The piston P2 was 
operated by a synchronous motor so as to cause the pressure to oscillate at 
15 c.p.s. and at an amplitude such that the corresponding oscillation in the 
order of interference was approximately equal to the half-width of the bright 
fringes. The piston P; was operated by hand to adjust the average pressure 
and was equipped with an accurate scale and vernier so that its displacement 
could be read. 

The etalon, of invar, was of an open cage type so that pressure changes 
and oscillations were transmitted equally to the inside and outside surfaces 
of the quartz plates. The average temperature in the chamber was controlled 
to within a few thousandths of a degree and was not subject to rapid drift. 

The output of the photomultiplier was fed to a battery powered pre- 
amplifier Aj, with an L-—C filter to remove components much different from 
15 c.p.s. The signal was then fed through the Q multiplier A,, which removed 
components near 15 c.p.s., and thence to an oscilloscope swept at 30 c.p.s. 
in phase with the oscillating piston. This arrangement gave a sensitive 
indication of 15 c.p.s. modulation of the light falling on the photomultiplier, 
which was in phase with the piston oscillation. The signal was observed as a 
separation of the positive and negative parts of a sinusoidal signal on the 


oscilloscope. 
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The principle of wavelength shift measurement is shown in Fig. 2. The 
curve represents the variation of total flux entering the phototube with 
variation of pressure in the etalon chamber. Y; and Y, represent the limits in 
the scan amplitude due to the oscillating piston. In general when Y; and FY2 





p> 


Fic. 2. Variation in flux falling on photomultiplier with pressure in etalon chamber. 


were not symmetrical with respect to the peak intensity there was a 15 c.p.s. 
component to the light modulation due to the difference in intensity at Y; 
and Y2 and this was indicated on the oscilloscope. On the other hand when 
the average pressure in the chamber was adjusted by means of P; so that 
the intensities at the limits of the scan were equal, as at Y,' and Y,’, there 
was no 15 c.p.s. component; since the higher frequency modulation due to 
sweep over the peak of the curve was largely filtered out by the electronic 
circuit, a minimum signal was observed on the oscilloscope. 

In order to measure a pressure shift, the piston P; was set to give minimum 
signal (on a bright fringe) under given conditions of the lamp; the pressure 
in the lamp was changed and the piston displacement required to restore 
the minimum signal was found. The change in order of interference (and 
hence the pressure shift) was then determined from a knowledge of the piston 
displacement required to change from one bright fringe to the next. 

This apparatus made possible very convenient and precise determinations 
of wavelength changes due to lamp operating conditions. In a typical series 
of settings on the null point at 1-minute intervals over a period of 30 minutes 
the standard deviation of single settings from the mean was less than 0.001 
of an order with a 62.5 mm. etalon. This corresponds to a proportional wave- 
length change of about 4:10° and included scatter due to temperature varia- 
tion, backlash and reading errors of P;, lamp fluctuations, etc. 

The above sensitivity was obtained for the Hg green and blue lines; it was 
somewhat less for the yellow lines and much less at the red end of the spec- 
trum, e.g. at the Cd line 6438 where the photomultiplier sensitivity is very 
low. Sensitivity could probably be increased by the use of a red-sensitive 
cell, by refrigeration, and by the use of a detection system of greater sensitivity 
to phase. However, the first two changes do not have as much benefit as one 
might expect because the noise, which is largely due to the steady light flux, 
is not greatly reduced in relation to the signal. In any case the apparatus as 
described was found adequate for the present investigation. 
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An important feature of the method of fringe scanning by the use of the 
central spot of the Haidinger fringes is its great advantage over methods which 
involve measurements on the outer rings in the elimination of systematic 
errors due to vignetting, lens distortion, asymmetrical energy distribution 
in the rings, etc. On asymmetrical spectral lines the wavelength setting will 
of course depend to a certain extent on the scan amplitude Y, Y2 in Fig. 2, 
but this uncertainty is inherent in the spectral line and must be taken into 
account for any method of measurement. 

Scanning by 15 c.p.s. oscillations is felt to have considerable advantage 
over d-c. type scanning, as described by Jacquinot and Dufour (1948) and 
others, in applications such as the present one where drift would introduce 
important errors. In its present form it does not permit study of line shape 
but an adaptation for this purpose is under development. It is also pointed 
out that the principle of scanning with the advantages mentioned in the 
preceding paragraph is not limited to photoelectric technique, but can be 
adapted to photography with the added advantage of recording many lines 
at once. 

PROCEDURE 

In order to make the pressure shift measurements a Meggers’ type Hg-198 
‘lamp containing about 3 mg. of Hg-198 of 98% isotopic purity was directly 
connected to a vacuum line and to sources of He, Ne, A, Kr, and Xe which 
could be metered out in small quantities. The pressure of gas in the lamp 
was read by means of a calibrated oil manometer. The lamp was separated 
from the rest of the system by a liquid air trap except in the case of Kr and 
Xe for which a CO, trap was used. 

Excitation was by means of a special 270 Mc./s. oscillator previously de- 
scribed by Smith (1957). The lamp was water cooled to operate at a tem- 
perature of 10°C. No systematic measurement was made of the effect of 
excitation or temperature on the wavelength. However, deliberate variation 
of these factors, which was much greater than any likely to occur during the 
pressure shift measurements, was found to produce no observable effect. 

The procedure to obtain a single point on the pressure shift — pressure 
curve was as follows: first the lamp was operated at well below 0.1 mm. gas 
pressure and three null settings of P; taken at 1-minute intervals; the lamp 
pressure was then raised to the value for which the shift was being measured 
and three more settings made at 1-minute intervals; finally, the lamp pressure 
was reduced to well below 0.1 mm. again and three more readings taken. 
This procedure was followed for all pressure values for which measurements 
were made according to a previously determined random schedule. 

For argon, measurements were made at all the points indicated in the 
graph in Fig. 3 on the lines 4047, 4358, 5461, 5770, and 5790 using both a 
62.5 mm. etalon and a 25 mm. etalon. Measurements were also made of the 
difference between a cold cathode Hg-198 lamp supplied by N.P.L. and the 
Meggers’ type lamp at near zero pressure. 

For the other gases points were obtained only for 0 mm., 5 mm., and 10 mm. 
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Fic. 3. Wave number shift of Hg-198 line 5461 with argon carrier gas pressure. 


lamp gas pressure, and only for one or two lines in each case because poor 
excitation or interfering lines made observation impractical. 
RESULTS 

A typical pressure shift curve is shown in Fig. 3, which shows the shift 
of the 5461 Hg-198 line with argon pressure. The variation is linear as expected 
and has the value 2.0X10-* cm~'!/mm. The probable error of this value cal- 
culated from the scatter of single points is +0.1 10-4 cm™ per mm., although 
the difference in the several independent determinations over a period of a 
few months amounted to about 0.21074 cm. 

The wave number difference between the Geissler type lamp supplied by 
N.P.L. and the low pressure Meggers’ type lamp was found to be 251074 
cm! for the 5461 line. As the N.P.L. lamp was said to be filled to 10 mm. 
pressure of A, the shift per mm. appears somewhat higher than that in- 
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dicated by the measurements on the electrodeless lamps. This may be due to 
the different method of excitation. 

The values for all the lines measured are given in Table I. It will be seen 
that except for He the shifts are all to the red and of the same order of magni- 








TABLE I 
d (A) Shift in cm™!X10~4 for 10 mm. (Torr) carrier gas pressure 
25.0 mm. 62.5 mm. 

etalon etalon He Ne A Kr Xe 
5791 26+1 

5791 17+1 
5770 30+8 

5770 19+1 
5461 25+1 

5461 +7.6 —I1 20+1 —21 
4358 29+2 

4358 +6.6 —15 23+1 — 20 
4047 21+1 

4047 21+1 


tude. The values for A are less than that obtained by Barrell, but together 
with the value for Ne are in agreement with those found by Kessler (1956) 
at the National Bureau of Standards. 

The interesting blue shift for He carrier gas suggests the possibility of 
making a lamp containing a mixture of He and A in order to have no net 
wavelength shift. This has been found to work, but adjustment of the mixture 
of gases must be done carefully and an envelope impervious to helium would 
be required if the lamp were to have a long life. Whether the two gases would 
clean up at the same rate has not been determined. 

The slight difference observed in the apparent shift with the two different 
etalon lengths is thought to be real and to result from the asymmetrical nature 
of the broadening. 

The magnitude of the shift is such that it must be taken into account when 
the usefulness of the Meggers’ lamp as a source for the primary standard is 
being considered. Lamps filled to 3 to 6 mm. pressure as were the earlier 
ones would have wavelength errors of the order of nearly one part in 10? 
which would likely change during the life of the lamp. This sort of error 
would be of importance in precise metrology. 

On the other hand lamps filled to only one or two tenths of a millimeter 
pressure as in Meggers’ later lamps could not be in error by more than one 
or two parts in 10° because of pressure shift, an amount which would likely 
be acceptable at present in a primary standard. 

It is of interest to note the corresponding shifts in two other sources proposed 
as the primary standard. The red line of Cd-114 has been found by Batar- 
chukova and Dubrovskii (1956) to suffer a red shift of 110-4 cm. per mm. 
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argon, which is used as the carrier gas. The Kr isotope lamps suffer a shift 
to the red due to self-broadening that follows the two-thirds law. It amounts 
to less than 10-° cm—' for the 5650 line at a Kr pressure equal to its vapor 
tension at the triple point temperature of Ne (which is the temperature of 
normal operation). Thus the pressure shift of the Hg-198 green line under 
present operating conditions is somewhat greater than for these other 
proposed standards. 
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ELECTRICAL RESISTIVITY OF GRAPHITE 
IRRADIATED WITH NEUTRONS! 


R. W. ATTREE AND O. DAHLINGER 


ABSTRACT 


Measurements have been made of the changes in the electrical resistance 
of an artificial graphite which are induced by neutron irradiation and subsequent 
annealing. The resistivity increases with irradiation; the increase in resistance 
is partially removed by annealing at 128° C. These changes have been inter- 
preted using Wallace’s theory of the electron structure of graphite. The number 
of atoms displaced by the irradiation has been calculated using the method of 
Seitz, and is shown to be consistent with the experimental results. 


INTRODUCTION 

The physical properties of graphite are markedly changed by irradiation 
of the material with fast particles. The properties which have been investi- 
gated include the thermal and electrical resistivity, Hall coefficient, Young’s 
modulus, thermoelectric power, density, and lattice spacing.? Of these, the 
change in electrical resistivity is the most convenient to measure, although 
the interpretation of the results is difficult. This difficulty arises from the 
nature of the conduction process. The electronic conductivity of a solid 
‘ depends both on the density of electrons which are available for conducting a 
current, i.e. the density which can be thermally excited above the Fermi 
level, and upon the effective path length which they traverse before being 
scattered by grain boundaries, thermal waves, or lattice imperfections. The 
displaced atoms and vacancies which result from irradiation affect the con- 
ductivity in two ways—they act as scattering centers, thus decreasing the 
mean free path of the electrons, and also trap electrons, thus changing the 
density of electrons available for conduction. 

The changes in resistivity may be partially removed by annealing at tem- 
peratures above the irradiation temperature. It may be expected that both 
the density of scattering centers and the density of trapped electrons will be 
altered by this treatment. 

We have measured the changes of resistance in graphite induced by neutron 
irradiation, and by subsequent annealing, and, using Wallace’s (1947) theory 
of the electronic structure of graphite, we have shown that a satisfactory inter- 
pretation of the results can be given in terms of the density of scattering 
centers and of conduction electrons. 


EXPERIMENTAL METHOD 


The specimens were cut from a single block of artificial (petroleum coke) 
graphite, perpendicular to the direction of extrusion. They were of circular 
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cross-section with dimensions shown in Fig. 1. Current and potential leads 
were attached using phosphor-bronze clips with jaws modified to clamp on 
the circumference of the specimens. 


prone eal 


Fic. 1. Dimensions of graphite samples. 


The resistance was determined as follows. Three or more measurements 
of the voltage and current (measured as voltage across a standard manganin 
resistor immersed in oil) were taken for several rotational positions of the 
specimen in the clips, and these were averaged. If the result for a position 
differed from the average by more than 1 part in 500, indicating poor contact 
between the clips and the specimens, readings were taken for more positions 
until consistency was obtained. The values for the resistances for each specimen 
were then averaged, those which indicated poor electrical contact being 
rejected. 

The resistance of two standards was determined in the same manner, 
whenever specimens were measured, and the specimen resistivities were 
corrected for any apparent change in the standards. These corrections were 
of the order of one per cent over a period of several months. 

Unirradiated specimens were annealed at 128° C. before use. All specimens 
were placed in the dry box surrounding the measuring equipment for several 
hours before measurement. All measurements were made in dry air at 29° C. 

Figure 2 illustrates the sequence of measurements in a typical series. The 
resistances of the three specimens (denoted by A, B, and C in Fig. 2) were 
measured, and the specimens were then irradiated in the same facility in the 
NRX reactor. After the irradiation, the resistances were measured. The 
ratios of the resistances after and before irradiation are denoted by the first 
point on the upper curve. (These ratios may be taken as the ratio of volume 
resistivities p/po, since dimensional changes are negligibly small.) Sample A 
was then annealed at 128° C.; the change in the ratio of resistances is in- 
dicated by the vertical broken line. The three samples were irradiated as 
before, and remeasured. Samples A and B were then annealed. It will be 
observed that at this stage the values of p/p are the same for A and B, within 
the error of the measurements, even though the sequence of annealing and 
irradiation differs. The samples were then irradiated, measured, annealed, 
and remeasured, as indicated in Fig. 2. In all, nine specimens, designated 
A to J, were measured, using three different irradiation facilities in the NRX 
reactor. The results are given in Table I. 
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IRRADIATION DOSE (NEUTRONS /cm*) 

Fic. 2. Resistivity changes in irradiated graphite. The upper curve shows the changes 
induced by irradiation with no annealing, while the lower curves show changes produced on 
reirradiating specimens which have been annealed at 128° C. Behavior on annealing is in- 
dicated by the vertical broken lines. When more than one sample was measured, the spread 
of the results is indicated. 


TABLE I 


Irradiation position ¢, neutrons/cm.? 


S-15-5 : al A Cc 
2.66 10'8 17 168 


6.40108 341 349 
.167) (a) 


11.07 10'8 : .413 (a) 004 
. 254) (db) .275) (a) 


E F 
3.38 X 10'8 205 203 


8.40 X 108 ‘ .418 408 
201) (a) 


13.65 X 10" : .459 (a) Oo 
312) (0) 326) (a) 


H J 
2.5010" 1.086 1.097 


5.07 XK 10'8 1.150 1.164 
(1.065) (a) 

9.00 X 108 1.194 (a) 1.257 
(1.108) (0) (1.130) (a) 


*Values enclosed in brackets are for specimens annealed after irradiation to the total dose 
indicated. Values for specimens annealed once following an initial irradiation are designated 
by (a); for those annealed twice, by (6); and for those annealed three times, by (c). The 
irradiation positions S-9-4, S-9-5, and S-15-5 are located in the graphite reflector outside the 
calandria of the NRX reactor. 
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Measurements were made during the course of annealing, which, in some 
cases, was continued for 10 days. It was found that after long annealing, 
further changes in the resistance were proportional to the reciprocal of the 
annealing time. This fact was used to obtain the values for resistance after 
infinite annealing times, shown in Fig. 3. All samples were annealed for at 
least 5 davs. 





2 4 6 8 10 12 14 x 10 
IRRADIATION (NEUTRONS/ cm?) 


Fic. 3. Ratio of the resistivity after irradiation and annealing at 128° C., pa, to the resistivity 
before irradiation, po, as a function of the integrated flux. The lines drawn are the ‘‘best’’ fit 
as determined by the method of least squares. The lines for irradiation in S-9-5 and S-15-5 
facilities coincide. 


Thermal neutron fluxes were measured with cobalt monitors. The tempera- 
ture of the graphite during irradiation was about 70-80° C. 


DISCUSSION 


One point is immediately evident from the data, that within the accuracy 
of the measurements a sample which has been irradiated, annealed, re- 
irradiated, and reannealed has the same resistivity as a sample which has 
first received the same total amount of irradiation and which has then been 
annealed. It is possible to imagine that the aggregates of imperfections which 
remain after irradiation and annealing will either serve as nuclei for the 
further aggregation of the imperfections produced by a second irradiation, 
or will be disrupted by the second irradiation. Each of these possibilities 
will lead to a state different from that arising from a single irradiation and 
annealing. Because no such difference is found, we conclude that, for the 
irradiations used, the imperfections introduced by the irradiation do not 
interact appreciably with the aggregates remaining after the annealing, nor 
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are these aggregates affected by further slight irradiations, within the accuracy 
of the measurements. This result is to be expected for the short irradiations 
used, but would not be expected for longer irradiations. 

Using Wallace’s (1947) theoretical relations for the energy dependence 
of the density of the electronic states in graphite, Hennig (1951) has derived 
an equation for the volume resistivity p: 

Ser 


=a [2kT In{1+exp(A/kT) } —A], 


(1) : 
p 
where 1/7 is the probability per unit time of scattering an electron, A is the 
difference between the Fermi energy and the energy at a corner of the first 
Brillouin zone, ¢c is the c-axis lattice spacing, and the other symbols have 
their usual meaning. Since the value of A depends primarily on the number 
of electrons per atom, and the fraction of these which are trapped, and 1/r 
depends on the density of scattering centers, it may be seen that equation 
(1) is a quantitative statement of the relation between the resistivity, the 
density of scattering sites, and the fraction of electrons trapped. The relation- 
ship between the number of electrons trapped per atom of the solid, m, and 
A has been given by Hennig (1951) as 


2.2m _ tA 4a 
RT 2NkT. 6 
if RT and A are expressed in electron volts. 

We shall now determine values of A, and hence m and 7, using these equations 
and the data observed. Let po be the initial resistivity of a specimen; p;(¢) 
the resistivity after irradiation in an integrated flux ¢; 1/7» the probability 
per unit time of scattering an electron wave, due to thermal effects and to 
lattice defects of all kinds which are present in the unirradiated graphite; 
1/7:(@) that due to imperfections introduced on irradiation; and 1/74(@) 
that due to imperfections, other than those in the original lattice, which 
remain after annealing. 

If we assume that the imperfections responsible for the scattering are 
independent, that is, that the displaced atoms and vacancies formed on 
irradiation or remaining after annealing are sufficiently far from other scatter- 
ing centers of all sorts that they do not overlap, and that the concentration 
of imperfections resulting from irradiation is proportional to the integrated 
flux ¢ received, we may define constants S; and S, such that 


(3a) 1/ri(o) = Sid, 1 /talb) = Sad 


and, further, 
(3d) 


We shall also assume that the number of electrons trapped per atom of carbon, 
m, is proportional to the integrated flux, and write 


(4) m = ud = Nod, 
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where is the number of electrons trapped per displaced atom and a is the 
cross section per atom for the displacement of atoms by neutrons. We shall 
neglect the change in c-axis spacing which appears in equation (1), as this 
can be shown to introduce a negligible correction to the results. 
Substitution of these terms in equation (1) leads to the expression 


pi(d) 2In2 
Po 


" = (14S: 106) Sint pexp(a/eF))CA/ET 


and an analogous expression for pq(¢) /po. 

Experimental values for p,(¢)/po are plotted in Fig. 3 as a function of ¢. 
It is apparent that the ratio of the resistivity of a specimen annealed at 128° 
after irradiation, to that of the initial unirradiated resistance, is proportional 
to the integrated flux received, and hence cannot depend on A. Regardless 
of the form of the dependence of the resistivity on A, this can only mean that 
A for such specimens is the same as that for unirradiated specimens. We 
may therefore conclude independently of Wallace’s theory that the imper- 
fections remaining after annealing do not trap appreciable numbers of elec- 
trons. Values for S,7) obtained from the slopes of the lines in Fig. 3 are given 
in Table IT. 











TABLE II 
VALUES OF THE PARAMETERS GIVING A “‘BEST’’ FIT IN EQUATIONS (3), (4), AND (5) 
Irradiation facility S-9-5 S-15-5 S-9-4 
Sito (cm.?/neutrons) 6.2X10- 6.1xX10™ 3.310" 
Saty (cm.2/neutrons) 2.2X10-" 2.2x«K10-" 1.210-% 
uw ~©6(cm.? electrons/atom neutrons) 3.5X10°% 2.8xX10-% 1.610~-% 
Sito/u (atoms/electrons) 1.7X10? 2.210? 2.010? 
*/ Soa 2.8 2.8 2:2 
n  (electrons/displaced atom) 1.4 1.1 Bul 





A sensitive test of the usefulness of Wallace’s theory is provided by com- 
paring the ratios S;/S, and S4ro/u obtained in the three irradiation facilities. 
These ratios are independent of the intensity of the flux, although they will 
depend on the form of the flux spectrum. Since the form of the spectrum is 
expected to be very similar in the three positions, constancy of these ratios 
will be a more stringent test than showing that the data for any one facility 
can be fitted by an equation of the form (5) by suitable choice of parameters. 

Values of S;7>) and yw which give the ‘‘best’’ fit to the data may be obtained 
for this purpose as follows: a value for u is assumed; from this and equations 
(2) and (4) the values of A(¢) are determined for the appropriate ¢. These 
calculated values and the corresponding experimental values of pi()/po 
are then used to determine, by means of a least-squares fit, the value of Sy. 
The value of » which gives the smallest root-mean-square deviation of the 
experimental points from a straight line, as required by equation (5), is 
selected as the “‘best’’ value. The sensitivity of the root-mean-square deviation 
to change in yp is shown in Fig. 4 for a typical case. The values of the param- 
eters determined in this way are given in Table II. 
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Fic. 4. Root-mean-square deviation of the data from the ‘‘best’’ straight line whose slope 
is Siro for data for S-9-5. w is expressed in cm.? electrons/atom neutrons, while d is, of 
course, dimensionless. 


It is very desirable to determine a value for the number of electrons trapped 
per displaced atom, n. This we can do from values obtained for pp = no, 
by calculating o, the cross section per atom for the displacement of atoms, 
using a method given by Seitz (1949). This method requires a knowledge 
of the flux spectrum, which we shall assume is given by 
goo dE 
= ae 


n(E)dE = EB? 


where @¢ is the thermal neutron flux as measured by cobalt monitors, FE is 
the energy of the neutron, and c is the ratio of the density of thermal neutrons 
as defined above to the density of neutrons in an energy interval In E. We then 
find, following Seitz’s (1949) procedure, that the fraction of the atoms dis- 
placed per unit flux, o, is 1.0X10-*'/c displaced atoms per atom/neutrons 
per square cm. According to activation experiments by Butler (1956) on gold 
and cobalt with and without cadmium shielding, values for c are 40 in S-9-5 
and S-15-5 positions, and ~66 in S-9-4. Corresponding values for n, the 
number of electrons trapped by each displaced atom, are given in Table IT. 

Considering the approximations involved, the agreement between the 
three values for S;79/u, S;/Sq, and n is surprisingly good. No reliance should 
be put on the numerical values for ”, as they may all be in error by a common 
factor of at least two. This error arises both from the approximations used 
in the calculations, and from the assumption that all the atoms displaced 
during irradiation remain as displaced atoms at room temperature. Austerman 
and Hove (1955) have shown that the resistivity changes in graphite irradiated 
at liquid helium temperatures decrease on annealing to room temperature, 
but it is not vet established that this is due to annealing or clustering of 
displaced atoms. We feel, however, that the values of so obtained are 
physically reasonable, and conclude from this that the method which Seitz 
gives for calculating the number of displaced atoms gives at worst a lower 
limit to the number produced. This is true because, although both the lattice 
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vacancy and the displaced atom may trap electrons, it is unlikely that more 
than three or four electrons would be trapped for each atom displaced. Since 
it is most probable that the atoms displaced by irradiation at temperatures 
above 0° C. are clustered, we cannot set a lower limit on the number of 
electrons trapped by each displaced atom. We can, however, conclude that 
the method of Seitz gives values for the number of displaced atoms which 
are not too low, for otherwise the values for » would be unreasonably large. 

In conclusion, then, we may state that, independent of the theory used, 
the imperfections left in graphite after short neutron irradiations, followed 
by an anneal at 128° C., do not trap electrons; that the equation (1) derived 
by Hennig, based on the Wallace tight binding theory for graphite, can be used 
to explain the observed facts; and that the method given by Seitz for estimating 
the average number of atoms displaced by an incident neutron gives in this 
case values which are not unreasonable, and which are not too small. 
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PHOTONEUTRON EMISSION FROM Th?, U*, Us, AND Pu*?! 


L. Katz, K. G. McNEILL,? M. LEBLANC,’ AND F. Brown? 


ABSTRACT 
_ Measurements have been made of the numbers of neutrons emitted by the 
fissile materials Th??, U2, U238, and Pu®® under irradiation by betatron X-rays 
of various energies up to 23 Mev. Yield curves and corresponding cross section 
curves are given for these nuclides. At low X-ray energies it is found that the 
yield curves differ from those found with non-fissile heavy elements, corresponding 
to the fact that even at zero X-ray energy the fissile materials suffer (spontaneous) 
fission. 

INTRODUCTION 

The measurement of neutron emission under betatron X-ray irradiation 
by the fissionable nuclides Th?, U?5, U*8, and Pu is of interest both as 
a continuation of the program of photoneutron investigation carried out in 
this laboratory and because these data may yield interesting results when 
analyzed in the light of investigations of photoneutron reactions in non- 
fissile heavy elements and photofission reactions in these nuclides. 

In the study of neutron emission from photoexcited fissile materials there 
can be no clear-cut distinction between photoneutron and photofission re- 
actions, as some excited nuclei may decay by the emission of one or more 
neutrons followed by the fissioning of the still-excited daughter nucleus. 
In fact the total photon capture cross section is 


Srotal = Cryin Oya Oyen tae Hoyt Oyng t+ oy2ngt tea 


Other reactions, involving the emission of charged particles, may be neglected 
owing to their inhibition by Coulomb forces. Following Gindler, Huizenga, 
and Schmitt (1956), it is convenient to use the terms: 

oyn = Cran Lea tee + v0y¢+ (v1) oyngt+ see 
and 


OyFr = Cy Cyap A Oyong eee 


to denote the effective cross sections for neutron emission (¢,y) and photo- 
fission (oy). oy, and oy, will be used, as above, to denote single neutron emission 
or simple fission. In the expression for o,y, v is the average number of neutrons 
emitted when a nucleus undergoes simple fission (roughly »v = 2.5), and 
should not be confused with the average number of neutrons that are obtained 
for every one fission that takes place in a bulk of irradiated material, that is: 
0yv/oyr. v Will probably itself be a function of X-ray energy. 

Measurements of the photofission (¢,) cross section of U*** have been 
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made by Katz et al. (1955) using chemical separation of the fission products, 
and similar work has been carried out at the University of Illinois (Gindler, 
Huizenga, and Schmitt 1956) in U8 and Th2”. The latter authors also meas- 
ured radiochemically the cross sections for the reactions U**5(y, 2)U"37 and Th?” 
(y, 2) Th?*!, Photofission experiments on a number of fissile materials, including 
U33, U3, Th*®, and Pu*’, are at present being carried out in this laboratory. 

Measurements of neutron emission yields from U have been made by 
Nathans and Halpern (1954) and from U and Th by Price and Kerst (1950). 
In Russia, Lazareva, Gavrilov, Valuev, Zarsepina, and Stavinsky (1955) 
have measured the yield curves and cross sections for neutron emission 
(y, V) for Th and U in the energy range 6-28 Mev. 

In the present experiments the (y, NW) yield curves and corresponding 
cross sections have been measured from 23 Mev. down to approximately 
5.5 Mev. in the cases of U** and Th?”, and down to approximately 6.5 Mev. 
in the cases of U*8 and Pu®?, 

EXPERIMENTAL PROCEDURE 

The paraffin-embedded BF; long-counter apparatus used in this experiment 
to detect the emitted neutrons has been previously described (Montalbetti, 
Katz, and Goldemberg 1953). In order to reduce neutron background, mainly 
produced by (y, 7) reactions in the lead of the betatron beam collimator, 
the higher energy work was carried out in a room separated from the betatron 
by a 7 foot concrete wall, through which a hole was drilled for the X-ray 
beam. At lower energies the betatron room background is lower, and in order 
to increase the measured yield in the range of 5-8 Mev., work was carried out 
in this room. As an additional precaution, the lead collimator was replaced by 
one of iron. 

In the original neutron detection apparatus only one BF; counter embedded 
in paraffin was used. This was so positioned that its response was insensitive 
to the initial energy of the neutrons. The addition of three other BF; counters 
has been found to increase the detection efficiency of the apparatus fivefold, 
and not to introduce any appreciable errors due to the fact that these counters 
were not specifically located in neutron energy insensitive positions (LeBlanc 
1954). Therefore the setup with four counters was generally used, but as an 
additional precaution only the energy insensitive counter was used during the 
normalization runs in which the neutron emissions from the fissile materials 
were compared with the neutron emission from copper. 

The samples were highly purified samples of Th, U, and Pu shaped in the 
form of disks. Each disk was embedded in polystyrene, which was placed in a 
thin cylindrical holder of Al. An aluminum cylinder of the same size as the 
others, but containing only polystyrene, was used to determine the reaction 
background from the holders when under irradiation. 

Table I lists the compositions and areal densities of the fissile samples used. 
Although all the compounds used in this work contained oxygen, the contri- 
bution of this element to the number of neutrons detected is negligible, 
amounting at the most to less than 0.1% of the total, owing to the compara- 
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TABLE I 
Nuclide Compound —_—— 
A a ThO, 1.26 
U 233 U3;08 0.684 
[238* U;0s 0.760 
Pu? PuO: 0.779 


*This sample depleted in U2*. 
{This sample contained a small amount of Pu*®. 


tively high threshold energy and low cross section of the O'*(y, 7)O' reaction. 
Calculation of the correction factor for absorption of the bremsstrahlung 
within the sample has been described by Montalbetti et a/. (1953). In the 
present experiment, the necessary correction for X-ray absorption in no case 
exceeded 3%. 
Activity leading to neutron emission may be induced in a sample by 
bombardments other than by bremsstrahlung. These are: 


(a) neutrons emitted by the sample itself causing fission in the sample, 





(6) neutrons in the betatron beam causing fission, 
(c) neutrons emitted by the sample holders causing fission. 


All these possibilities were considered, and it was found that none of these 
sources could cause 1% of the measured number of neutrons. 

Not all emitted neutrons are detected by the BF; counters and the absolute 
efficiency of the apparatus was unknown. Use was therefore made of the fact 
that the number of neutrons emitted by natural Cu when irradiated with a 
22 Mev. maximum energy X-ray beam is known by activation methods to be 
2.71X10° n./100 r./mole.* A sample of Cu of areal density 0.911 g./cm.? 
was irradiated at 22 Mev., and the measured yield compared with the measured 
yields from the fissile samples under identical conditions. The absolute 22 Mev. 
yield for these samples was thus obtained, and yield versus energy curves for 
the samples were normalized at this energy. Bremsstrahlung absorption in 
the Cu necessitated a correction of 1% to the measured copper yield. 

The X-ray beam intensity was normally measured by an ionization chamber 
whose output was integrated. The X-ray energy response of this chamber was 
checked against a Victoreen chamber placed in an 8 in. lucite block. No 
absolute beam measurements had to be made, as all results were, as stated 
above, normalized against the 22 Mev. Cu yield. 


RESULTS AND DISCUSSION 
The yield curves of the four fissile materials studied, in neutron yield per 
100 roentgens per mole, with respect to betatron operating energy are shown at 
the tops of Figs. 1, 2, 3, and 4. The low energy ends of the yield curves for 
*The yield from natural copper used in the paper by Montalbetti, Katz, and Goldemberg 
(1953) was 3.01108 neutrons per 100 roentgens per mole. Measurements by Johns in our 


laboratory have shown our Victoreen monitor to read 10% too low at energies above a few 
Mev. The yield has thus been correspondingly reduced. 
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Fic. 1. Yield and cross section curves for the total neutron emission, (y, N), from Th? 
under betatron X-ray irradiation. 

Fic. 2. (y, N) yield and cross section curves for U*, 

Fic. 3. (y, N) yield and cross section curves for U8. 

Fic. 4. (y, N) yield and cross section curves for Pu, 
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Fic. 5. The low energy section of the U28, Th®?, and Bi? neutron yield curves plotted on 
a logarithmic scale. Compared with the previous figures, the energy scale is expanded. 


U3 and Th? are shown on an expanded scale in Fig. 5. These curves are 
believed to be accurate within the statistical accuracy of the individual 
points. This accuracy varies from about 2% at 22 Mev. to 3% at 12 Mev. and 
to about 30% at 6 Mev. 

To find the cross section curves corresponding to the measured yields the 
“inverted matrix’? method of analysis of Leiss and Penfold (to be published) 
was used. This is a development of the ‘‘photon-difference’’ method of analysis. 
The results of these calculations are shown at the bottom of Figs. 1, 2, 3, and 4. 

To some extent the smoothing of the yield curves, necessary before the 
cross sections can be calculated, is arbitrary and subjective. As a result there 
may be some distortion of the actual shapes of the cross sections. It is regarded 
that possible errors of 25% should be attributed to any one point of the cross 
section curve. However, such parameters as the energy of the maximum 
cross section, and the integrated cross section, fo(E)dE, are not very sensitive 
to the smoothing process or to the actual shape of the bremsstrahlung spec- 
trum chosen. A greater degree of confidence may thus be placed in these 
parameters. 

In Table II are given the values found in the present experiment of the 
reaction parameters: maximum cross sections, ¢max, the energies corresponding 
to these maxima, E(omax), the half widths of the giant resonances, I, the 
integrated cross sections, [°odE, and the 22 Mev. yields. For comparison, in 
Table III are listed corresponding values found by other authors, where 


these are available. 
From the tables it will be apparent that there is reasonable agreement 
between the results of the present work and those of the other authors quoted. 
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TABLE II 
22 Mev. yield, omax, Eline): r, SordE, 
Nuclide —_n./100 r./mole barns Mev. Mev. Mev-barns 
Th? 5.00 109 0.99 14.2 6.0 7.15 
U238 7.54X109 1.29 15.2 6.4 9.74 
U3 8.88109 1.67 14.0 6.0 11.2 
Pu®9 9.51X10° 1.58 13.6 6.3 11.6 
Cu (63 and 65) 2.71108 (normalization figure) 
TABLE III 
2 Mev. vield, 
n./100 r./mole Gian E(Gwsx); F, \ondE, 
Nuclide x< 10° barns Mev. Mev. Mev-barns 
Th2s2 a: \lay* 0.80+0.10 (c) 14.5 (c) 5.6 (c) 6.6+0.6 (c) 
[238 5:3" (a) 0.98 (b) 13.8 (b) 6.6 (6) 7.15 (b) 
5.75 (b) 1.18+0.15 (c) 14.9 (c) 6.8 (c) 12.9+1.0 (c) 


*(a) Price and Kerst 1950. (b) Nathans and Halpern 1954. (c) Lazareva et al. 1955. 


From Fig. 5 it will be seen that the initial portions of the uranium and 
thorium yield curves may be divided into three parts. In the high energy 
and low energy parts, the logarithm of the yield is roughly proportional to the 
photon energy, whilst the third and middle region is one of transition between 
the other two. This transition occurs around 6-7 Mev. The upper region from 
7 Mev. upwards lies above the (y, ”) threshold (6.0 Mev. for U** (Huizenga 
et al. 1951) and 6.4 Mev. for Th? (Magnusson et al. 1951) and also above 
what is often termed the fission threshold.* This region, indeed the whole curve 
above 7 Mev., shows a neutron yield — photon energy dependence which is 
typical of the photon absorption process. We have shown in Fig. 5 the initial 
portion of the bismuth photoneutron excitation curve (Chidley and Katz 
1956). This curve is typical for non-fissile heavy elements; for energies more 
than 1 Mev. above the (vy, ) threshold it corresponds to the total photon 
absorption curve since neutron emission is the dominant mode of de-excitation. 

The U* and Th?” curves above 7 Mev. and the Bi?®’ curve above the 
Bi(y, 2) threshold are clearly similar. In this region then we may suppose that 
the probabilities of fission and neutron emission in those U and Th nuclei 
which have absorbed photons do not vary rapidly with photon energy. This 
agrees with the fact that the ratio of (y, f) to (y, m) in U** is roughly constant 
over a wide photon energy range (8-22 Mev.) (Duffield and Huizenga 1953). 
The slope of the curves above 7 Mev. is thus governed chiefly by the increase 
in the photon absorption cross section with photon energy. It must be 
emphasized however that one should expect the neutron yield per absorbed 
photon from U and Th to decrease somewhat above the (y, 2) threshold. 
This is because below the (y, 2) threshold for each photon absorbed we have 
2.5 neutrons (fission), while above this threshold some of the absorbed photons 
result in only one neutron (y, 2) reaction. A rough estimate of the change in 








*This is a somewhat vague term. In the present context it may be taken to mean the energy 
above which fission accounts for a significant fraction of all de-excitation modes. 
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slope in the U and Th curves of Fig. 5 as the (vy, m) threshold is crossed shows 
that it cannot all be accounted for by the presence of (y, ~) competition on 
the higher energy side. This may indicate a considerable barrier effect against 
fission. 

As the photon energy falls below the (y, 2) threshold the bismuth neutron 
yield necessarily drops extremely rapidly to zero yield. The uranium and thor- 
ium curves pass through a transition region which is longer than the sharp 
‘“‘knee”’ in the Bi curve and reach a steep but not vertical region extending 
from 6 or 6.5 Mev. downwards in which the neutron yield is due entirely to 
fission. In this energy range we must suppose that the probability of fission 
following photon absorption is strongly energy dependent; probably fission 
barrier effects have now become important. The steep slope of this curve is 
therefore due partly to the decrease in photon absorption with decreasing 
energy and partly to the decrease in the probability of de-excitation by fission 
as the excitation energy is lowered. 

Useful data on Pu*® could not be obtained below 6.5 Mev. owing to the 
high neutron background caused by the spontaneous fission of Pu*?®. 
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THE GRAVITATIONAL FIELDS OF ELECTRIC AND MAGNETIC 
DIPOLES! 


GERALD E. TAUBER 


ABSTRACT 
Einstein’s field equations are solved for cylindrical symmetry in the presence 
of electric and magnetic singularities (dipoles). 


1. INTRODUCTION 


Amongst the few known exact solutions of Einstein's field equations are the 
ones representing the gravitational field of a charged mass or several charged 
masses distributed with rotational symmetry (Weyl 1917; Curzon 1925; 
Papapetrou 1947). It would be of some interest to find solutions corresponding 
to more general situations, say the fields due to rotating magnetic dipoles 
which occur in some astrophysical problems. Unfortunately the non-linear 
character of the field equations makes it extremely difficult to find exact 
solutions of the resulting system of equations. 

Recently Bonnor (1954) has shown that some known electrostatic solutions 
of the field equations can be adapted to the magnetostatic case by the introduc- 
tion of the conjugate field tensors F*‘’. However, if both electric and magnetic 
sources are present his solutions are not applicable. The property of cylindrical 
symmetry on the other hand permits us to include both types of sources in 
the equations. Moreover, the resulting equations for the electromagnetic 
potentials are of a very simple structure and can be solved by an extension 
of the work of Curzon (loc. cit.). It is the purpose of this note to discuss this 
extension of Curzon’s work and to obtain the solution of the field equations 
for electric and magnetic dipoles. 

2. THE FIELD EQUATIONS 

In canonical co-ordinates the line element for a field with axial symmetry 

is given by 


(1) ds? = —e*[(dx')? + (dx?)?] — e7°(x?)?(dx?)? + e°(dx*)?, 


where x! is the axis of symmetry, x? is a radial co-ordinate, x* the angular co- 
ordinate, and x‘ time-like. Both \ and p are functions of x! and x’ only. Ina 
region free of matter Maxwell’s equations are (cf. Tolman 1934) 


it i. ie. d(\/—¢ F”’) 
2 4 2 SS ee a se 
(2) Fs de” a” ax” . 


where A, are the components of the four-potential and F,, the electromagnetic 
field tensor. In accordance with the symmetry requirement of (1) A, are 
functions of x' and x? only. All non-zero components of the electromagnetic 
field tensor can be obtained from two components of the potential: 


'1Manuscript received November 30, 1956. 
Contribution from the Department of Physics, Western Reserve University, Cleveland, 
Ohio. 
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(3) y = A,and ¢ = 


With these conditions the surviving Maxwell’s equations (2) then become 


, a a(e = 
(2a) ax rts - ote 6:) = 0, 


(26) Siete A+ a(e" x po) = 0. 

The field equations in the presence of magnetic fields are given by (Tolman, 
loc. cit.) 

(4) Ri—4 gh R= —«T?, Te = — FF a tighk Fas, 


where « denotes the gravitational constant and R% is the Ricci tensor. For 
cylindrical symmetry they have already been given by Curzon* (1925): 
p 


9 No 9 9 9 
(4a) Aurbdeet er’) = — 6 Un (dr — bs!) +Ke (Mv), 
Az 2p2 e 2 2 ~p 2 2 
(46) Aut Avett p27 —3— "8? = «7a (br — 2) — KO" (Hr — ve"), 
x x (<<) 
(4c) RO PA —2k . didrzt2ke “Yip: 
Pip2 x“ x < (x?)? 1P2T oke 1¥2, 


(4d) putpstes =k oy? (go: +2) +xe°(vir +y2°). 


In equations (2) and (4) and all following ones, the subscript | refers to dif- 
ferentiation with respect to x' and 2 to differentiation with respect to x. 

In addition to these equations, the absence of a rotation term g34 vields a 
connection between the components of the four-potential 
(4e) 14 = T3 = divit doe = 0. 
It has been shown by Bonnor (1953) that in the case of an electrostatic field 
it is sufficient to solve equations (2b) and (4d). The presence of another com- 
ponent of the four-potential will not materially change the situation as we 
shall show in the following sections. 


3. THE ELECTRIC DIPOLE 
In the absence of any magnetic fields one can set ¢ = 0 and thus obtain 


the case previously considered by Curzon et al. Assuming that ¥ only depends 
on p, (2b) together with (4d) gives an equation for y, 


(5) Wop + «ec? vy." es Vp = 0, 
which has as a solution 
(6) y = A(e’*—-1), yW—O as p—0, 


where the constant A is given by 


A = V/2/x. 


*The difference between the present equations and those given by Curzon is the absence 
of terms depending on ¢ in his equations. 
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When this solution is substituted into (4d) this equation becomes 


(7) pir + pre + p2/x? = 3 (px? + pr”), 
which can be reduced to the axial harmonic 

(7a) O11 + G22 + o2/x* = 0 

by setting p = —2 Inc. 


Curzon’s solution in terms of two singularities on the axis of symmetry 
is given by 


(8) eu k=2.5, K=1 as ol, p-0O, 


where 7; and r2 are bipolar co-ordinates of a point P referred to two origins 
A and B on the x'-axis (see Fig. 1) 





Fic. 1. Bipolar co-ordinates of point singularities in cylindrical co-ordinates. 7; and rz are 
the bipolar co-ordinates of point P from the singularities, ry and @ polar co-ordinates, a the 
distance from the origin. 


ry? = (x! — a)? + (x?)?, 
(9) 9 9 2\9 
P= (x! + a)? + (x*)?, 


re 


and ¢« and e are constants of integration. Substituting (6) into the remaining 
equations (+) gives 

(10) A+ p=0, 

taking \ = 0 when p = 0, so that one finally obtains the following solution 
for two charged point masses along the axis of symmetry: 


a ve (i-g-g)"aa}. avers 


(12) A= —p = 2in 1-2-2]. 


"1 


t> ito 


It is now simple to see how from (11) and (12) an “electric dipole’ solution 
can be constructed. First of all, we shall assume that the two charges e; and 
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€, are of opposite sign but equal magnitude and thus put e; = —e = e. 
If one expands 7; and 7, in terms of the distance r (>a) and aximuthal angle 
6 (see Fig. 1), 

r;?> = r? + a? — 2arcos8@, 
r. = r? + a? + 2ar cos 8, 
and lets a go to zero and at the same time e go to infinity, so that 


p = 2ae 


remains finite, the potential is now given by 


(13) ane, 
’ | 
where 
P = /2/« b, 
while the gravitational potentials can be written as 
(14) =e? =1-— £, cos 6 + £, cos’6. 
" | 


It is seen that this solution is quite similar to the one corresponding to a charged 
mass-particle, which for a potential given by y = ¢/r has as a solution of the 
field equations (Tolman 1934) 

" Kee 4 Ke 


(15) e=e’=1-2 = oP 
2r 






2r°° 





In both cases there are two terms, one of the order of the potential and one 
of its square, the terms having opposite sign. 










THE MAGNETIC DIPOLE 
If we have only magnetic fields, then it is possible to set y = 0 and to 
consider only the solution of Maxwell's equation (2a). At first sight it would 
not seem possible to apply the method of Curzon, as the form of this equation 
differs from the one previously considered (26) and therefore also from the 
field equation (42). However, if we define a new “‘potential’’ such that 


4. 












pe e 


(16) %, = a, gd. and &, = ~3 1 





equation (2a) is simply the compatibility condition @., = y». 
From the fact that we also demand that ¢2 = dx an equation for ® is 







obtained, 
(17) 







Pip + Pe» + b2/x? = piPy + poo. 





It is seen that this is exactly of the form of (20) and hence will have as a solution 


(18) @ = A’e!* + B’ 







(assuming, of course, that @ is only a function of p again). The rest of the 





calculation goes as before and one finally obtains for the two-particle solution 
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(19) = a 


r ro 


where again m, and mz are constants of integration. 
In terms of this solution the magnetic field components are given by 


'S .2)2 (x2)? 
(20a) Fs, = ¢, = a Ve a ) ad s 
K Yr) Yo 
2 , Sle calat pte tape! |. 


tie 


Yr) le 
In order to obtain the vector-potential corresponding to a magnetic dipole 
we shall take m; = mz = m and let a go to zero and m go to infinity so that 


ll 


(206) F3. = 2 


pw = 2am 
remains finite. Expanding in terms of the quantity a/r we get for the fields 


/2 Bu (x2)? 
(21a) oi = “WV : ute) cos 8, 


iat 2 px” 3ux'x” 
(21d) do = ae” ce cos 6], 
K 


and finally for the vector-potential 
x*)?  sin’@ 
(22) @¢=M =>, 
r r 


where 
M = /2/« u. 
The appropriate solution of the field equations then again becomes 


2 


(23) e=e’=1- a cos o+5 7 cos’6. 
r 2r 


5. THE ELECTRIC AND MAGNETIC DIPOLE 

Finally we wish to consider the case when neither Y nor @ vanishes. 
The simplest possibility, although perhaps not of most physical interest, is 
to assume the presence of stationary electric and magnetic sources, viz. 
magnetic and electric dipoles. Outside the singularities, the field equations 
for empty space and the Maxwell's equation without current and charge 
distribution hold. 

From the result of the last two sections one finds that the scalar and 
“alternate’’ vector-potential must be related through 


(24) ’=yt+F, 

where ¢ and F are arbitrary constants. This will also automatically satisty 
equation (2e), which connects the two components of the four-potential. 
As a solution of Maxwell's equations, one again obtains 


(25) y = A”velP 
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where, now, the constant A” is given by 


The solution of the field equations results in 


2 
r 2 d d 2 
(26) e =e °= 1—-3cos 0+=4 cos 8, 
r 2r 
where we have defined 
d = lim 2am. 


a +0 
moo 


The corresponding scalar potential is given by 


iS D cos 0 na ea 
2 2 = an eae 
(27) y 2 ; D=d x -+e%)’ 


and the appropriate vector-potential becomes 
pyr 
M sin’6 


2 ee eee 
(28) ¢ 3 ; 


” ae. 
M= ei o- . 


It is thus seen that c measures the ratio of the electric to magnetic dipoles 
and that the total gravitational field is due to a superposition of the two. 
A simple way of seeing this is to express d in terms of Dor M and thus write 


(29) d = 4/ — D= / me ae M. 


6. CONCLUSION 


where we have defined 


Thus we have been able to find solutions of Einstein's field equations in 
empty space in the presence of magnetic and (or) electric fields. Although 
the solutions are based on rather special assumptions they describe actual 
physical situations, which may, however, play not too important a part in 
actual models. In connection with the field of the magnetic dipole it would 
be of great interest to see what would happen if it were to rotate. Unfortunately 
the line-element describing rotation is more complicated and so far no solutions 
have been obtained, although at present work is in progress along these 
lines. 
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ADIABATIC OSCILLATIONS IN LIQUID HELIUM II' 


F. D. MANCHESTER?’ AND J. B. BRowN 


ABSTRACT 

Oscillations have been observed in an experimental arrangement consisting 
of an adiabatic container placed in a helium II bath and connected with the 
bath by means of a “superleak’’. A ‘‘thermal pulse”’ introduced into the con- 
tainer caused the liquid level to oscillate. Containers of different geometries 
and employing two different types of superleak were successfully used to observe 
oscillations. With one of these the temperature dependence of the frequency 
was measured between 1.39° and 2.065° K. and found to be of the form and 
magnitude expected from the calculations of Robinson (1951). Other features 
of the oscillations have been investigated. An analogy between the oscillations 
in an adiabatic container in liquid helium II and those of a gas in a Helmholtz 
resonator is suggested and discussed. 


INTRODUCTION 

In their pioneering paper on the superfluid properties of liquid helium IT, 
Allen and Misener (1939) reported the observation of oscillations of the liquid 
level in a container which was partly immersed in a bath of the liquid and 
connected with the bath by means of a superleak, i.e., a flow channel through 
which only superfluid helium could move easily. These oscillations were 
later investigated by Atkins (1950), who used the helium II film as the super- 
leak linking the container with the bath. Atkins provided an analysis of the 
oscillation mechanism taking into account the variations in thickness of the 
helium II film and from this was able to use measurements on the oscillation 
period to obtain information on the variation of the thickness of the film 
with height and with temperature. 

Shortly after Atkins’ work, Robinson (1951) produced a general analysis 
of the oscillatory flow of liquid between a container and a bath, which showed 
that oscillations could be divided into two groups. One group, the isothermal 
oscillations, were the type already reported by Allen and Misener and by 
Atkins; the other group, which Robinson called adiabatic oscillations, had 
not hitherto been observed. The present paper deals with an experimental 
investigation of these adiabatic oscillations and gives a more detailed account 
of experiments which were briefly reported on earlier (Manchester 1955 a, 6). 
The study of these oscillations is of interest because it provides an opportunity 
for obtaining an experimental check on the equations of motion of liquid 
helium II in the region of small velocities when the equations are applied 
to the flow of liquid through a superleak. 


OUTLINE OF THE THEORY OF OSCILLATIONS IN LIQUID HELIUM II 


Isothermal oscillations have been fully discussed by Atkins (1950) and 
Robinson (1951). These oscillations occur in a system consisting of a con- 
tainer connected to a helium II bath by means of a superleak and under 
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conditions where there is very good thermal contact between container and 
bath. Their mechanism may be viewed as the periodic interchange of energy 
between the kinetic energy of the superfluid moving in the superleak and 
the potential energy of the difference in liquid levels. The type of arrange- 
ment in which these oscillations have been observed and their typical form 
are shown in Fig. I(a). Isothermal oscillations have been observed by Allen 





TIME —> 


“(b) 


Fic. 1. Representation of (a) isothermal and (b) adiabatic oscillations. The isothermal 
case shows both the film and a channel system acting as superleaks. 


ana Misener (1939), Atkins (1950), and Picus (1954) but only in the case 
of Atkins’ experiments were measures taken to ensure good isothermal 
conditions. 

The case of the adiabatic oscillations may be considered by referring to 
Fig. 1(6). Instead of giving good thermal contact between the container 
and bath as in the isothermal case, the apparatus is designed to reduce the 
thermal contact to a minimum. Under these conditions temperature differences 
can be sustained between container and bath when liquid flows from one to 
the other. For example, when liquid is flowing from the container to the 
bath, the temperature of the container liquid rises because of the mechano- 
caloric effect, setting up a temperature gradient in such a direction that the 
thermomechanical effect tends to force liquid back into the container. As 
the superfluid moving in the superleak possesses sufficient inertia, this situation 
gives rise to oscillations. 

A theoretical expression for the frequency of ideal isothermal oscillations 
was introduced by Allen and Misener and a derivation was later provided 
by Atkins for the case of the helium II film as superleak. Dash (1954) and 
Picus (1954) have produced analyses of isothermal oscillations dealing with 
particular aspects of the helium II film in the role of superleak. From the 
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analysis due to Robinson (1951) the angular frequency, w;, of the isothermal 
oscillations is given by 


(1) Yo, = (e-e) 
i p la ’ 


where p; and p are the superfluid and bulk liquid densities respectively, g 
is the acceleration due to gravity, and as is shown in Fig. 2(a), o is the cross- 
sectional area of the superleak, / its length, and a@ the cross-sectional area 


(a) 





o 


A a 
‘et 

! 
‘ 


(b) 


Fic. 2. The adiabatic container, (a), and the Helmholtz resonator, (b), showing geometrical 
parameters. 


of the observational capillary. Atkins found experimentally that the frequency 
of the oscillations decreased with increasing temperature as expected from (1). 
The frequency of ideal adiabatic oscillations w,, for the container of Fig. 2(a), 


1s 


(2) We = w,(1+a)}, 


ek ae €) 
on eve Cl pS ae 


In the above, S and C are the entropy and specific heat respectively of liquid 
helium II, (dp/d7)vap is the gradient of the vapor pressure curve, and 7 
is the absolute temperature. For the derivation of this relation the reader is 
referred to Robinson's paper. Robinson has shown that the equations of 
motion describing the flow of superfluid have solutions which can be divided 
into three regions, according to the degree of thermal linkage between the 
container and bath. These are: the adiabatic oscillatory region for zero or 


where 


eA OO NR A nw en 





486 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1957 


very low thermal linkage, an aperiodic region for intermediate values, and the 
isothermal oscillatory region for good thermal contact. The general solution 
for x, the position of the meniscus level in the container, may be written for 
an oscillatory region as 


(3) x = Ae”"* cos(pt /wa— )+Be 0", 


For the practical case of a container with some small but finite thermal 
linkage to the bath, the theoretical picture of the adiabatic oscillations is as 
follows. If the liquid in the container is suddenly displaced from its equili- 
brium position to some new position the liquid level will commence to oscillate 
with angular frequency p/w,, the oscillations being damped according to the 
magnitude of b. During this process and after the oscillations die out, the 
mean liquid level returns toward its initial position at a rate determined by 
the second term of the right-hand side of (3). For a particular temperature 
the p’s and b’s of (3) are dependent on L, the index of thermal linkage between 
container and bath (L = 0 and LZ = & correspond to ideal adiabatic and 
ideal isothermal oscillations respectively), and therefore the frequency, 
damping, and aperiodic decrease observed in a particular experiment should 
depend on the amount of thermal isolation achieved. The temperature depen- 
dence of the ideal isothermal and ideal adiabatic oscillation frequencies, for 
the particular geometry of one of the containers used in the experiments to be 
described, is shown by the smooth curves of Fig. 5. 

Equation (2) may be expressed in a form which is convenient for the 
calculation of numerical values: 


am, (w; +hwp’)}, 


ae) 


a, 


using equation (1) and the relation 
Ps Tz 
mn = 
pp C 


for the velocity of second sound (Landau 1941). It is shown in the appendix 
that wr may be identified with the natural frequency of an ideal adiabatic 
resonator. 

The method of calculation used here is very insensitive to the values chosen 
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for the entropy and specific heat, which is an advantage, because the values 
of these quantities are somewhat uncertain at the present time (Brewer, 
Edwards, and Mendelssohn 1955). Geometrical constants of the apparatus are 
given in Table I and the data used in plotting the smooth curves of Fig. 5 
are given in Table II. 








TABLE I 
GEOMETRICAL CONSTANTS OF THE APPARATUS 
Optical flats Wire-filled-tube 
apparatus apparatus 

Volume of chamber A 4.0 cm.3 6.6 cm.3 
Observational capillary diameter 1.14 mm. 1.10 mm. 
Length of superleak 0.68 cm. 11.7 cm. 
Cross section of superleak 2X10-3 cm.? 7.2X10-4 cm.? 

TABLE II 


THEORETICAL VALUES OF ai, wr, AND wa FOR THE WIRE-FILLED TUBE APPARATUS 


Wi, WR, Wa, 
PK. k radians/sec. radians/sec. radians/sec. 
1.0 0.946 2.50 0.499 2.55 
Ss 0.928 2.49 0.684 2.58 
1.2 0.923 2.47 0.952 2.63 
LS 0.920 2.45 1.27 2.74 
1.4 0.919 2.41 1.66 2.89 
1.5 0.920 2.36 2.91 3.11 
1.6 0.921 2.27 2.64 3.40 
ba 0.923 2:17 3.11 3.69 
1.8 0.925 2.05 3.51 3.95 
1.9 0.928 1.86 3.81 4.12 
2.0 0.932 1.60 3.91 4.10 
2.15 — 0.79 2.32 — 
2.19 a 0 0 0 





To obtain the calculated values of w;, w,, and wz the following experimental 
data were used: values of p,/p and p;/p from Peshkov (1946), second sound 
velocities from smoothed values quoted by Band and Meyer (1948) and by 
Maurer and Herlin (1949), the entropy of liquid helium II from Kramers 
et al. (1952), and the vapor pressure relation for liquid helium II from Bleaney 
and Simon (1939). 

EXPERIMENTAL 

Two versions of an adiabatic container were used in the experiments, one 
using the annular gap between two optically flat glass disks to provide a 
radial flow path which acted as a superleak, and later a container using a 
wire-filled tube superleak, similar to those used by Allen and Misener (1939), 
and Brown and Mendelssohn (1947). Some typical results obtained with 
the apparatus using optical flats have already been published (Manchester 
1955a). A schematic diagram of the two types of containers and also the 
essentials of the auxiliary apparatus are shown in Fig. 3. 

With both types of container the external experimental arrangements 
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Fic. 3. Schematic diagram showing the general form of the apparatus and the two types 
of experimental beaker used. A is the main volume containing the heater, B is the observational 
capillary, and C is designed to minimize any pressure fluctuations above the liquid oscillating 
in the capillary. 


were very similar. The experimental container, henceforth referred to as the 
‘beaker’, with its surrounding vacuum jacket, was enclosed in a large glass 
envelope to which liquid could be admitted from the surrounding liquid 
helium dewar. The lower end of this glass envelope was provided with a large 
copper-glass seal and a copper bottom to provide good thermal contact be- 
tween the liquid in the envelope and the liquid being pumped in the main 
helium dewar. Arrangement was made for the beaker to be raised or lowered 
so that its position could be adjusted with respect to the liquid level in the 
envelope, the “bath level’. The temperature of the liquid helium bath was 
stabilized to 10~* ° K. or better by means of a bridge system of the type de- 
scribed by Boyle and Brown (1954). The main volume A of the experimental 
beaker was fitted with a heater. 

In the case of the optical flats superleak the separation of the flats could 
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not be measured very reliably; however the general appearance of white 
light interference fringes indicated that the separation of the flats was 
approximately 10~‘ cm., although measurements on oscillations and helium II 
flow rates obtained in successive experiments indicated that the separation 
varied from one experiment to the next. Because of these variations in the 
optical flats superleak and also because a lower frequency of oscillation could 
be obtained with a longer superleak (see equations (1) and (2)), a wire-filled 
tube was used. 

The wire-filled tube was prepared in the following way. A large number of 
fine constantan wires (1200) were arranged to thread a cupronickel tube. 
The tube was then passed through successively smaller holes in a steel die 
plate to draw it down (Allen and Misener 1939). A photomicrograph taken of 
a cross section of such a tube showed that the wires were arranged in a regular 
pattern of closely packed hexagons (Allen and Misener 1939), the separation 
between parallel sides of neighboring hexagons being the channels available 
for gas or liquid flow. The mean length of the sides of the hexagons was found 
by direct measurement with a microscope to be 24.6 uw. For wires of hexagonal 
cross section, the number of parallel flow channels in the tube is three per wire 
or 3.1 for the finite number of wires in the particular tube considered here. 
With the mean channel breadth measured, the depth was determined by careful 
gas flow measurements and the application of Poiseuille’s relation for flow 
between parallel planes. The value obtained was 7.8 X 10~ cm. for the average 
depth and 7.2X10~‘ cm.? for the total area of cross section of the superleak. 
This was the value used in computing the theoretical values of we. As in the 
present case the mean free path of the atoms of the helium gas used in the 
flow measurements is of the same order of magnitude as the dimensions 
of the flow channels, Knudsen’s relation for gas flow was also used and it gave 
a value of 5X10~> cm. for the channel depth. As this value is in moderate 
agreement with that given by Poiseuille’s relation, it is possible that the 
flow was occurring in the transition region between the two types of flow. 

The wire-filled tube was fastened by means of Kovar seals along the axis 
of a glass cone which fitted into a ground-glass socket and formed the bottom 
of the beaker, as shown in Fig. 3. The cone was also evacuated to provide 
thermal insulation. Glycerin was used to seal the ground-glass joint. 


EXPERIMENTAL RESULTS 

The general pattern of oscillation behavior as established in a series of 
experiments may be summarized as follows. 

If the condition of temperature and level equilibrium between beaker and 
bath is disturbed, the meniscus in the observational capillary will move to 
some new position and start to oscillate after the position has been reached. 
The oscillations continue as the mean level of the meniscus returns toward 
the initial equilibrium position, but as they are fairly strongly damped they 
die out long before that position is reached. The general form of the oscillations 
observed at several temperatures is shown in Fig. 4. It was possible to initiate 
oscillations in several ways (Manchester 1955a) but the method used in all 
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Fic. 4. General form of the adiabatic oscillations. 


quantitative work was to discharge a condenser through the heater to release 
8.1 10 ergs in the beaker. 

Some measurements of adiabatic oscillations were obtained with the 
optical flats apparatus (Manchester 1955a) using a stroboscopic flash lamp 
method of photographing the movement of the meniscus. Such measure- 
ments served to establish that the frequency of oscillation was close to that 
predicted by theory but no useful data were obtained on the temperature 
dependence of the oscillation frequency, because of the lack of reproducible 
conditions with the optical flats superleak. 

Experiments were then carried out with a similar experimental beaker 
fitted with a wire-filled tube superleak. The two principal advantages of such 
a superleak are that its dimensions can be measured with moderate precision 
and that they remain fairly constant throughout a series of experiments. 
Because of the different geometry employed (cf. equations (1) and (2)) the 
frequency of the adiabatic oscillations was much lower with the wire-filled 
tube. As a result the motion of the meniscus could be followed with the eye, 
so the stroboscopic flash lamp and camera were replaced with a very simple 
system using a stop watch, a cathetometer, and a low wattage neon lamp. 
The observed temperature dependence of the oscillation frequency is dis- 
played in Fig. 5 and Table III. The points plotted in Fig. 5 are the experimental 
values of the oscillation frequency, the smooth curves being the theoretical 
values of w; and w, computed as described previously. Points obtained on a 
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Fic. 5. The temperature dependence of the observed oscillation frequency. 


TABLE III 


EXPERIMENTAL VALUES OF w OBTAINED WITH THE WIRE-FILLED TUBE APPARATUS IN THE 
COURSE OF FIVE DIFFERENT EXPERIMENTS 


w, 7. a, 
radians /sec. radians/sec. 


1.39 2.37 | 760 3.25 
1.432 2.58 | 825 3.39 
1.500 73 826 3.84 
1.559 906 3.65 
1.640 | 992 3.9: 
1.700 | 2.065 


particular experimental run are indicated by a common symbol. Temperatures 
were measured using the 1949 vapor pressure -— temperature scale. The 
measured periodic times are considered reliable to a little better than a tenth 
of a second. The experimental points agree in form with the temperature 
dependence of the theoretical curve; the difference between this curve and 
one fitting the points is presumably due to an uncertainty in the knowledge 
of the geometry of the apparatus. Because of some blockages which will 
occur to an unknown extent in the channels of the wire-filled tube it is not 
practical to estimate the uncertainty in the value for the superleak cross- 
sectional area and for this reason the curves have been presented in the form 
of Fig. 5. The most straightforward way to check the geometry factor would 
be to carry out isothermal oscillation measurements using the same superleak 


and observational capillary. 
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It would be desirable to obtain more experimental values in the region 
above 1.9° K. in order to strengthen the evidence for the curve having a 
maximum. However such points are difficult to obtain, as the amplitude of 
oscillation decreases and the damping increases with rising temperature. 
Also it was more difficult to obtain control over the stability of the meniscus 
when endeavoring to take readings at a temperature in the region near 2° K. 
This is because, in the region of lower superfluid velocity, the time between 
a countering change of the bath temperature and its effect on the meniscus 
was prolonged, while the sensitivity to temperature imbalance was greater 
because of the increase of the magnitude of the thermomechanical effect with 
temperature. 

Measurement of the damping of the oscillations could not be carried out 
very satisfactorily because of the superposition of the damping on the fall 
of the mean meniscus level and because of the small number of periods avail- 
able for measurement. Also minute variations in bath temperature tended 
to produce variations in amplitude which had an appreciable effect when 
diminution in amplitude was being measured. Such measurements as were 
made on the damping did not give a value for the thermal linkage between 
bath and container (the L value, see equation (3)) which agreed at all with 
that obtained from the data on the aperiodic decrease of the mean meniscus 
level. However, further examination of this question should be postponed 
until better experimental data on the damping are available. 

The rate of return of the mean meniscus level to the original equilibrium 
position was measured, that is, the aperiodic decrease represented by the 
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Fic. 6. The aperiodic fall of the meniscus after oscillations ceased. The curve for 1.992° K. 
is for a triggering pulse of 24,000 ergs. 





MANCHESTER AND BROWN: ADIABATIC OSCILLATIONS 493 


second term in equation (3). Some results typical of those obtained are shown 
in Fig. 6. The bending which appears in such curves, and which occurs in 
all to an extent not necessarily apparent in Fig. 6, is thought to be due to a 
slight warming of the bath during the period of measurements, presumably 
because the liquid level in the outer helium dewar was slowly falling. This 
warming could occur because the temperature-stabilizing bridge controlled 
the temperature of the outer dewar liquid directly and the helium in the 
glass envelope by thermal contact through the copper. 

The relaxation times for the return were determined by measuring the 
slopes of the straight line portions of the curves such as those shown in Fig. 6. 
In Fig. 7 the points are the values of these times plotted against temperature, 
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Fic. 7. The temperature dependence of the relaxation time for the fall of the meniscus after 
oscillatins ceased. 


the smooth curve is proportional to 7®. It can be shown that the relaxation 
time may be written approximately as 


_ (pa .) 
(4) r= (@)( K : 


where A is an effective coefficient of thermal conduction between the beaker 
and bath. If it is assumed that the major contribution to the thermal con- 
duction is through the superleak (see below) then A will vary as 7° for a con- 
stant temperature difference as found experimentally by Mever and Mellink 
(1947), for a superleak of 0.5 4 width. Using the approximation that the 
entropy of liquid helium II is proportional to 7°* over the temperature 
range of interest, the relaxation time is given as proportional to 7*. A more 
direct comparison between theory and experiment cannot be made because 
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there is no satisfactory quantitative treatment for heat transport through 
superleaks of the dimensions used in these experiments. 

Calculations were made to determine the magnitude of the heat flux into 
the container along possible paths such as the electrical leads and the metal 
of the wire-filled tube. None of these seemed to provide a thermal link of the 
magnitude required for agreement with the experimental observations. 
The expected heat flux through the wire-filled tube superleak calculated by 
means of the relation of London and Zilsel (1948) was too small by a factor 
of 10-. This at least is in the same direction and of about the same magnitude 
as the discrepancies noted by London and Zilsel for superleaks of the dimen- 
sions considered here. 

The temperature dependence of rates of fluid flow occurring during the 
oscillations is shown in Fig. 8. These are simply the average velocities between 
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Fic. 8. The temperature dependence of flow rates. Upper curve: the flow rate up to the 
oscillation point. Lower curve: the flow rate during the first backswing of the oscillation. 
turning points. Both of the curves have a form similar to that usually associated 
with the temperature variation of critical velocities. However, in the present 
case the velocities are very low; that for the superfluid velocity at 1.39° K. 
is indicated on the graph and is much lower than published values of critical 
velocities (see Winkel, Delsing, and Poll 1955). At the time of writing no 
opportunity has occurred for measuring critical velocities in steady flow 
experiments with the same wire-filled tube. It was suggested by London 
(1951) that the study of adiabatic oscillations might yield information on 
conditions under which there was a critical flow velocity. The present experi- 
ments do not provide such information but it is possible that useful data on 
flow rates may be obtained from a modified experiment to be discussed 


below. 
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The equations for the adiabatic oscillations are solvable for the temperature 
difference as well as the liquid level difference (see Robinson’s paper) although 
from a knowledge of the magnitude of the fountain effect and the amplitude 
of the liquid level oscillations, the expected temperature oscillations would 
be of very small amplitude, 10-°° K. or less. In spite of this small amplitude, 
it was decided to try to detect the temperature oscillations because it was 
thought possible that their recording might have been more convenient. 
No temperature oscillations were observed in several attempts, and effort 
in this direction was then abandoned in favor of concentrating on measure- 
ment of the liquid level oscillations. The limitations on the temperature 
measurements were the extreme sensitivity required of the thermometers 
together with the necessity for a fast response. 

DISCUSSION 

The general behavior of the oscillations agrees with the theoretical picture 
presented by Robinson (1951) while the measured frequency in the range 
1.39° K. to 2.065° K. has the same form of temperature dependence and the 
experimental values are fairly close to those computed theoretically for 
the apparatus. 

The observed agreement with the theoretically predicted values of fre- 
quency gives support to the assumption made in the theory that the fluid 
composition in the superleak is that of the bulk liquid and that it is the super- 
fluid fraction which moves when flow occurs. As the cross section of the 
superleak is essentially independent of temperature over the range considered, 
the evidence given here supports the assumption made by Atkins (1950) 
concerning the substitution of Peshkov’s p,/p values (Peshkov 1946) in the 
isothermal oscillation equations. 

If the geometry factor could be determined more precisely a better experi- 
mental check could be made on the theoretically expected values for the 
damping, through the quantity p/w, (a/we of Robinson’s paper). This might 
be more easily achieved by using a different approach suggested by the 
resonator analogue discussed in the appendix. This involves the use of a 
source of second sound placed in the liquid outside the adiabatic container 
and designed to operate at frequencies in the region of its resonance. Tem- 
perature oscillations generated by this source should produce movement of 
the meniscus in the capillary of the adiabatic container and enable the reso- 
nance curve of the container to be traced by visual observation of the amplitude 
of the meniscus motion. Measurements of the damping carried out in this 
manner should give much more information than was obtained in the present 
experiments. 
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APPENDIX. THE RESONATOR ANALOGUE 
Equation (2) may be simplified by making the approximations that 
a = aTS*/gVC 
and that @ is much larger than unity. Then 


Nps oO TS" , 


Wy = a 
tell € 
and making use of the relation for #2, the second sound velocity (Landau 
1941), 


_ TS 
"le 


Ue 
the expression for w, becomes 
(I) 


Equation (I) gives the frequency for ideal adiabatic oscillations assuming 
that w, > w; and neglecting the effect of the vapor pressure over the liquid 
in the container. For reasons which will be given below the expression for 
w, in the form of equation (1) will be denoted by wz. It is possible to derive 
the relation for we from an exact treatment of a simplified case which takes 
only thermal effects into account and neglects gravitation. For higher tem- 
peratures in the liquid helium II region, where the thermomechanical effect 
is relatively very large, the temperature dependence of wz is close to that of 
w, as given by equation (2), but with decreasing temperature the value of 
wr falls to zero at 0° K. whereas w, tends to w;, the isothermal value. 

It is interesting to compare equation (I) with the expression given by 
standard acoustical theory for the resonant frequency of a Helmholtz resonator 
in a gaseous medium. For such a resonator, in the region of long wavelengths, 
the frequency of resonance is given by (see for instance Stewart and Lindsay 
1930) 


wy = c(A/I'V)}, 


where c is the velocity of sound in the gas, 4 the cross-sectional area of the 
resonator neck or opening, /’ the effective length of the neck, and V’ the 
volume of the resonator, Fig. 2(b). There is good reason, therefore, for regard- 
ing the ideal adiabatic container in the second sound field as analogous to the 
Helmholtz resonator in the first sound field, both being considered for the 
long wavelength case only. 

Considered from this point of view of association with second sound, the 
adiabatic container may be thought of as either a resonator or an oscillator. 
As an oscillator emitting radiation it could be regarded as a source similar to 
that used by Peshkov (1948), who placed a superleak in front of a mechanical 
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transducer emitting first sound so that relative motion of super and normal 
fluid was produced and the device acted as a source of second sound. 

It should be noted that the approximation w, > w; used above is affected 
by the geometry of the adiabatic container. In the experiments described 
in this paper the geometry of the container was chosen to give a low value of 
w_ so that experimental observation would be easier. Thus in this particular 
case the approximation w, > w; is not a good one but experimental conditions 
could easily be chosen to make it much better. 
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NOTES 








ON THE THEORY OF MULTIPLE MESON PRODUCTION! 
YASUSHI TAKAHASHI? 


In a previous paper (Takahashi 1956a), a general expression for the S- 
matrix for multiple meson production was derived, and the relation between 
Fermi’s thermodynamic approach and the field theoretical treatment dis- 
cussed, by determining the most probable value of the meson production 
probability, neglecting meson correlation. 

On the other hand, it has been pointed out that the Fermi model is too 
simple to explain experimental data available so far (Bethe and de Hoffman 
1955; Landau 1951). This point has been discussed in detail by Landau and a 
Japanese group (Koba 1956), who stress the importance of the meson—meson 
interaction. 

In this note, the meson—meson correlation is discussed starting from the 
general expression of the S-matrix (Takahashi 1956b). According to the pre- 
vious paper (Takahashi 1956)), the probability for the ~-meson production 
is generally given as 


1 = dk 2 
o = ideas JU aay Mote BOF, 


i 


ith,...2 = 2 es BS, 


(permutation) 


ott ste) = fo. PTT (Geo exp(—aterd} 
x(®) , T*(O(x1), .. . , Oxy) ) Bo). 


Since the function |¢,(fi,...,f,)|? in (2) is a positive symmetric function 
with respect to the f’s, it can be expanded in the following way: 


(3) | dp (fi, Seats f,)|° a TTou(t.) IT oa(ts t,) IT oa(ts e;, 9) Peg 


(2)* 


where (7, 7) signifies the product over all possible combinations of ¢ and j 
and (i, j, k) the product over all possible combinations of 7, j7, k, and so on. 
The expansion of a symmetric function in the form (3) is always possible but 
is not always unique. This will be discussed in a detailed paper later. We 
suppose, at any rate, that the go(f;, £;), os(fi, f;, f),... are determined as 
close to unity as possible, and we will neglect, for the sake of simplicity, the 
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terms which include more than two arguments. If we express (3) in terms 
of the meson occupation number, we get 


(4) ldn(fa,- + f,)|° = nl TT {or(to))" Tf etn fox) pire) 


x I] I] { bo (Eiry, Bea) } 2", 


where 7; is the number of mesons with the momentum fi, 22 the number of 
mesons with the momentum f,), and so on. 
The total probability is 


(5) Zz dW, (m1, mo, « 


where 
(6) dW,(m, m2,... atten)" ‘TI { be (Fir, tn er* 


X ITI {¢2(%, fo) }"™ 


8 r#8 


1 
Q(F¢s)) — Kis $1(f(s)). 


We can easily determine the most probable dW,(m, m2,...), subject to the 
energy conservation law, by taking a variation with respect to n, of the 
quantity 
(8) I = log dW,(m, m2, ...) + B(AE—DOKis) ms). 
Thus, 
(9) o=-> oF in, “ei én, log n,+ x én, log fe 
Tr r 
—B YE Kinin t+ DD in, . n, log{ $2(tin, fo) }.- 


We have, therefore, 


Is a 
(10) ny = Af) exp(— 6K») TT] { b2(Fen» fear) } 
= FOF) exp(— 8K) I] exp (Ars ns) 


1 = = 
ss Q(f..)) exp(— BK y)) exp(A, 2), 


where we put 


$2(Frr), f(s) = exp Arsy 


(11) Dw / DS te 
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The last factor’of (10) takes care of the effect coming from the meson—meson 
interaction and effectively modifies the ‘spectrum volume’’ Q(f). This is 
encouraging in view of the fact that raising the volume @ is suggested 
to bring the ratio of two-meson to one-mescn events into agreement with 
experiment. 

The relation between equation (10) and Landau’s theory will be discussed 
in the next paper. 

The author would like to thank Dr. G. Herzberg and Dr. T. Y. Wu for their 
kind hospitality. He is also indebted to Dr. G. R. Allcock for his stimulating 
discussion. 


BetHE, H. A. and pE HorrMAN, F. 1955. Mesons and fields, Vol. II (Row, Peterson & 
Co., Evanston, IIl.). 
Kosa, Z. 1956. Private communication. See also the Proceedings of the Rochester Con- 
ference, 1956. 
Lanpau, L. D. 1951. Izvest. Akad. Nauk. S.S.S.R. 17, 51. 
TAKAHASHI, Y. 1956a. Can. J. Phys. 34, 378. 
19566. Nuovo cimento (10), 4, 531. 


RECEIVED AuGustT 23, 1956. 
DIVISION OF PuRE Puysics, 
NATIONAL RESEARCH COUNCIL, 
Ottawa, CANADA. 


A SIMPLE FORMULA FOR EVALUATING THE MADELUNG CONSTANT OF AN 
NaCl-TYPE CRYSTAL 


J. K. MACKENZIE 


In a recent paper, Benson (1956) gave a new formula for the Madelung 
constant of a crystal of the NaCl type. The purpose of this note is to give a 
direct proof of this formula. 

The Madelung constant ay is given by the triple summation 


(1) —ay = dow (—1) tht 8+’ +h)", 


11, lo, lg 


where the prime on the summation means that the term with /; = /; = /; = 
0 is omitted. By means of a physical model with repulsive forces, Benson 
showed that 


(2) ay = 48x01 vp exp(— xp’) /[1+exp(— xp')]’, 


where v, is the number of ways p = /;?+/.2 can be represented as the sum of 
squares of positive odd integers. He also showed that equation (2) gives the 
correct value of ay to at least nine significant figures.* 


*The last two digits of the first entry in the second column of his Table I should read 63 
instead of 54. With this change 10 significant figures are obtained correctly. 





It is clear that 


(3) -ay =3 D%, (-1)" 88 +h +h)”, 


11, lo, lg 


, eo 
(4) =6r! ee (agp 1" uw exp[— (dy +le +15" )uldu, 

1, lo, Ug 0 
on using the integral which defines the gamma function. Interchange of the 
order of the summation and the integration can be justified and further, since 
all terms in equation (4) with /; = 0 are zero, the summation can be extended 
over all values of /;, /2, 73. Thus 


(5) —ay = 127? >) (—1)" 1 exp(—1,°u)f(u) u'du, 
ty 0 


where 
f(u) y nisl (—1)"*" exp[— (le?+15")u], 


lg 


le 
(1/u) > da oo exp[— (ly?+E15") 3” /4u], 

lo, lgodd 
on using a known result (Whittaker and Watson 1935). When the series for 
f(u) given by equation (7) is substituted into equation (5), interchange of 
the order of summation and integration can again be justified and the integrals 
can be evaluated by using the result (see Watson 1944) 


(8) J u"' exp[—(a’u+b6"/u)]du = 2(b/a)" K,(2ab). 
0 


Finally, when A,(x) is put equal to (#/2x)'/e-* and the summation over /; 
performed, the result given in equation (2) follows. 

Similar methods can be used to evaluate other lattice sums but the resulting 
series do not always converge more rapidly than the series which are already 
known. This is so, for example, in the case of the series obtained by the present 
method for the Madelung constant for a CsCl-type lattice. 
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THERMAL CONTACT IN A PARAMAGNETIC SALT! 
J.S. DuGcpace, D. K. C. MACDONALD, AND A. A. M. CROXxON 


In recent measurements of the electrical resistance of metals below 1° K. 
(Dugdale and MacDonald 1957) we wished to measure the thermal contact 
resistance between a gold wire and a paramagnetic salt which had been 
crystallized onto it. To do this, an experiment was made with two “pills” 
of the salt (iron ammonium alum) each of which was crystallized onto a speci- 
men of the gold wire. The two gold wires were joined by a copper wire whose 
thermal resistance was estimated to be considerably less than the thermal 
contact resistance to be measured. In one of the pills a rhodium-—platinum 
resistance element was embedded to act as an electrical heater. Both salts 
were then demagnetized from the same field at about 1° K. in the normal way; 
after demagnetization, the temperatures of the two pills were measured (by 
their magnetic susceptibilities) at regular intervals. From time to time 
known quantities of heat were supplied to one salt by means of the heater, 
and from the temperature-time graphs of the two salts a value of the thermal 
contact resistance was deduced. 

However, a peculiar phenomenon was also observed. As soon as the heater 
was switched on in the “hot” pill, the temperature of the ‘‘cold” pill was 
found to fall and it was not until later that the flux of heat from the “‘hot”’ 
pill caused the temperature of the cold pill to start to rise in the expected 
manner. It was suggested that an adsorbed helium film showing super- 
fluidity might be responsible for this behavior,* and we were then concerned 
whether the presence of such a film might also play a significant role in the 
observed normal thermal conduction. If this were so, our estimate of the 
thermal contact resistance between wire and paramagnetic salt might be 
seriously in error. We have, therefore, carried out further experiments with a 
double pill arrangement to examine this question. 

We assumed that any heat flow through the helium film would depend 
essentially on the perimeter of the link joining the two pills. We therefore made 
two experiments with the same pills and gold wires but with different links. 
These had the same length and over-all perimeter but quite different thermal 
resistances. In the first experiment connection between the pills was made 
by a heavy copper wire of 0.1 in. diameter (some of the results are shown in 
Fig. 1(a)). In the second experiment the link was made of manganin. It con- 
sisted of 20 manganin wires of the same length as the copper and of 0.005 in. 
diameter giving thus the same total perimeter; its thermal resistance however 
was estimated to be at least a hundred times bigger than that of the copper. 
The results with the manganin link are shown in Fig. 1(6). It is clear that in 
this case the two pills remain effectively isolated even when the temperature 
of the ‘‘hot”’ pill has risen as high as 0.5° K. 

From these experiments we deduce that with the copper connecting rod 
the total thermal conductance is about 50 ergs/sec.deg., and with the manganin 


1Issued as N.R.C. No. 4237. 
*We are grateful to Dr. F. D. Manchester for making this suggestion. 
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Fic. 1. (a2) Temperature-time curves for pills connected by heavy copper wire. 

(b) Temperature-time curves for pills connected with manganin wires. 
connection only about 3 ergs/sec. deg. in the comparable temperature range. 
Therefore, even if we attribute all of the latter heat flow to the helium film 
it is still small compared with that when the copper link is used. Consequently 
it appears that the thermal conductance of the contact between each salt pill 
and the gold wire is of the order of 100 ergs/sec. deg. at these temperatures, 
which is even better than we thought originally. 

It may be remarked that in the present experiments the anomalous cooling 
of the ‘‘cold”’ pill is much less noticeable than in the earlier work. We believe 
that this may be due to the achievement of a better vacuum in the present 
experiments, so reducing the formation and/or superfluid flow of the helium 


film. 
DvuGpDALE, J. S. and MacDonaLp, D. K. C. 1957. Can. J. Phys. 35, 271. 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words an, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Structure de la molécule CH.-CH—C=N par diffraction des électrons* 


Nous avons obtenu, au moyen du diffracteur électronique de I’Université de Montréal 
(Grégoire et Rouault 1951) des photographies de diffraction d’électrons de 50 kv. par les 
vapeurs du cyanure de vinyle CH>=CH—C=N. Un faisceau sensiblement paralléle d’électrons 
rencontre un faisceau moléculaire de gaz. Les électrons diffractés impressionnent une émulsion 
photographique placée dans le plan focal d’une lentille magnétique. La répartition de I’intensité 
I des électrons diffractés en fonction de s = 42(sin 6/2) /A (@ = angle de diffraction, \ = longueur 
d’onde associée aux électrons) a été tracée au moyen d’un microphotométre Leeds and North- 
rup, en faisant tourner la photographie autour du centre des anneaux de diffraction pendant 
l’enregistrement (Bélair et Rouault 1951). La structure de la molécule a été étudiée par la 
méthode de comparaison des courbes expérimentale et théorique représentant la fonction 
s3J (Rouault 1940). 


La structure des composés a double liaisons conjuguées et des composés de C=N et l'étude 
du spectre hertzien (Wilcox, Goldstein, et Simons 1954) du cyanure de vinyle rendaient 
trés probable la structure plane de la fig. 1 avec les caractéristiques suivantes: 


C=C = 1.38A C—C = 146A C=N = 116A 
C—H =1.09A ZC=C—C = 120° ZH—C—C = 120° 


La fig. 2 reproduit la courbe théorique correspondant a ce modéle et la courbe expérimentale 
déduite de l'étude de la photographie. L’intensité J des électrons diffractés par un gaz con- 
stitue de ces molécules et donnée par la formule: 
4 sin slj; 
Sf] = 2iXjbi9; ee oe 


She; 


gi(s) = Zi =e) 


sh : 


F;(s) = facteur de diffusion pour rayons X de l’atome /, 
l;; = distance des 2 atomes numéros i et 7. 


La courbe II de la fig. 2 représente cette fonction s3J. La courbe expérimentale (courbe III) 
a une forme générale un peu différente due trés probablement a une erreur dans la détermination 
de la courbe de sensibilité de l’émulsion utilisée; mais I’influence sur les abscisses des maxima 
et minima de s°J est faible. Si la forme choisie pour la molécule (rapport des différentes 
distances interatomiques) est bonne, on peut calculer une distance particuliére /, C==C par 
exemple, en écrivant pour chaque point caractéristique, maximum ou minimum, Sexp lexp = 
Sth/th; les résultats de ce calcul sont présentés dans le tableau I; les différentes valeurs de ce 
tableau sont en bon accord et ne montrent pas de dérive systématique en fonction de s dans 
une méme photographie. On peut en conclure (Rouault 1940) que le modéle choisi est en 
accord avec nos mesures. La moyenne des nombres du tableau I est 1.38 A, confirmant non 
seulement la forme de la molécule, mais aussi sa grandeur. 

Le résultat ci-dessus est confirmé par le calcul de la fonction s*J pour un modéle pour lequel 
les distances des liaisons sont inchangées sauf la distance C=C qui a été prise égale a 1.35 A; 


*Subventionné par le Conseil National de Recherches du Canada. 
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LETTERS 


CyH3N 


LT —— ~— COURBE THEORIQUE (intensite atomique) 


_e— )=— COURBE THEORIQUE (intensité totaie) 


ee evewtecewee POINTS EXPERIMENT AUK 


TABLEAU I 
Distance C=C 


10 juin, série 
10 juin, série no 4 10 juin, série no 2 no 3 


Plaque Plaque Plaque Plaque Plaque Plaque Plaque Plaque Plaque Plaque 
nol no 2 no 3 no 4 no 1 no 2 no 3 no 4 no 1 no 2 


1: 1.378 1.375 1.378 1.384 1.375 1.370 1.380 1.384 1.370 
1.376 1.380 1.380 1.372 1.383 < 1.388 1.380 1.376 1.384 

Ma 376 1.378 1.376 1.376 1.387 1.378 1.382 1.376 

is 372. 1.375 375 

M 37: 375 378 (1.391) 1.374 

he 382 380 «1.375 

VV, 373 .e 386 1.382 


Moy. 1.374 1.378 1.377 37? 379 ‘ 380 1.378 ‘ 377 


Moyenne générale: C=C — 1.38 A. 


le tableau des distances C==C calculées comme ci-dessus non seulement montre des résultats 
deux fois plus dispersés, mais une augmentation systématique des valeurs de C=C déduite 
des maxima successifs d'une méme photographie lorsqu’on se déplace vers les s croissants 
(fig. 3). 

L’écart quadratique moyen déduit du méme tableau I est ¥ 3? ~ 5X10™, soit en valeur 
relative 4X 107-3. Mais a cause de I'incertitude sur la forme de la molécule, incertitude difficile 
a chiffrer, on ne peut pas dire que la distance C=C soit connue a mieux que 0.01 A pres. 
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A {distances 
c=C_ 
1.39] Colculées 






modele 


C=C 130% 
137. i 
3 modeie 
C=C 135k 





578 T7I2 734 885 948 12.04 1272 


Fic. 3. 
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